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Theoretische und experimentelle Untersuchungen 
an Tandem-Pumpengittern starker Umlenkung * 


Von Hideo Ohashi 


E 1. Einleitung und Aufgabenstellung. Beim Entwurf eines Axialverdichters kann manchmal 
der Fall eintreten, daf die Umlenkung in einer Schaufelreihe so gro8 sein muB, daB ein einreihiges 
= ritter diese starke Umlenkung nur mit sehr groBen Verlusten verwirklichen kann, z. B. beim Nach- 
s eitrad eines Hochdruckaxialverdichters oder beim Leitrad eines Uberschallaxialverdichters. In 


GitterI 


Fs 
i 
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b c 


Abb. 1, Bezeichnungen am Schaufelgitter. 
a) Geschwindigkeitsdreieck; b) Anordnung als einreihiges Gitter; c) Anordnung als Tandem-Pumpengitter. 


iz 


E solchen Fallen kann man bekanntlich statt einer Schaufelreihe ein aus zwei Schaufelreihen bestehen- 
) des Tandem-Pumpengitter benutzen. Durch diese Aufteilung der gesamten Umlenkung auf die 
| beiden Gitter wird die Belastung der einzelnen Gitter etwa halbiert, und dadurch kénnen die Ver- 
uste trotz des komplizierteren Strémungsvorganges erheblich vermindert werden (vgl. Abb. 1). 
Bei der Anwendung des Tandem-Pumpengitters ist es jedoch bis heute tblich gewesen, das Tandem- 
Gittersystem einfach als zwei selbstandige Gitter zu betrachten und die Ergebnisse fir Einzel- 
chaufelgitter auf die beiden Gitter einzeln anzuwenden ohne Beriicksichtigung der gegenseitigen 
Interferenz. Diese vereinfachte Behandlung kann bei gewissen Tandem-Anordnungen zu groBen 


; * Gekiirzte Fassung der von der Fakultat fiir Maschinenwesen der Technischen Hochschule Braunschweig 
) genehmigten gleichnamigen Dissertation. Berichter: Prof. Dr. H. Schlichting; Mitberichter: Prof. Dr. H. Peter- 
mann. Die Untersuchungen wurden im Institut fiir Stroémungsmechanik der Technischen Hochschule 
| Braunschweig unter der Leitung von Prof. Dr. H. Schlichting durchgefiihrt und durch die Deutsche Forschungs- 
emeinschaft und ein Stipendium des Verfassers vom Deutschen Akademischen Austauschdienst erméglicht. 


— Bei der Abfassung dieser gekiirzten Fassung hat Herr Dr.-Ing. K. Gersten mitgewirkt. 
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Abweichungen von der wirklichen Strémung fiihren. AuBerdem muB man dabei auf die Még- 
lichkeit verzichten, durch Ausnutzen des Interferenzeinflusses eventuell sogar geringere Verluste 
zu erreichen als ohne Interferenz. 

Die vorliegende Arbeit soll die Frage beantworten, mit welcher Anordnung ein Tandem-Pumpen- 
gitter eine gegebene Umlenkung mit geringstem Verlust verwirklichen kann, und zwar unter Be- 
riicksichtigung der gegenseitigen Interferenz. Eine bestimmte Tandem-Gitteranordnung laBt sich 
durch acht geometrische Parameter kennzeichnen, namlich durch die Profilformen, die Schaufel- 
winkel und die Teilungsverhaltnisse der beiden Gitter und durch den Abstand und die Versetzung 
zwischen den beiden Gittern. 

Um eine optimale Tandem-Gitteranordnung fiir ein gegebenes Geschwindigkeitsdreieck zu be- 
stimmen, mu8 man zunichst den Einflu8 der acht Parameter auf die Verluste einzeln untersuchen 
und danach die zur optimalen Anordnung gehérenden acht Parameter auswahlen. 

Die Untersuchungen sind hierbei zweckmaBig in die beiden folgenden Teile aufzuspalten: — 
a) Untersuchungen iiber die optimalen Tandem-Gitteranordnungen mit unendlich groBem Abstand 
zwischen den beiden Gittern (Anderung der sechs Parameter fiir die Geometrie der beiden Gitter), 
b) Untersuchungen iiber die Anderung der Gittereigenschaften mit dem Abstand und der Ver- 
setzung bei festen Anordnungen der beiden Gitter (Anderung der zwei Parameter fiir die gegen- 
seitige Stellung der beiden Gitter). 

Der erste Teil wird hierbei auf Grund der schon vorhandenen systematischen Ergebnisse fiir 
ebene Schaufelgitter mit Profilen NACA 0010, 4410 und 8410 behandelt!. Es ergeben sich daraus _ 
die optimalen Tandem-Gitteranordnungen fiir ein beliebig gegebenes Geschwindigkeitsdreieck ohne | 
Beriicksichtigung der gegenseitigen Interferenz. 

Im zweiten Teil, der sich mit der aerodynamischen Interferenz zwischen den beiden Gittern | 
befaBt, werden die Strémungserscheinungen an Tandem-Pumpengittern bei verschiedenen Stellun- 
gen der beiden Gitter zueinander auf experimentellem und theoretischem Weg untersucht. Auf 
Grund der Ergebnisse der beiden Teile kann das in der Praxis vorliegende Entwurfsproblem eines 
optimalen Tandem-Pumpengitters in einfacher Weise gelést werden. 

' Erganzend wird der InterferenzeinfluB der Dellen vom ersten Gitter auf die Verluste des zweiten | 
Gitters, welcher in Tandem-Pumpengittern sehr ausgepragt ist, jedoch mit der Potentialtheorie | 
nicht erfaBt werden kann, in vereinfachter Weise experimentell untersucht. 

Uber den Interferenzeinflu8 des Tandem-Schaufelgitters liegen bereits einige experimentelle | 
und theoretische Arbeiten vor. Experimentell hat K. Fickert? ein Pumpengitter aus Spalt- 
fliigeln untersucht, und zwar mit verschiedener Spaltgeometrie. H. Marcinowski* hat ein Tandem- 
Pumpengitter mit verschiedenen Abstanden und Versetzungen gemessen und einige charakteri- 
stische Eigenschaften des Tandem-Pumpengitters herausgestellt. AuBerdem hat H. E. Sheets* ein 
Pumpengitter aus Spaltfliigeln in einem einstufigen Axialverdichter verwandt und damit einen 
hohen Wirkungsgrad und eine groBe Férderhohe erreicht. 

Den rein potentialtheoretischen InterferenzeinfluB haben W. E. Spraglin® und S. Abe® mit 
Hilfe der konformen Abbildung sowie S. Miiftiioglu’?’ mit der Singularitatenmethode untersucht. 
Natiirlich kann dabei der in der wirklichen Strémung sehr stark auftretende ReibungseinfluB nicht 
mit erfaBt werden. 


: 
: 


2. Bezeichnungen (vgl. Abb. 1, 11, 20 und 26). Es sollen bedeuten 


l die Schaufeltiefe, 

t die Gitterteilung, 

sn die maximale Schaufeldicke, 

z=x+1y die komplexe Koordinate (x in Richtung Schaufelsehne, y senkrecht dazu), 
rund s krummlinige, orthogonale Koordinaten senkrecht und langs der Skelettlinie, 


* L. Speidel und N. Scholz, Untersuchungen uber die Strémungsverluste in ebenen Schaufelgittern, VDI- 
Forschungsheft 464 (1957). 

2 K. Fickert, Forsch. Ing.-Wesen. 16 (1949/50) S. 141. | 

° H. Marcinowski, Messungen an axialen Schaufelgittern fiir verzégerte Strémung, Bericht Nr. VL 13105, 
J. M. Voith (1946). Vgl. auch ,,Das Versuchswesen der Maschinenfabrik J. M. Voith‘ S. 24 (1949). 

4 H, E. Sheets, Transactions ASME, 78 (1956) S. 1683. 

° W. E. Spraglin, Flow through cascades in tandem, NACA TN 2393 (1951). 

®° S. Abe, On the mutual interference of neighbouring two cascades composed of hydrofoil or aerofoil 
profils, Reports of the Institute of High Speed Mechanics, Téhoku Univ., Vol. 7 (1956). 


s erence Strémungsuntersuchungen an Tandem-Schaufelgittern, Jahrbuch 1957 der WGL., 
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fy. die Ordinaten der Skelettlinie, 
lx die Dickenordinaten der Schaufelkontur, 
i a den Abstand zwischen der Hinterkante des ersten Gitters und der Vorderkante 
des zweiten Gitters, senkrecht zur Gitterfront, 
h die Versetzung der Vorderkante des zweiten Gitters von der Bahnlinie der Nach- 
. laufdelle des ersten Gitters, parallel zur Gitterfront, 
3% die komplexe Koordinate des zweiten Nullpunktes beziiglich des ersten Koordi- 
natensystems, 
8b) und A, die Versetzung der Nullpunkte der beiden Koordinatensysteme senkrecht und 
o parallel zur Gitterfront, 
'D Bs den Schaufelwinkel, 
i ; und fp, den Zu- und Abstrémwinkel, 
'& Pz den Zwischenstrémwinkel, 
‘B46 die Umlenkung, 
A den Staffelungswinkel, 
B die charakteristische Breite einer Delle, 
‘WY W, und W, die Geschwindigkeit weit vor und weit hinter dem Tandem-Pumpengitter, 
W, die Geschwindigkeit zwischen den beiden Gittern bei sehr groBem Abstand, 
eV oo die Grundstrémungsgeschwindigkeit (vektorielle Mittelgeschwindigkeit von Zu- 
und Abstrémung des einzelnen Gitters), 
U; und Vp die Geschwindigkeitskomponenten der Translationsstrémung, 
u und v die Geschwindigkeitskomponenten der durch Singularitaten induzierten Stré- 
mungsgeschwindigkeit, 
w=u—iv die komplexe induzierte Strémungsgeschwindigkeit, 


Wy, W, und Wz die Strémungsgeschwindigkeit in MeBebene A, V und Z, 
g den Gesamtdruck, 

Pp den statischen Druck, 

Q die Quellstarke, 

a]. die Zirkulation, 

q die Quellbelegung, 

y die Wirbelbelegung, 


i = Ag/ S W; den Verlustbeiwert, 


Re = WI/y die Reynoldszahl, 
v den Energieverlustbeiwert [Definition nach Gleichung (6.1)], 
K = F/F, die Uberdeckungsziffer (vgl. Abb. 26). 


Als Indizes sollen in folgender Bedeutung benutzt werden: 


zur optimalen Tandem-Gitteranordnung gehdorig, 

0 zur Strémung ohne gegenseitige Interferenz gehérig, 
j=Iund II _— zum ersten und zweiten Gitter gehérig, 

u, a fiir Werte in Richtung parallel baw. senkrecht zur Gitterfront, 
fir Quellbelegung, 

fir Wirbelbelegung, 

fiir Werte auf der Schaufelkontur, 

fiir Werte auf der Skelettlinie, 

fiir Werte der Einzelschaufel, 

fir Werte des Einzelgitters, 

fiir Werte des Restgitters, d. h. fiir Werte aller iibrigen Schaufeln des Einzel- 
gitters auBer der betrachteten Kinzelschaufel, 

fiir Werte von der anderen Schaufelreihe. 


YW ram AR 


3. Optimale Gitteranordnungen des Tandem-Pumpengitters ohne Beriicksichtigung der gegen- 
seitigen Interferenz. Gesucht sind die optimalen Tandem-Gitteranordnungen, die ein gegebenes 
Geschwindigkeitsdreieck mit geringstem Verlust verwirklichen, wobei der Abstand zwischen den 
beiden Gittern als unendlich groB angesehen wird und infolgedessen die aerodynamische Inter- 
ferenz nicht in Betracht gezogen ist. 

Zur Bestimmung der optimalen Tandem-Gitteranordnungen tritt zunachst die Frage auf, wie 
die gesamte Umlenkung auf die beiden Gitter verteilt werden soll, um den geringsten Gesamt- 


14* 
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verlust zu erreichen. Die Aufteilung der Umlenkung wird hierbei durch das Belastungsverhaltnis 
AW.4:/AW,, gekennzeichnet (vgl. Abb. 2). 
Wenn das Belastungsverhaltnis einmal festgelegt wird, reduziert sich das Problem, die optimale 
Tandem-Gitteranordnung zu bestimmen, im wesentlichen darauf, zwel optimale einreihige Gitter- 
anordnungen zu bestimmen, wobei die Umlenkung Af; = 6, —B in Gitter I, die Umlenkung 
ABir = By — B; in Gitter IT einzeln verwirklicht werden soll. 


+ =149° 
Curteyn Fy =160° é bur bn iy 7 | 


O# 


AW,7/AW, 005 


baCim aa TE Bee one) (4 
Curly fy=120°U. 110° 
Sa o . 005) 
Pee eee | S ee — == 
a7 aN ? cy Cea) Ea tf 
(z+ Syz)min AW, 7/4, . 
a ae ee et 


08 OH 05 G6 Q7 € 03 OY 05 06 Q7 


Abb. 2. Anderung der Gesamtverlustbeiwerte des Tandem-Pumpengitters mit dem Belastungsverhiltnis bei festen Zu- und 
Abstrémwinkeln (ohne Beriicksichtigung der gegenseitigen Interferenz). Benutzte Profile: NACA 0010, 4410 und 8410. 


Gesamtverlustbeiwerte oy I + oy rae 
2 
--------- Optimales Belastungsverhaltnis beim Minimalverlust. 


Das Problem, fiir ein gegebenes Geschwindigkeitsdreieck eine optimale einreihige Gitteranord- | 
nung zu ermitteln, ist bereits verschiedentlich behandelt worden! ?:3. Um das optimale Belastungs- | 
verhaltnis (AW,1/MW,,)* zu bestimmen, bei dem das Tandem-Pumpengitter das gegebene Ge- | 
schwindigkeitsdreieck mit dem geringsten Verlust verwirklichen kann, wurden hierbei zunachst— 
zwei fiir sich allein optimale einreihige Gitteranordnungen bei verschiedenem Belastungsverhaltnis- 
bei festem Zu- und Abstrémwinkel ermittelt. Die zu den optimalen Gitteranordnungen gehérenden 
Gesamtverlustbeiwerte Cy,; + Cy, wurden dann in Abb. 2a bis 2e iiber dem Belastungsverhaltnis 
AW,,;/AW,, aufgetragen. Hierbei liegen die Ergebnisse der ebenen Schaufelgitter von Profilen 
NACA 0010, 4410 und 8410 zugrunde, die im Institut fiir Strémungsmechanik der Technischen 
Hochschule Braunschweig in systematischer Weise rechnerisch gewonnen wurden. 


1 L, J. Herrig, J. C. Emery and J. R. Erwin, Systematic two-dimensional cascade tests of NACA 65-Series 
compressor blades at low speeds. NACA TN 3916 (1957). 

2 A. R. Howell, The present basis of axial flow compressor design, ARC R & M 2095 (1942). 

3 B, Eckert, Axialkompressoren und Radialkompressoren. Berlin/Géttingen/Heidelberg 1953. 

4 Siehe FuBnote 1 von Seite 202. 
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Bei groBeren Zustr6mwinkeln, z. B. bei 6, = 160° in Abb. 2a, andert sich der Gesamtverlust- 
beiwert stark mit dem Belastungsverhaltnis, und das optimale Belastungsverhaltnis (Minimal- 
verlust) ist dabei deutlich erkennbar, wahrend bei kleineren Zustrémwinkeln, 6, < 140°, sich der Ge- 


# samtverlustbeiwert praktisch kaum mit dem Belastungsverhaltnis andert. Deshalb hat die richtige 


Wahl des Belastungsverhdltnisses nur bei gréBeren Zustrémwinkeln 6, > 140° Bedeutung. Die 
weiteren Ergebnisse sind in Abb. 3 bis 8 in einer Form wiedergegeben, wie man sie unmittelbar 
fir die Auswahl der optimalen Tandem-Gitteranordnungen verwenden kann. Abb. 3 dient der 


| Auswahl fiir das optimale Belastungsverhiltnis und fiir die Schaufelprofile der beiden Gitter beim 


ye 


70° 


60° 


40° 


30° 


20° 


110° 720° 730° 740° 150° 760° 110° 


Abb. 3. Auswahldiagramm fiir die Schaufelprofile des Gitters I und II und fiir 
das optimale Belastungsverhiltnis. 

Optimales Belastungsverhaltnis (AW, 1/4 Wy)*; 

Abgrenzung der Schaufelprofile. 


Bereich | A 


Profil des Gitters I NACA 8410 | NACA 4410 | NACA 4410 | NACA 0010 
Profil des Gitters II NACA 8410 | NACA 8410 | NACA 4410 | NACA 4410 


© — Anwendungsbeispiel nach Abb. 9. 


Minimalverlust fiir ein beliebig gegebenes Geschwindigkeitsdreieck, wobei sich das optimale Be- 
lastungsverhaltnis von den Linien konstanter Belastungsverhiltnisse ablesen l48t. Linien kon- 
stanter optimaler Belastungsverhaltnisse sind nur in dem Bereich mit gréBerem Zustrémwinkel 
eingetragen. In dem iibrigen Bereich spielt das Belastungsverhiltnis fiir den Gesamtverlustbeiwert 
fast keine Rolle. Fiir eine giinstige Aufteilung der gesamten Umlenkung scheint in diesem Bereich 
das Belastungsverhaltnis von (AW,;/AW,,)* = 0,5 am besten zu sein. 

In Abb. 4 bzw. Abb. 5 sind die zu den optimalen Anordnungen gehérenden Teilungsverhaltnisse 
des ersten bzw. des zweiten Gitters fiir ein beliebig gegebenes Geschwindigkeitsdreieck angegeben. 
Abb. 6 baw. Abb. 7 dient ferner der Bestimmung der Schaufelwinkel des ersten bzw. des zweiten 
Gitters fiir die optimalen Tandem-Gitteranordnungen. Der Schaufelwinkel wird hierbei mittelbar 
aus den Winkeliibertreibungen der beiden Gitter Afs, und Af,,, die von den beiden Diagrammen 
direkt abzulesen sind, bestimmt mit Hilfe der Beziehungen 


Bs, = 16 + Bs) + ABs; fiir Gitter I 


und 


Bsn =f; + Apsn fiir Gitter IT. 
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Abb. 4. Auswahldiagramm fiir das Teilungsverhaltnis 
des Gitters I, (t/l)*, bei der optimalen Tandem-Gitteranordnung. 
Teilungsverhaltnis des Gitters I; 
Abgrenzung des Schaufelprofiles. 
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Abb. 6. Auswahldiagramm fiir den Schaufelwinkel 


des Gitters I, Bey bei der optimalen Tandem-Gitteranordnung. 
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Abb. 5. Auswahldiagramm fiir das Teilungsverhaltnis 
des Gitters IT, (/D E> bei der optimalen Tandem-Gitteranordnung. 


Teilungsverhaltnis des Gitters II; 
Abgrenzung des Schaufelprofiles. 
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Abb. 7. Auswahldiagramm fiir den Schaufelwinkel 
des Gitters II, Bary bei der optimalen Tandem-Gitteranordnung. 
en 
eat ae Bs ate 4B oo 
Winkeliibertreibung AB 
Abgrenzung des Schaufelprofiles. 
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| Abb. 8 zeigt Linien konstanter Gesamtverlustbeiwerte des Tandem-Pumpengitters beim Minimal- 
| verlust, wobei als Reynoldszahl Re = W,l/y = 5- 10° angenommen wurde. In Bereichen, in 
denen die Schaufelbelastung nicht allzu gro8 ist, kann man statt des Tandem-Pumpengitters auch 
ein gewohnliches einreihiges Gitter benutzen. Die strichpunktierten Linien in Abb. 8 zeigen Zu- 
} satzlich auch die Verlustbeiwerte des einreihigen Gitters beim Minimalverlust. Aus dem Vergleich 
der beiden Kurvenscharen laBt sich weiter eine gestrichelte Abgrenzungslinie eintragen. Oberhalb 
dieser Linie hat ein Tandem-Pumpengitter einen geringeren Verlustbeiwert als ein entsprechendes 
} einreihiges Gitter, wahrend unterhalb dieser Linie ein einreihiges Gitter giinstiger ist. Aus dem 
| Verlauf der Abgrenzungslinie erkennt man, 
daB der Bereich, in dem das Tandem-Pum- 
pengitter hinsichtlich der Verluste giinstiger 
ist als das einreihige Gitter, wider Erwarten 
groB ist und auch Falle geringer Umlenkung 


oNiE 


V 
Profil: NACA 8410 (E)=a9 ( 
I 


20° 


110° 720° 130° 140° 150° 160° 170° b 

Abb. 9. Anwendungsbeispiel fiir eine optimale Tandem-Gitter- 
anordnung ohne Beriicksichtigung der gegenseitigen Inter- 
ferenz. — a) Gegebenes Geschwindigkeitsdreieck: 6B, = 160° 


Verlustbeiwerte als Tanden-Pumpengitter (C777 +¢y,11)* und f, = 100°, b) Gesuchte Gitteranordnung mit Hilfe von 
Abb. 3 bis 8, zugehériger Gesamtverlustbeiwert 


Abb. 8. Verlustbeiwerte des Tandem-Pumpengitters beim Minimalverlust 
und Vergleich mit einreihigem Gitter. 


—— .—_———. Verlustbeiwerte als einreihiges Gitter Cy) z. @. 


. dia Sea Abgrenzung der giinstigeren Gitteranordnung, oberhalb Gritty w* = 070. 
der Kurve Tandem-Pumpengitter, unterhalb der Kurve > > 
einreihiges Pumpengitter giinstiger. 


mit einschlieBt, so daB das Tandem-Pumpengitter in Wirklichkeit einen sehr grofen Anwendungs- 
bereich besitzt. In Abb. 8 sind noch zwei weitere Begrenzungslinien eingetragen, die den Formeln 


car = se =19-baw. 15 (3.1) 
entsprechen, wobei cy, den auf die Grundgeschwindigkeit W. bezogenen Auftriebsbeiwert der 
Gitterschaufel bedeutet. Herkémmlicherweise wird die obere Begrenzungslinie cy Vt a 1,5 als 
die Belastungsgrenze eines einreihigen Gitters angesehen?*. Wie man sieht, umschhieBt diese 
- Linie auch noch einen breiten Bereich, bei dem das Tandem-Pumpengitter schon giinstigere Werte 
liefert. Nach den hier gewonnenen Ergebnissen wiirde sich die Kurve cy, 1/t = 1,2 als Belastungs- 
grenze besser eignen. Oberhalb dieser Grenze ist es jedenfalls empfehlenswert, das Tandem- 
Pumpengitter zu verwenden, um damit die Méglichkeit der Verbesserung des Wirkungsgrades aus- 


zunutzen. Rites 
Abb. 9 stellt ein Anwendungsbeispiel fiir eine groBe Umlenkung von Af = 60° mit Hilfe der 


hier angegebenen Auswahldiagramme dar. 


1 B. Eckert, Axialkompressoren und Radialkompressoren. Berlin/Géttingen/Heidelberg 1953. 
2 Siehe FuBnoten 1 und 2 von Seite 204. 
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Solche Untersuchungen iiber Tandem-Pumpengitter starker Umlenkung kann man auch mit 
den Gitterergebnissen anderer Profilserien durchfiihren. Es ist zu vermuten, da die vorliegenden 
Ergebnisse weitgehend unabhangig von der untersuchten Profilserie sind und ziemlich allgemeine 
Kigenschaften des Tandem-Pumpengitters wiedergeben. 


4, Experimentelle Untersuchungen iiber den Interferenzeinflu8 des Tandem-Pumpengitters bei 
verschiedenen relativen Stellungen der beiden Gitter. Um die Anderung der Gittereigenschaften 
des Tandem-Pumpengitters, namlich die Anderung der Umlenkung, der Verluste und der Druck- 
verteilung, bei festen Anordnungen der beiden Gitter experimentell klarzustellen, wurden Wind- 
kanalmessungen im kleinen Gitterwindkanal des Institutes durchgefiihrt. Abb. 10 stellt die MeB- 
strecke des Kanales dar. Die beiden Gitter bestehen aus je 11 Schaufeln von Profilen NACA 8410 
(8% Wélbung bei 40% Wolbungsriicklage und 10% Maximaldicke). Die Schaufeln haben eine 
Profiltiefe von 1 = 60 mm und eine Schaufelhéhe von 


300 mm. Sie wurden aus Stahl hergestellt und besitzen x ys 
bei der verwendeten Anblasegeschwindigkeit hydrau- “Ss \ SG 
lisch glatte Oberflache. Ferner sind zwei der Schaufeln Pro=nge WACHeHD 


in ihrem Mittelschnitt mit insgesamt 26 Druckmefstellen 
versehen, von denen je 13 auf einer Schaufelseite an- 
geordnet sind. Gitter I 


Abb. 10. Einbau des Tandem-Pumpengitters im Versuchsstand. 1 Verstellbarer Abb, 11. Gemessene Gitteranordnung des Tandem-Pum- 
Boden des Zustrémkanales; 2 Absaugung zum Gebliise; 3 Drehscheibe; 4 Wand- pengitters. Bei unendlich groBem Abstand a/l = o er- 
bohrung zur Messung des statischen Drucks p,; 5 Erster Schaufeltriger fir méglicht die Anordnung einen Abstrémwinkel f,, = 90° 
Gitter I aus 11 NACA 8410 Schaufeln; 6 Zweiter Schaufeltrager fiir Gitter II beim Zustrémwinkel 6,, = 140° mit einem Belastungs- 
aus 11 NACA 8410 Schaufeln, parallel und senkrecht zur Gitterfront verstelibar; verhaltnis (4W,, 7/4 W,) = 0,5. (Fubzeiger 0 entspricht 


7 MeBebene fiir die Nachlaufmessungen, MeBabstand e = 30 mm (0,5 1). dem Fall ohne Interferenz.) Bei den Messungen sind nur 


der Abstand a und die Versetzung h variiert worden. 


Die Gitteranordnungen beider Gitter sind hierbei so bestimmt, daB die beiden Gitter bei un- 
endlich groBem Abstand, namlich ohne Interferenz, ein gegebenes Geschwindigkeitsdreieck von 
Bio = 140° und (3) = 90° verwirklichen (vgl. Abb. 11). Fiir die Aufteilung der gesamten Um- 
lenkung wird als optimales Belastungsverhaltnis (AW,;/AW,,)* = 0,5 auf Grund der Ergebnisse 
in Abb. 3 ausgewahlt. Infolgedessen soll das erste Gitter die Umlenkung Af,) = 19,2° bei Zu- 
strémwinkel /,, = 140°, das zweite Gitter die Umlenkung Afio = 30,8° bei Zustrémwinkel 
Peo = 120,8° als Einzelgitter erméglichen. In Hinblick auf die praktischen Anwendungen wurde 
fiir beide Gitter die gleiche Teilung, und zwar t/l = 0,75 gewahlt. Von den Voruntersuchungen 
sind die zum Verwirklichen der gegebenen Umlenkung als Einzelgitter erforderlichen Schaufel- 
winkel und die zugehérigen Verlustbeiwerte ohne Interferenz bestimmt worden. Die gemessene 
Tandem-Gitteranordnung ist in Abb. 11 dargestellt. 

Zur Festlegung der gegenseitigen Stellung der beiden Gitter dienen hierbei zwei Parameter, 
und zwar der Abstand a und die Versetzung h. Der Abstand ist die Spaltweite zwischen den 
beiden Gittern und wurde bei den Messungen auf a/l = 0,15; 0,30; 0,50 und 0,80 eingestellt. 
Die Versetzung h wird hierbei nach der geometrischen Lage zwischen der Delle der ersten 
Schaufel und der Vorderkante der zweiten Schaufel festgelegt. Wenn die Delle von der Hinterkante 
der Schaufel I, deren Bahnlinie naherungsweise als eine gerade Linie mit Abstrémwinkel 120,8° 
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(== Bzo) angenommen sei, direkt die Vorderkante der Schaufel II trifft, betragt die Versetzung un- 
abhangig vom Abstand immer h = 0. Die weitere Versetzung wird nach der Verschiebung der 
Schaufel II in Richtung der Saugseite gegeniiber der Versetzung h = 0 angegeben und hat dabei 
positives Vorzeichen. Die Messungen sind bei Versetzungen von h/t = 0; 0,2; 0,4; 0,6 und 0,8 
durchgefiihrt worden. Im iibrigen stimmt die Versetzung h/t = 1,0 wieder mit der Versetzung 
h/t = 0 tberein. 

Die Schaufeln der beiden Gitter befinden sich jeweils in einem Schaufeltrager. Der erste Schaufel- 


' trager ist mit der Drehscheibe des Kanales fest verbunden, wahrend der zweite Schaufeltrager, im 


Abstand und in der Versetzung verstellbar, auf dem ersten Schaufeltrager aufgeschraubt wird. 
Vor dem Eintritt ins Tandem-Pumpengitter wird am oberen und unteren Kanalende die Wand- 
grenzschicht der Zustrémung durch Schlitze abgesaugt. Es wurden Nachlaufmessungen in einer 
Entfernung von 50% der Schaufeltiefe im Mittelschnitt hinter dem Gitter II und Druck- 
verteilungsmessungen lings der Schaufelkontur der beiden Gitter durchgefiihrt. Die Nachlauf- 
messungen dienen der Bestimmung des Abstrémwinkels 6, und der Verlustbeiwerte der beiden 
Gitter Cy,1 und Cy, 1, die sich hierbei nach den Formeln ermitteln lassen: 
t 


Se ee FD) 
Sra sia ee dy , (4.1) 
oe g Me 
it) 
A t 
1 _— 
i ; one : seer 83(¥) dy. (4.2) 
3 Mi 2g Va 


Bei Versetzung h/t = 0,4 und 0,6 kann man bei der Nachlaufmessung hinter dem Gitter IT 
deutlich zwei getrennte Dellen von den beiden Gittern beobachten. In diesem Fall kann man ohne 
weiteres die Verlustbeiwerte der beiden Gitter getrennt ermitteln. Bei Versetzung h/t = 0,2 und 
0,8 iiberlagern sich die beiden Dellen an ihren Randern. Aber es ist dabei auch méglich, diese 
Dellen mit ausreichender Genauigkeit voneinander zu trennen und die Verlustbeiwerte der beiden 
Gitter einzeln zu ermitteln, da die Form der Dellen an ihrem Rand aus den anderen Messungen 
geniigend genau geschatzt werden kann, und der zusammenfallende Verlust von den beiden Gittern, 


hag {402 
° Verseteung 
88 : t 
98 mee G2 o4 06 08 0 G2 


Abb. 12. Anderung des Abstrémwinkels 6, und des Gesamtzirkulationsverhialtnisses I'/I", 
mit der Versetzung bei verschiedenen Abstanden. 


Ag, + Ag,» in die einzelnen Anteile aufgespalten werden kann. Dagegen kommen die beiden 
Dellen bei Versetzung h/t = 0 vollkommen zusammen und man kann daraus nur den Gesamt- 
verlust erhalten. Nach den Mefergebnissen andert sich der Verlust des Gitters I beinahe sinus- 
férmig mit der Versetzung. Deshalb ist es nicht schwierig, von den bereits ermittelten Ergebnissen 
bei Versetzung h/t = 0,2; 0,4; 0,6 und 0,8 den Verlust des Gitters I bei Versetzung h/t = 0 zu 
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interpolieren. Der Verlust des Gitters II bei Versetzung h/t = 0 ist hierbei als der Restteil des 


Gesamtverlustes zu verstehen. 
Bei allen Messungen des Tandem-Pumpengitters betrug der Staudruck vor dem Gitter I 


200 mm WS, entsprechend einer Austrittsgeschwindigkeit des Gitters I von W, = 43 m/s und 
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Abb. 13. Verlustbeiwerte des Tandem-Pumpengitters bei verschiedenen Abstiinden in Abhangigkeit von der Versetzung. 
Rey = W,. UD S i Tf . 108; 
Reqy = W,. Up = SS + 105. 


des Gitters II von W; = 37 m/s. Die auf die Austrittsgeschwindigkeit bezogene Reynolds-Zahl 
betragt daher 


W,1 


I = == Ile oll fur Gitter I, 


W, 
Rey = —- = 1,5 - 108 fir Gitter II. 


v 
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Auf Grund der Nachlaufmessungen hinter dem Gitter II zeigt Abb. 12 die Anderung des Ab- 
strémwinkels Bz mit der Versetzung bei verschiedenen Abstinden. Die rechte Nebenskala dient 
der Bestimmung der entsprechenden Anderung der Gesamtzirkulation nach der Formel 


LP AWu _ 8 By — ete B, 
Py ~~ AWuo ctg By — ctg Byy ° 


(4.3) 


Die Abweichung des Abstrémwinkels von 8, = 90° ist reiner InterferenzeinfluB zwischen den beiden 
Gittern. Diese Ergebnisse zeigen, daB der Abstrémwinkel im wesentlichen unabhangig vom Ab- 
stand, dagegen sehr stark abhaingig von der Versetzung ist. Bei Versetzung h/t = 0, bei der 
also die Delle der Schaufel I direkt die Vorderkante der Schaufel II trifft, vermindert sich die 
Umlenkung oder die Zirkulation erheblich. Diese Erscheinung kann darauf zuriickgefiihrt werden, 
daB die Grenzschicht auf der Schaufelkontur des Gitters II wegen des Delleneinflusses eher zur 
§ Ablésung neigt und dadurch zu einer Abnahme 

- der Zirkulation fiihrt. Bei Versetzung h/t = 0,2 
q und 0,4 bei denen also eine Delle der Schaufel I 
i sich direkt an der Druckseite der Schaufel IT 
befindet, nimmt die Gesamtzirkulation im allge- 
meinen um etwa 1% zu. Bei Versetzung h/t = 0,8 
verschwindet der InterferenzeinfluB vollkommen, 
und die Abstrémungsrichtung bzw. die Gesamt- 
zirkulation stimmt mit dem Fall ohne Interferenz 
genau iiberein. Die Schwankung des Abstrém- 
winkels klingt mit zunehmendem Abstand all- 
mahlich ab, aber selbst bei dem verhaltnismaBig 
groBen Abstand a/l = 0,8 ist die Schwankung 
noch ziemlich grof. 


In Abb. 13 ist die Anderung der Verluste des 


~Tandem-Pumpengitters dargestellt. Verlustbei- ' 
werte des Gitters I, des Gitters II und der beiden ' 
Versetaung®: 0 02 Of O06 
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fi VI min AN | 
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Gittern zusammen sind bei verschiedenen Ab- \ 

a =F 5 Ue LO ARE Wie ty 
standen iiber der Versetzung aufgetragen, wobei Schautel IT 

eae Vergleich auch die Verlustbeiwerte ohne Abb. 14. Gegenseitige Lage beim InterferenzeinfluB des Gitters II 


Interferenz angegeben sind, die von der Ver- auf Gitter I. 
: Ge . . ° Ort der Vorderkante der Schaufel II, bei dem das 
setzung nicht abhangen. Die Verlustbeiwerte des Cliter E bel gegolenons Abstand elven Merial! 
| ersten und zwei itters si in I = verlust erreicht; 
7 a % ren. Gitt sind in ihrem Ver —-— @—-— Ort der Vorderkante der Schaufel II, bei dem 
lauf deutlich voneinander verschieden, woraus sich das Gitter I bei gegebenem Abstand einen Mini- 


: : . ; Ivedastrercoicat! 
wesentliche Erkenntnisse itiber die Entstehung parce ieee 


des Interferenzeinflusses ergeben. Die Verlust- 

beiwerte des Gitters I andern sich sinusférmig um einen Mittelwert, der mit dem Verlustbeiwert 
ohne Interferenz gut iibereinstimmt, und die Phase der Verlustkurve verschiebt sich mit zu- 
nehmendem Abstand allmahlich in Richtung der Saugseite. 


Abb. 14 stellt die Orte der Vorderkante der Schaufel II dar, bei denen Gitter I bei gegebenem 
Abstand durch den Interferenzeinflu8B vom Gitter IT den Minimal- bzw. Maximalverlust erreicht. 
Die Orte liegen praktisch auf einer Geraden, infolgedessen andert sich die Phase linear mit dem 
Abstand. Bei sehr kleinem Abstand ist der Verlust des Gitters I bei Versetzung h/t = 0 maximal 
und bei Versetzung h/t = 0,5 minimal. 

Aus der regelmafBigen sinusformigen Anderung der Verluste des Gitters I und auch aus den 
gemessenen Druckverteilungen an Gitter I in Abb. 15 beim Abstand a// = 0,15 kann man den 
InterferenzeinfluB des Gitters II auf das Gitter I erkennen: Gitter II beeinfluBt zunachst die 
Druckverteilung des Gitters I je nach der gegenseitigen Stellung beider Gitter, und diese Anderung 
der Druckverteilung verursacht dann die Schwankung der Verluste des Gitters I. Zum Beispiel 
hat die Druckverteilung des Gitters I bei Versetzung h/t = 0,2 wegen des verschlechternden Inter- 
ferenzeinflusses des Gitters II in der Nahe der Hinterkante starken Druckanstieg, was einen 
Zuwachs der Verluste zur Folge hat. Demgegeniiber verlauft die Druckverteilung bei Versetzung 
h/t = 0,6 und 0,8 ahnlich derjenigen ohne Interferenz, und infolgedessen besitzt Gitter I dafiir 
kleinere Verluste. Wie schon oben erwdhnt, kann also die Strémung in Gitter I potential- 
theoretisch erklart werden. 
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Bei Gitter II andern sich die Verluste mit der Versetzung derart, daB die Verluste bei jedem 
Abstand zwischen Versetzung h/t = 0 und 0,2 erheblich (um 50% bei a/l = 0,15) zunehmen, 
wihrend sie bei den iibrigen Versetzungen etwa konstant bleiben. Aus dem Ergebnis, daf die 
Verluste immer zwischen Versetzung h/t = 0 und 0,2 zunehmen, d. h. die Phase der Schwankung 
sich nicht mit dem Abstand andert, ist die Folgerung zu ziehen, daf die Verlustdellen von den ersten 
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Abb. 15. Druckverteilungen des Gitters I beim Abstand a/l = 0,15 und verschiedenen Versetzungen h/t. 


Schaufeln fiir die Verluste des Gitters II eine entscheidende Rolle spielen. Der Verlauf der Ver- 
luste des Gitters II laBt sich nicht durch eine potentialtheoretische Betrachtung erfassen, bei 
der also die Anderung der Verluste nur durch die Anderung der Druckverteilung hervorgerufen 
wird, sondern nur durch eine grenzschichttheoretische Betrachtung, welche aufzeigt, wie eine 
Delle die Grenzschichtbildung der nachstehenden Schaufel und damit die Lage der Ablésungspunkte 
beeinfluBt. 

Die Strémung in einer Delle enthalt eine kleinere kinetische Energie als die Strémung aufer- 
halb der Delle. Deshalb wird beim Uberwinden des von der AuSenstrémung aufgepragten Druck- 
anstieges die Strémung in der Delle eher als die gesunde Strémung ihre kinetische Energie er- 


‘ 
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schépfen und infolgedessen eher ablésen, so daB sich daraus ein gewisser Zuwachs der Verluste 
ergibt. Dieser ausgepragte InterferenzeinfluB einer Delle in der Zustrémung wird in Ziff.6 noch 
ausfihrlicher behandelt. 

Demgegeniiber erreicht Gitter IT noch kleinere Verluste als ohne Interferenz bei Versetzung 
h/t = 0,4; 0,6 und 0,8, bei denen die Schaufel IT sich in gesunder ungestérter Strémung befindet. 
In Abb. 16 sind die Druckverteilungen des Gitters II bei Abstand a/l = 0,15 bei verschiedenen 
Versetzungen dargestellt. Auer bei Versetzung h/t =0 weichen die Druckverteilungen kaum 
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Abb. 16. Druckverteilungen des Gitters II beim Abstand a/l = 0,15 und verschiedenen Versetzungen h/t. 


voneinander ab. Aus der Druckverteilung bei Versetzung h/t = 0 ist zu beobachten, da der Druck 
auf der Saugseite in der Nahe der Hinterkante infolge der Ablésung der Grenzschicht nicht weiter 
ansteigt. Aus dieser Ablésung der Grenzschicht folgt unmittelbar der erhebliche Zuwachs der 
Verluste. Alle anderen Druckverteilungen bei Versetzung h/t = 0,2; 0,4; 0,6 und 0,8 sind der 
Druckverteilung ohne Interferenz ahnlich. Die Ergebnisse der Druckverteilungsmessungen an 
-Gitter II liefern also qualitativ dieselbe SchluBfolgerung wie die Ergebnisse der Nachlaufmes- 
sungen. 

in Abb. 17 ist der Gesamtverlustbeiwert Cy,; + Cy, tber dem Zirkulationsverhaltnis ["/I", 
aufgetragen. Bei festem Abstand wird dabei durch Anderung der Versetzung eine geschlossene 
Kurve durchlaufen. Die Werte der Versetzung sind an die einzelnen Punkte angeschrieben. Mit 
zunehmendem Abstand schrumpft die geschlossene Kurve allmahlich zusammen und ist schlieBlich 
bei unendlich groBem Abstand nur noch ein Punkt von J'/J’) = 1,0 und (Cy,1 +¢y,11)p = 9,136. 
Das Diagramm stellt den Vorgang des gegenseitigen Interferenzeinflusses sehr tibersichtlich dar. 
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Daraus ist fiir den Entwurf des Tandem-Pumpengitters die Folgerung zu ziehen, daB die Ver- 
setzung h/t = 0,8, bei der die beiden Schaufeln sich wie ein Spaltfliigel mit einer Spaltweite 
von /, der Teilung anordnen, und die Schaufel IT sich in der nicht abgelosten gesunden Strémung 
an der Druckseite der Schaufel I befindet, bei einem in der Praxis gewéhnlich vorkommenden 
Abstand von a/l = 0,15 bis 0,5 in folgender Hinsicht am giinstigsten ist: erstens erreicht der Ge- 
samtverlust dabei den Minimalwert, der etwa 10° kleiner ist als der Verlust ohne Interferenz, und 
zweitens stimmt die Gesamtzirkulation mit derjenigen ohne Interferenz gut iiberein, so daB keine 
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Abb. 17. Gesamtverlustbeiwert. des Tandem-Pumpengitters Cy 7 + Cy,11 ber dem Zirkulationsverhaltnis T/T, fiir verschiedene Abstande 
und Versetzungen. 


Die Zahlen an den Kurvenpunkten geben die entsprechende Versetzung an. 


Korrektur fir den Abstrémwinkel nétig ist. Der Abstand ist hierbei fast gleichgiiltig und kann 
beliebig ausgewahlt werden, solange er zwischen. a/] = 0,15 und 0,5 liegt. 


Nach den in Ziff. 3 gewonnenen Ergebnissen haben die optimalen Anordnungen der beiden 
Gitter ohne Beriicksichtigung der gegenseitigen Interferenz nicht unbedingt die gleichen Teilungs- 
verhaltnisse. Wenn die Teilungen der beiden Gitter verschieden sind, andert sich die relative Lage 
der beiden Schaufeln zueinander nicht periodisch mit der Teilung, und es kann sich der Interferenz- 
einflu8 an einer Stelle, bei der die beiden Schaufeln sich wie ein Spaltfliigel anordnen, giinstig und 
an anderer Stelle, bei der sich die Schaufel II im Nachlauf der Schaufel I befindet, ungiinstig aus- 
wirken. Um solches mit der Teilung nicht periodisches Strémungsfeld zu vermeiden und nur die 


giinstige Wirkung des Interferenzeinflusses auszunutzen, sollten die Teilungen der beiden Gitter 
méglichst gleich sein. 


In Abb. 18 ist ein optimales Tandem-Pumpengitter unter Beriicksichtigung der gegenseitigen 
Interferenz fiir ein Geschwindigkeitsdreieck von 6, = 160° und Af = 60° dargestellt, wobei die 


Profilformen, Teilungsverhaltnisse und Schaufelwinkel der beiden Gitter bereits nach dem Beispiel 
in Abb. 9 ohne Beriicksichtigung der Interferenz bestimmt wurden. 
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5. Ein Berechnungsverfahren zum Entwurf eines optimalen Tandem-Pumpengitters fiir ein 
gegebenes Geschwindigkeitsdreieck. Neben den experimentellen Untersuchungen wurde das Ent- 
wurfsproblem eines optimalen Tandem-Pumpengitters fiir ein gegebenes Geschwindigkeitsdreieck 
auch auf theoretischem Wege unter Beriicksichtigung der gegenseitigen Interferenz untersucht. 
Diese Untersuchungen haben ebenfalls das Ziel, die giinstigste gegenseitige Stellung der beiden 
Gitter zu bekommen. Zunichst wurde hierbei ein potentialtheoretisches Rechenverfahren zum 
Entwurf eines Tandem-Pumpengitters unter Beriicksichtigung der gegenseitigen Interferenz (erste 
| Hauptaufgabe des Tandem-Schaufelgitters) entwickelt, und zwar als Erweiterung des Verfahrens, 
welches NV. Scholz! fiir den Entwurf eines einreihigen Schaufelgitters mit Hilfe der Singularitaten- 
/ methode angegeben hat. Weiter wurde die Anderung der Verluste des berechneten Tandem- 
1 Pumpengitters mit der Versetzung auf Grund eines grenzschichttheoretischen Rechenverfahrens2, 
ermittelt. Daraus laBt sich die optimale Versetzung (Minimalverlust) theoretisch bestimmen. 


a) Potentialtheoretisches Berechnungsverfahren zum Entwurf eines Tandem- 
- Pumpengitters. Das Berechnungsverfahren legt die Singularitatenmethode zugrunde, bei der 
eine Schaufel durch eine passend gewahlte Wirbel- und Quellbelegung auf der Skelettlinie ersetzt 
wird, und sich das Strémungsfeld aus der Uberlagerung der Translationsstrémung und der 
resultierenden, von den Singularitétenbelegungen induzierten Geschwindigkeiten ergibt. Die 
gegenseitige Interferenz des Tandem-Pumpengitters wird dadurch erfaB8t, daB® sich die indu- 
zierte Geschwindigkeit von der Singularitatenbelegung der anderen Schaufelreihe mit der gegen- 
seitigen Stellung der beiden Gitter zueinander andert. Infolgedessen besteht die Aufgabe im 
wesentlichen in der Berechnung der induzierten Geschwindigkeiten von der anderen Schaufelreihe. 
| Bei der Berechnung sind Geschwindigkeitsdreieck, Schaufeltiefe und Gitterteilung der beiden Gitter 
sowie Abstand und Versetzung zwischen den beiden Gittern vorgegeben. Im iibrigen wurde hierbei 
nur der fiir die Praxis hauptsachlich interessierende Fall behandelt, bei dem die beiden Gitter die 
gleiche Teilung besitzen. 


Abb. 19 zeigt das gegebene Geschwindigkeitsdreieck fiir die Berechnung. Jedes der beiden 
Gitter besitzt ein eigenes Koordinatensystem (vgl. Abb. 20). Die Nullpunkte der beiden Koor- 


i -dinatensysteme befinden sich im Mittelpunkt der Schaufelsehne und haben voneinander den Ab- 


stand b, bzw. h, senkrecht bzw. parallel zur Schaufelfront. Die «-Richtung des Koordinatensystems 
Ibzw. IList parallelzur Grundstrémung W,,, 
des Gitters I bzw. W.,,, des Gitters IT, und 
es gilt 


dj =Boy—Z G=lund Il). (5.1) 


GitterI 


» 
ee igne 
Sy 


Gittertront 


Abb. 19. Geschwindigkeitsdreieck fiir die Berechnung Abb. 20. Koordinatensysteme der beiden Gitter. 
des Tandem-Pumpengitters. 


4 


Damit sind die zwei Koordinatensysteme I und II fiir die beiden Gitter festgelegt, und die Orts- 
vektoren der beiden Systeme sind nach der Transformationsformel 


Zl = 2p + Zr e—%d mit 6 =A; —Ay (5.2) 
miteinander verbunden. 


1 N. Scholz, Stromungsuntersuchungen an Schaufelgittern, Teil I: Ein Berechnungsverfahren zum Ent- 
wurf von Schaufelgitterprofilen, VDI-Forschungsheft 442 (1954). 
2 Siehe FuBnote 1 von Seite 202. 
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Zur Vereinfachung werden die Singularitaten statt auf der Skelettlinie auf der Schaufel- 
sehne angeordnet. Hierdurch wird die Anwendbarkeit dieses Verfahrens auf Profile mit maBiger 
Wolbung beschrankt. Die passend gewahlte Wirbelbelegung y(x) auf der Schaufelsehne erzwingt 
die Umstrémung der Skelettlinie, wahrend eine zusitzliche Quellbelegung q(x) die Schaufel auf die 
vorgegebene Dickenverteilung aufblast. Aus der Wirbelbelegung folgt die Gleichung fiir die Ge- 


samtzirkulation 


1j/2 
r= J, v(x) dx, (j=1, 1), (5.3) 
und die Quellbelegung muB die Bedingung fiir eine geschlossene Schaufelkontur 
1j/2 
0,= J qj(x;)dx,;=90 (J =I, 0) (5.4) 
ae 


erfiillen. Die Komponenten der induzierten Geschwindigkeiten u und v lassen sich als Summe ver- 
schiedener induzierter Geschwindigkeiten folgendermaBen darstellen: 


uj = (Ujny + Ujng) + (Ujry + Ujrg) 1 (Ujry + Uj 74) > (5.5a) 

05 = (Yjzy + Eq) + jry + Rg) + Ory + % 54)» (5.5b) 

mit j = I und II, wobei die Indizes in ihrer Reihenfolge anzeigen, erstens auf welcher Schaufel, 

zweitens von welcher Schaufel und drittens von welcher Singularitatenbelegung die Induktion erfolgt. 

1. Induzierte Geschwindigkeit von der Einzelschaufel w;_. Die induzierte Geschwindigkeit im 

Punkt z von der Singularitatenbelegung der zu dem Punkt gehérenden Schaufelreihe we erhalt 
man bekanntlich zu 


ta; 452? id; 
Tye) = “Sr | laos) + én ete [2 | —2)] ae (5.6) 
13/2 


Wenn der Aufpunkt sich direkt auf der Sehne der betrachteten Schaufelreihe befindet, wird der 
Integrand am Aufpunkt unendlich. Das Integral ist dann als Cauchyscher Hauptwert zu verstehen. 
Laf&t man in (5.6) die GréBe z > x, also y > 0 gehen, so liefert der Grenziibergang fiir die komplexe 
Geschwindigkeit auf der «-Achse 
1j/2 : 
7 sles) — i aye) Nene 
Wiley) = APSE LS | [alej) +t yi(xf)] Ctg | — (&} — 3) 
—ij/2 
Dabei gilt das positive Vorzeichen fiir die obere Schaufelseite und das negative Vorzeichen fiir die 
untere Schaufelseite. 
Im Grenziibergang t—> 0 liefert (5.7) die induzierte Geschwindigkeit der Einzelschaufel. Man 

erhalt unter Benutzung der Reihenentwicklung des hyperbolischen Kotangens 


dx;. (5.7) 


1j/2 
oa yj(xj) — 1 qy(x;) 1 ‘ / . , d. f 
CC ee ee Ue | late) + inj) Se (5.8) 
yp : 
also in Komponenten : 
14/2 
as __ 1 vi(xj) 1 qy(x} , 
js = Uppy + Hyeg = tO 4 et (5.8a) 
—lj/2 
1;/2 
a4 i . yj(x7) ‘3 qj(xj 
te Yep tise — oy | ee del ch eae (5.8b) 
—lj/2 


mit. j == 1 und If. 


2. Induzierte Geschwindigheit vom Restgitter w;z. Wenn man von der induzierten Geschwindig- 
keit eines ganzen Gitters w;¢ nach (5.7) die induzierte Geschwindigkeit der Einzelschaufel Wz 
nach (5.8) abzieht, erhalt man daraus die induzierte Geschwindigkeit des Restgitters w;p auf der 
«-Achse der betrachteten Schaufelreihe, wodurch das ganze Integral zweckmabig in zwei Anteile, 
mit und ohne Singularitat, aufgespalten werden kann: 

one ‘i? id; 
Ha Wo — We = ay | (ales) + ir) {ta —%5)| 


—lj;/2 


t 
Wy ev (xy — x2) 
mt et Aj(xj x5) 


dx; . (5.9) 


e 
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1 Das in der geschweiften Klammer von (5.9) stehende Glied ist nun eine regulére Funktion mit den 


Variablen A; und (x; — x))/t. Fiihrt man die folgende EinfluBfunktion F ein: 


$ ag i Xj — xy t 
F,=R, +il, =e Cig [xe 4 i) ey (5.10) 
wobei 


muh, < OFS: a 
cos Aj Sin (2 I ! cos iy) sin A; sin (2 n i sin i) 


R; = Be eer as: c Tpar crm (5.10a) 
Co} (2 ca eee cos ii) — cos (2 nm I ; ! sin i) ge 
und 
sin Aj Sin (2 nm pigs! cos Aj — cos Aj sin ( I 2 —t sin Aj 
I; =z - 7 (5.10b) 
Xj — Xj — x; z 
Coj (22 ! cos iy) —e0s Gees sin i) 


so lassen sich die Komponenten der induzierten Geschwindigkeit des Restgitters im eigenen 
Koordinatensystem wie folgt ausdriicken: 


1;/2 17/2 
1 / / ik / ¢ 
pe Way Fyag = — | abs) Tyo | qy(x}) R, dx, (5.11a) 
—l;/2 —lj/2 
1j/2 17/2 
1 7 fs 1 4 tf 
Oin = tyny + ng = —5- | v(x) Ry dx — = | a(x) I; dx, (5.11b) 
—ijj2 —Ij/2 


= mit 7 — Asuna LI, 


_ 3. Induzierte Geschwindigkeit von der anderen Schaufelreihe w;;. Aus der Grundgleichung des 
| induzierten Geschwindigkeitsfeldes nach (5.6) ergibt sich die vom Gitter I auf die Schaufelsehne 
_ des Gitters II induzierte Geschwindigkeit zunachst beziiglich des Koordinatensystems I zu 


Iy/2 


= ey : met Ay 
Wy. 92) = 5 | a(x.) +iy,(x,)] Ctg (a, — x) | dx, . (5.12) 
—i/2 
Mit Hilfe der Translationsgleichungen 
org ee Ot = (by = 1g) ete 2, C3"? 
und = 
w(z,) = w(z,,) e? 
| lautet (5.12) im eigenen Koordinatensystem II 
wu (x11) = wy y(21) i 
1/2 
iA ; 
ae | [a (&)) + i yy(x)] Cty = (by + iby + xy e' —x,6t) dey. (5.13) 


—iyJ2 
In gleicher Weise kann man die vom Gitter II auf die Schaufelsehne des Gitters I induzierte Ge- 
schwindigkeit wy ,(x1) beziiglich des eigenen Koordinatensystems I ableiten: 
3 Iyy/2 
Sy k < . ° 
wy y(%) = Se | [au (% 1) + 17 (%)] Ctg = (by + i hy + ay et — xy e'71) day. (5.13) 
—hyq/2 


Das hyperbolische Glied der Integranden in (5.13 a) und (5.13) ist wieder eine regulare Funktion 
auf der Schaufelsehne und es ergibt sich eine neue EinfluBfunktion F'; mit den Variablen bo, ho, 


Ay Aq» %, und a: 
F,=R, +i 1, =Ctg = (by + i hy + a1 e841 — xy e*“1) , (5.14) 
15 
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wobei gilt 
es Gin2 Xx (5.14 a) 
J Gof 2X—cos2Y’ 


~ Coj 2X — cos 2 Ney 
ae == (by + x, cos A,, — x, cos Ay) > 
ee = (hy + %,, sinA,, — x, sind,) . 


Mit Hilfe dieser neu definierten EinfluBfunktion sind die induzierten Geschwindigkeiten in (5.13 a) 
und (5.13b) umzuschreiben in 


I1r/2 
aye, ; ; 
Wr s(%) = uy — ims =; * (an (%) + 8 Yp(%q)] (Ry +t Lz) day» (5.15 a) 
pele. 
ee 
wn s(%q) = Uy — tons = =; [ay(xy) + i yy(%)] (Ry +t Iz) dx. (5.15 b) 
—Iy/2 


Fiir die numerische Auswertung mu8 man die Gleichungen weiter in die Komponenten zerlegen. 
Hierbei gilt fiir Gitter I 


in A A 
urs y(%) = 5 | Yan) Ry(Xp %q) dey + SS | Val%y) Lyla qq) yy, (5-16 a) 


Il II 
in A A 
mrs g(%) = ors | Gly) Ly (% p> %) bq, — ss : | Qr(%) Kyla, %q) AX q » (5.16b) 
ri if 
: in / A 
P15 y(%) = — = : | Vly) Ly %q) Fey + = - | Vnl%) Ry(ep %y) eq» (9-16 ¢) 
ri i 
in A cos A 
ee i dnl) Rylp Xn) dey + 2S | Qnl%n) Iyltp %q) dq (5-164) 
II Il 
und fir Gitter II 
in A A 
unt Jy(ty) = — 2m | yy(%) Ry (ys &q) de, — 954 | ygleq) Ty (2s eq) dx, (5.1 7a) 
I I 
in A A 
Unt g(%y) = — = a i Qn(%q) Ly (%q, yy) day + “S _ | qu(%q) Ry(xq, xy) dy, (5.17b) 
i i 
in A A 
i Sy(%qy) = + = Fi | V1(%) Ly (%p %q) dx — “F = | Val) Ry(xq, xy) day , (5.17) 
i i 
in A A 
vist) =— 0 | gu(%) Ry(%p %q) oxy — | qy(x,) L(x X41) dx, « (5.17) 
1 I 


Damit sind alle Komponenten der induzierten Geschwindigkeiten in (5.5) in einer praktisch 
brauchbaren Form dargestellt. Die numerischen Auswertungen der Integralglieder kénnen nach 
der Simpsonschen Regel durchgefiihrt werden } 2 3, 


Die Schaufelform ergibt sich aus der Uberlagerung einer Wélbungsverteilung y,(x) der Skelett- 
linie und einer Dickenverteilung rx(x), und zwar in der Art, daB die Kontur aus den Enveloppen 
der eingeschriebenen Kreise mit den Halbmessern rg(x) um die Aufpunkte auf der Skelettlinie 
entsteht. 


1 H. Ohashi, Theoretische und experimentelle Untersuchungen an Tandem-Pumpengittern starker Um- 
lenkung. Dissertation der TH Braunschweig (1958). 


* F. Riegels und H. Wittich, Zur Berechnung der Druckverteilung von Profilen, Jahrbuch Deutsch. Luft- 
fahrtforsch. 1942, S. I 120. 


° T. Theodorsen and J. E. Garrick, General potential theory of arbitrary wing section, NACA TR 452 (1933). 
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Die Skelettlinie y, ergibt sich als Stromlinie, indem die betrachtete Schaufel nur mit einer 
| Wirbelbelegung versehen wird, alle iibrigen Schaufeln, auch die von der anderen Schaufelreihe, 
+ aber die vollstandige Singularitatenbelegung erhalten. Hierbei gilt fiir die Neigung der Skelettlinie 


dys\ — Vrjt+ r+ ust jEy 
(a: ) baa Urj + ujyr+ ujs : (5.18) 
Hieraus erhalt man durch Integrieren iiber die x-Achse fiir die Skelettlinie 
=H] 
é ee Vrj t+ yr + YjJ + Ey 7) 
Ysj ae i Ur; i ujR ze ujJ dx; . (5.19) 


—Ij/2 


Zur Ermittlung der Dickenverteilung der Schaufelkonturen wird ein krummliniges Ko- 
ordinatensystem s und r langs der Skelettlinie bzw. senkrecht dazu eingefiihrt (vgl. Abb. 21). 
/ Nach der Kontinuitatsgleichung muB die von der 
% Schaufelnase bis zur Stelle s auf Grund der Quell- 
at belegung q(s) ausgetretene Fliissigkeitsmenge inner- 
Bb halb der Kontur bleiben. Daraus ergibt sich als Be- 
il dingung fiir die Dickenkoordinate die Gleichung 


i Oey 
1 


s | qj(s)) ds; = | Ws (spr;) dr;. (6.20) 
0 0 


Abb. 21. Zur Berechnung der Schaufelkontur. 


| Ersetzt man W, (s;,r) naherungsweise durch die Geschwindigkeit auf der Skelettlinie W, j{Sj)> 
© so erhalt man 


Ss. ej 
1 V2 , 1 TA f 
FI Was. | 98) 45) = ape, i qj(xj) 4x; « (5.21) 
0 —1j/2 


i Hierbei ist die Skelettgeschwindigkeit Ws ;(s;) aus dem induzierten Geschwindigkeitsfeld der voll- 
)§) standigen Singularitatenbelegung so anzuschreiben: 


Ws (sj) = Wg,;,(s;) + Ujzq(S;) 
1 ; ee: 
E ji+¥2 Upj + Ur t Ys t+ Ujnq(%}) + ¥sj Voz + Rn + Yt 2), ; (5.22) 
sf 

» wobei W,;,(s;) die Skelettgeschwindigkeit fiir den Fall bedeutet, bei dem die betrachtete Schaufel 
.@ nur mit einer Wirbelbelegung versehen ist. 
_ Die Geschwindigkeitsverteilung auf der Schaufelkontur Wx, wird auf gleiche Weise wie beim 
a) Hinzelgitter! aus dem gesamten induzierten Geschwindigkeitsfeld abgeleitet. Die Konturgeschwin- 


Ws; UjE 
Wepre IPL Y j 1). 5.23 
BLS Vis "Ki 2 Wsj Vl + y¥3% oe 


s) Aus der Konturgeschwindigkeit kann man leicht nach der Bernoullischen Gleichung die statische 
Druckverteilung px auf der Kontur erhalten: 


tts a ie a tal (5.24) 
“ We a ( a 


‘a, Aus den hier angegebenen Gleichungen sind die Schaufelkonturen und Konturgeschwindigkeits- 
yi -verteilungen der beiden Gitter fiir gegebene Ausgangsdaten zu ermitteln. Hierbei erfordert natur- 
ot) gemaB der Rechenschritt fiir die Ermittlung der induzierten Geschwindigkeit von der anderen 
© Schaufelreihe den gréBten Aufwand im ganzen Rechenverlaut. 

b) Berechnungsbeispiele. Fiir die Berechnungsbeispiele an Tandem-Pumpengittern sind 
‘3 die Schaufelkonturen, Geschwindigkeitsverteilungen sowie daran anschlieBend die zugehérigen 
8) Verlustbeiwerte fiir die in der Tabelle 1 angegebenen Auslegungsdaten ermittelt worden. Das 


1 Siehe FuBnote 1 von Seite 215. 
Gye 
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optimale Belastungsverhaltnis ist nach Abb. 3 als (AW,,/AW,)* = 90,5 angenommen. Dadurch 
wird die Zwischenstrémung und die Zirkulation festgelegt: 


8, =120,8°, ,/1W,=Iy/l W, = 0,596 . 


In diesem Fall stimmt die Translationsstromung W, mit der Zwischenstrémung W, tiberein. Die 
x-Richtung der beiden Koordinatensysteme wird parallel zur Grundstrémung W..; des Gitters I 
bzw. W.,, des Gitters IT festgelegt. Um eine giinstige Druckverteilung um die Schaufel zu gewahr- 
leisten, wurden die Wirbel- und Quellbelegungen wie in 

Abb. 22 ausgewahlt. 
| Die Rechnungen sind bei drei Anordnungen durch- 
| | gefiihrt worden, namlich bei Abstand a/l = 0,1 und 


15 7 ie 


Versetzung h/t =0, bei Abstand a/l = 0,1 und Ver- 
setzung h/t = 0,5 und bei unendlich groBem Abstand 
a/l = co. Die letzte Anordnung ohne Interferenz dient 
hierbei als Vergleichsbeispiel fir den InterferenzeinfluB. 
; Die Rechenergebnisse zeigen, daB die Anderung der 

‘eee = Skelettlinie mit der Versetzung schon bei verhaltnismaBig 
ie kleinem Abstand a/! = 0,1 recht geringfiigig ist, und zwar 
| andert sich der Schaufelwinkel mit der Versetzung bei 
f 4) Gitter I héchstens um 0,1°, bei Gitter IT um 0,5°. Auch 
I | die Anderung der Dickenverteilung der beiden Gitter 


| ist dabei nicht nennenswert. Aus diesen rechnerischen 
O5,—~* + 


2, af 


40r. 


Tabelle 1, Auslegungsdaten der Rechenbeispiele fiir den Entwurf 


eines T'andem- Pumpengitters 


= ale Bx 140° Ay 41,8° 
{ B, 120,8° Any 16,6° 
\ B; 90° 6 25,22 


AB 50° all 0,10 
ABy 19,2° bo/l 0,952 


aS Abu 30,8° 2 ie = 0,567 
5 


‘ hy [h 
ie Ly Boy 131,8° le = 0, 0,067 


Ons 0 05 Boy 106,6° Lyy/t 1,00 


0) 
Abb, 22. Wirbelbelegung y und Quellbelegung q langs rt Wo 1,192 dmax|t 0,10 
der Schaufelsehne bei den Beispielrechnungen. I, n/l Wa 0,596 Wa Syl m/s 


Ergebnissen kann man wohl die Folgerung ziehen, daB der EinfluB der gegenseitigen Interferenz 
auf die Schaufelkonturen schon bei kleinem Abstand a/l = 0,1 so gering ist, daB man beim Entwurf 
eines ‘Tandem-Pumpengitters kaum darauf Riicksicht zu nehmen braucht, da die Abweichung der 
Konturen von denjenigen ohne Interferenz nur so gro ist wie die Genauigkeit bei der gewohnlichen 
Schaufelherstellung. 

Abb. 23 zeigt die Konturgeschwindigkeitsverteilungen beider Gitter. Bei Gitter I ist die An- 
derung der Konturgeschwindigkeit in der Nahe der Hinterkante ziemlich groB, dagegen andert sich 
die Konturgeschwindigkeit des Gitters II in der Nahe der Vorderkante recht wenig. Es ist nahe- 
liegend, da die Abweichung der Konturgeschwindigkeiten nur in dem hinteren Teil des Gitters I 
und in dem vorderen Teil des Gitters II entsteht, weil diese Teile der anderen Schaufelreihe be- 
sonders nahe stehen und am starksten in deren EinfluBgebiet liegen. 

Anhand der oben ermittelten Konturgeschwindigkeiten wurden die Verlustbeiwerte der beiden 
Gitter nach dem grenzschichttheoretischen Verfahren von L. Speidel und N. Scholz! unter der An- 
nahme berechnet, daf die Grenzschicht vollturbulent ist und die Reynoldszahl Re, = Wy = 1,7: 1a 
fiir Gitter I, Re, = Wl/y =1,5- 10° fir Gitter IT betragt. Die Ergebnisse sind in Abb. 24 dargestellt. 
Bei Gitter I andert sich der Verlustbeiwert periodisch mit der Versetzung und weicht maximal um 1,5% 
vom Mittelwert ab. Der Minimalverlust tritt bei Versetzung h/t = 0,5 auf und der Maximalverlust 


1 Siehe FuBnote 1 von Seite 202. 
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Abb. 23. Ermittelte Konturgeschwindigkeitsverteilung. — a) Auf der Schaufelkontur des Gitters I, b) Auf der Schaufelkontur des Gitters IL 
bei Abstand a/l = 0,1 und Versetzung h/t = 0; 

bei Abstand a/! = 0,1 und Versetzung h/t = 0,5; 

bei Abstand a/l = co, ohne Interferenz. 


Die ermittelten Schaufelkonturen bei a/] = 0,1 und Versetzung h/t = 0 sind mit eingezeichnet, s) Skelettlinie, k) von Skelettlinie und 
Dickenverteilung aufgebaute Schaufelkontur. 


cy 


O15 
Crt bva)y 
0,10 
005 
Versetzung ® 
, 0 05 10 


Abb. 24, Anderung der gerechneten Verlustbeiwerte des Tandem-Pumpengitters mit der Versetzung 
a: l/y = 1,7- 105, Reyy = Wy-l/y = 1,5 - 10°. 
Verlustbeiwerte bei Abstand a/I = 0,1; 
Verlustbeiwerte bei Abstand a/]l = 0. 


Rey = W. 
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bei Versetzung h/t = 0. Der Mittelwert der Verluste stimmt mit dem Verlustbeiwert ohne Inter- 
ferenz gut iiberein. Die Verlustbeiwerte des Gitters II andern sich demgegeniiber kaum, weil die 
Anderung der Konturgeschwindigkeit vernachlassigbar gering ist. Der resultierende Verlustbeiwert 
von den beiden Gittern andert sich ebenfalls periodisch mit Abweichung von 5% vom Mittelwert, 
und die Phase der Schwankung stimmt mit derjenigen des Gitters I iiberein. 

Diese rechnerischen Ergebnisse beweisen deutlich, daB die Schwankung der Verluste des Git- 
ters I im wesentlichen durch eine potential-theoretische Betrachtung erfaBt werden kann, weil die 
rechnerische Anderung der Verluste des Gitters I mit der Versetzung hinsichtlich der Starke und 
der Phase der Schwankung mit den MeBergebnissen sehr gut tibereinstimmt. Demgegeniiber kann 
die in Wirklichkeit auftretende starke Anderung der Verluste des Gitters II nicht mit der Poten- 
tialtheorie, sondern nur mit einer grenzschicht-theoretischen Betrachtung erklart werden. 

Nach den hier gewonnenen Ergebnissen kann man im allgemeinen beim Entwurf der Schaufel- 
konturen des Tandem-Pumpengitters mit maBigem Abstand (a/! > 0,10) die Gitter mit guter 
Naherung als zwei an sich selbstandige Gitter betrachten. Dariiberhinaus kann man durch Auswahl 
der giinstigsten Versetzung fir den Verlust ein Minimum erhalten. 


6. Interferenzeinflu8 einer Delle in der Zustrémung auf die Verluste der nachstehenden Schaufel. 
Um den bei den Messungen an Tandem-Pumpengittern deutlich beobachteten starken Interferenz- 
einflu8 der Nachlaufdellen vom Gitter I auf die Schaufelverluste des Gitters II noch naher zu unter- 
suchen, wurde ein einfacher Modellversuch durchgefiihrt. Um den Vorgang dieses Einflusses der 
Delle zu verstehen, kénnte man sich zunachst vorstellen, daB eine Delle in der Zustrémung die 
Druckverteilung der nachstehenden Schaufel beeinfluBt, und daB sich daraus eine gewisse Anderung 
der Verluste ergibt. Die Ergebnisse der Druckverteilungsmessungen zeigen jedoch, daB sich 
die Druckverteilung durch eine Delle sehr geringfiigig andert und dadurch keine groBe Anderung der 
Verluste zu erwarten ist. Infolgedessen muB eine Delle in der Zustrémung unmittelbar den Vorgang 
der Grenzschichtbildung, vor allem die Ablésungspunkte an der Schaufel, stark beeinflussen und 
die entsprechende Anderung der Verluste zur Folge haben. 

Die Strémung in einer Delle enthalt gewiB eine kleinere kinetische Energie als die Strémung 
auBerhalb der Delle. Deshalb wird beim Uberwinden des von der AuBenstrémung aufgepragten 
Druckanstieges die Strémung in der Delle eher als die ungestérte Strémung ihre kinetische Energie 
erschépfen und infolgedessen eher ablésen, so daB daraus ein gewisser Zuwachs der Verluste folgt. 
Dieser Vorgang wird hierbei Interferenzeffekt der Delle genannt. 


MeBebene: Z 


tx 
ALS. 


Abb. 25. Versuchsanordnung fiir die experimentelle Bestimmung des Interferenzeinflusses einer Delle in der Zustrémung auf die Verluste 
der nachstehenden Schaufel. — J Reduzierdiise, 2 Einlaufteil, 3 Diffusorteil, 4 Austrittsteil, 5 Storkérper, verstellbar und austauschbar, 
6 MeBschaufel NACA 0010, 7 MeBsonde fir den statischen Druck Py 8 Nachlaufsonde fiir den Gesamtdruck & 4 und statischen Druck Py: 


a — Abstand zwischen dem Stérkérper und der Mefschaufel, 
h — seitliche Versetzung des Stérkérpers. 


Andererseits hat die Strémung in einer Delle eine kleinere mittlere Geschwindigkeit und einen 
entsprechend kleineren mittleren Staudruck als die ungestérte Strémung. Weil sich der Reibungs- 
und Druckwiderstand einer Schaufel in einem ahnlichen Strémungsfeld proportional zum Staudruck 
andert, kann eine Delle direkt vor der Schaufel die Verminderung des Widerstandes bzw. des Ver- 
lustes verursachen. Dies soll als Schatteneffekt der Delle bezeichnet werden!. 


_1 P. Ruden, Windkanalmessungen uber den Windschatteneinflu8 auf Rechtecktragfliigel mit symme- 
trischem Profil, Jahrbuch Deutsch. Luftfahrtforsch. 1940, S. I 204. 
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In Wirklichkeit treten beide Effekte gleichzeitig auf, und zwar in einem bestimmten Verhiiltnis, 
welches von der Starke des Druckanstieges langs der Schaufel und von der Anordnung der auf- 
treffenden Delle abhangig ist?. 

Die Versuchsanordnung ist in Abb. 25 dargestellt. Es wurde in vereinfachter Weise eine Einzel- 
schaufel vom Profil NACA 0010 in einen Diffusor mit 8° Offnungswinkel eingebaut, dessen Druck- 
anstieg etwa demjenigen eines gewohnlichen Pumpengitters entspricht. Weil das Problem nur den 
Vorgang der Grenzschichtbildung auf der Schaufelkontur betrifft, geniigt es, den Druckgradienten 
langs der Schaufel in dieser Weise festzulegen, um den EinfluB der Delle zu untersuchen. Bei den 
Untersuchungen wurden die gewiinschten Dellen in der Zustrémung mit Hilfe von zwei verschie- 
denen Stérkérpern erzeugt. Beide Stérkérper sind leicht austauschbar und kénnen parallel zur 
Schaufel in verschiedenen Abstanden und seitlichen Versetzungen angeschraubt werden. 

Die Nachlaufmessungen sind zunachst in MeBebene V (Vorderkante der MeBschaufel) ohne 
MeBschaufel nur mit den Stérkérpern durchgefiihrt worden, um zu bestimmen, wie die Nachlaufe 
der Stérkérper in der Vorderkantenebene aussehen. Anhand dieser Nachlaufmessungen kann man 
die Geometrie der Dellen in Abhangigkeit von der Anordnung der Stérkérper bestimmen und gleich- 
zeitig die Verlustbeiwerte der Stérkérper allein ermitteln. Dann sind die systematischen Nachlauf- 
messungen in MeBebene A (halbe Schaufeltiefe hinter der MeSschaufel) bei verschiedenen Anord- 
nungen der Stérkérper durchgefiihrt worden. 

Der Verlustbeiwert ¢y ist hierbei als Energieverlust nach der Gleichung 


pa Energieverlust Wy) gz—asly) dy (6.1) 
Y ~~ Strahlenergie in Z-Ebene W, 0 we i] 
pee 


in Delle 


definiert. Diese Definition hat den Vorteil, daB sie auch bei der verzégerten Grundstrémung fir 
einen bestimmten Energieverlust unabhangig von der MeBebene immer den gleichen Verlust- 
beiwert liefert. Weil der Nachlauf in MeBebene A die beiden Verlustdellen vom Stérkérper und 
von der MeBschaufel enthalt, ergeben sich die Verlustbeiwerte der MeBschaufel als Differenz der 
Gesamtverluste und der Stérkérperverluste nach der Formel 


rf Way) &z—8Aay) ty Wry) 8&z—srly) dy (6.2) 
cys Wz te | A Wz 0 we Le , 
yy We “oy Z 


Abb. 26. Charakteristische Parameter einer Nachlaufdelle und Zusammenhang mit der MeBschaufel. 
1) Stérk6rper, 2) MeBschaufel, 3) MeBebene V, 4) Nachlaufdelle des Stérkérpers ohne Mefschaufel, 5) Einflufdreieck der auftreffenden Delle. 
WwW. . |W: Tiefe der Delle; 


cigs a Charakteristische Breite der Delle; 


: Fliche des EinfluBdreiecks mit Grundlinie d 


° x : 
: von der charakteristischen Breite B bedeckte Flaiche des Dreiecks. 


yao 


Bei allen Messungen wurde die Zustrémgeschwindigkeit Wz, = 44 m/s eingehalten. Die auf die 
Zustrémgeschwindigkeit bezogene Reynoldszahl der Schaufel betrug daher 


Ree 5.8108: 
v 
In diesen Untersuchungen laBt sich die Geometrie der auftreffenden Dellen durch ihre Tiefe 
und eine charakteristische Breite darstellen (vgl. Abb. 26). Die Tiefe der Dellen ist hierbei dem 


1 K. Bursig, Experimentelle und theoretische Untersuchungen iiber die Stromung um ein Tragfliigelprofil 
bei Druckabfall und Druckanstieg und bei ungleichformiger Zustrémung, Diplom-Arbeit Nr. 53 des Instituts 
fiir Strémungsmechanik der Techn. Hochsch. Braunschweig (1958). 
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folgenden Verhaltnis entsprechend quantitativ festgelegt: 
Minimalgeschwindigkeit der Delle 


Zustromgeschwindigkeit _ We 


Die charakteristische Breite der Delle B 1a8t sich als die Breite eines angenommenen rechteckigen 
Loches in der Zustrémgeschwindigkeit definieren, welches die gleiche Tiefe und die gleiche Ver- 


StorkorperI — Storkorper Ir 
Stromungsrichtung 


Abb. 27. Anderung der charakteristischen Parameter der Dellen mit dem Abstand a. 
a) Anderung der Tiefe der Dellen bei Stérkérper I und II, b) Anderung der charakteristischen Breite der Dellen bei Stérkérper I und II. 
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Abb. 28. Anderung der Verlustbeiwerte mit der seitlichen 
Versetzung bei verschiedenen Abstinden fiir Stérkérper I. 


Re = W, -l/v = 5,8 - 108. 
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Abb. 29. Anderung der Verlustbeiwerte mit der seitlichen 
Versetzung bei verschiedenen Abstinden fiir Stérkérper II. 


Re = -I/y = 5,8 - 108 
e = W-l/v = 5,8 - 108, 
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drangungsdicke wie die urspriingliche Delle besitzt. Die charakteristische Breite ergibt sich dann 
nach der Formel 


1 


B=— 2 
Vy Wy min 


{ [W, — WAy)] dy. (6.3) 
in Delle 


Dabei bedeutet Wy die ungestérte Zustrémgeschwindigkeit auBerhalb der Delle in MeBebene V. 
Die Anderung der Tiefe und der charakteristischen Breite der Dellen ist in Abhangigkeit von der 
Anordnung der Stérkérper in Abb. 27a und b dargestellt. 

Aus den systematischen Nachlaufmessungen in MeBebene A wurden die Verlustbeiwerte er- 
mittelt, die die Verluste der MeBschaufel Cj, und des Stérkérpers Cys, enthalten. Sie wurden 
fiir die Stérkérper I und II in Abb. 28 bzw. 29 iiber der seitlichen Versetzung bei verschiedenen 
Abstanden der Stérkérper aufgetragen. Die ausgezogenen geraden Linien geben dabei die Verlust- 
beiwerte der MeBschaufel bzw. der Stérkérper wieder, die ohne Interferenzeinflu8 einzeln gemessen 
wurden. Sie werden fiir die Messungen mit Interferenz zum Vergleich herangezogen. 

Wie bei Tandem-Pumpengittern nehmen die Schaufelverluste auch hier unter dem EinfluB einer 
Delle in der Zustrémung stark zu. Es ist daraus deutlich zu erkennen, daB eine die Schaufelmitte 
treffende Delle in der Zustrémung die Schaufelverluste um etwa 35 bis 40°% vergroéBert, wobei die 


; Verluste kaum von der Tiefe und Breite der Delle abhangen. Der Interferenzeinflu8 einer Delle 


vermindert sich allmahlich mit zunehmender seitlicher Versetzung; und wie man sich leicht vor- 
stellen kann, verschwindet der InterferenzeinfluB einer Delle mit einer schmaleren charakteri- 
stischen Breite bei kleinerer seitlicher Versetzung. 

Der Interferenzeinflu8 einer Delle ist praktisch abgeklungen, wenn die Delle so weit versetzt 
wird, daB die charakteristische Breite der Delle sich auBerhalb der Maximaldicke der Schaufel befindet. 

Die Verlustkurve bei Abstand a/l = 0,05 mit Stérkérper II in Abb. 29 zeigt zwischen h/d,,,, 
= 0,50 und 0,75 einen starken negativen InterferenzeinfluB der Delle. Das beudetet aber nicht, 
daB eine Delle bei dieser Anordnung einen giinstigen Einflu8 auf die Schaufelverluste ‘ausiibt, 
sondern nur, daB der Stérkérper dabei auf die MeSschaufel wie ein Vorfliigel wirkt und daB der 
Verlust des Stérkérpers selbst wegen der beschleunigten Strémung durch den Spalt abnimmt. 
Daraus ergibt sich eine Verminderung des Gesamtverlustes. 

Der Interferenzeinflu8 hangt von drei Parametern ab, namlich von der Tiefe der Delle, der 
Breite der Delle und der seitlichen Versetzung. Das verursacht die Schwierigkeit, den EinfluB 
dieser drei Parameter in einem Diagramm darzustellen. Dazu ist ein Versuch unternommen wor- 
den, aus der charakteristischen Breite und der seitlichen Versetzung einen einzigen neuen Parameter 


06 


"06 G7 
Abb. 30. Relative Anderung der Schaufelverluste in Abhangigkeit von der Tiefe qy min/4z = (W,, min| 2)" 


und der Uberdeckungsziffer K der auftreffenden Delle. 


—P. 
4 = 0,195;  Reynoldszahl: Re = W, «Uy = 5,8 . 108. 


@ 2 
— W 
2 Z 


Druckanstieg des Diffusors: 


zu bilden und die Ergebnisse in Abhangigkeit von nur zwei Parametern in einem allgemeingiiltigen 
Diagramm darzustellen. Ein neuer Parameter — Uberdeckungsziffer K — wird mit Hilfe der 
Darstellung in Abb. 26 durch die Definition K = F/F, festgelegt. Samtliche MeBergebnisse sind 
in Abb. 30 dargestellt, und zwar sind die relativen Anderungen der Verluste hierbei als Linien 
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konstanter Interferenzeinfliisse in Abhangigkeit von der Tiefe qymin/¢z = (Wy min/Wz)? und der 
Uberdeckungsziffer K eingetragen. Bei jeder Tiefe der Delle nehmen die Schaufelverluste mit der 
Uberdeckungsziffer zu, und der Einflu8 der Uberdeckungsziffer auf das Anwachsen der Verluste ist 
um so groRer, je tiefer die auftreffende Delle ist. 

Um den in wirklichen Tandem-Pumpengittern auftretenden InterferenzeinfluB der Delle quanti- 
tativ analysieren zu kénnen, muB man solche Untersuchungen bei verschiedenen Druckanstiegen langs 
der Schaufel durchfiihren, da der Druckgradient einer Schaufel im Gitterverband auf der Saug- und 
Druckseite verschieden ist. Das KinfluBdiagramm in Abb. 30 gilt nur fiir den gemessenen Druckanstieg. 

Parallel zu den experimentellen Untersuchungen ist auf verschiedene Weise versucht worden, 
den Interferenzeinflu8 mit einer grenzschicht-theoretischen Betrachtung zu erklaren. Diese theore- 
tischen Behandlungen haben aber alle ergeben, daB der Schatteneffekt viel starker als der Inter- 
ferenzeffekt auftritt und die resultierenden Schaufelverluste unter dem Kinflu8 einer Delle stark 
absinken, wahrend die gemessenen Verluste stets um etwa 35 bis 40°% gréRer sind. 

Dieser groBe Unterschied zwischen der Messung und der Theorie ist vermutlich darauf zuriick- 
zufiihren, daB ein Stérkérper nicht nur ein Nachlaufprofil hinter sich erzeugt, sondern auch gleich- 
zeitig sehr ausgepragte Turbulenz oder sogar eine Art intermittierende Stérung der Strémung wie 
in einer WirbelstraBe in seinem Nachlauf hervorruft. Solche starke Turbulenz oder Stérung kann 
natiirlich den Vorgang der Grenzschichtbildung stark beeinflussen und einen groBen Interferenz- 
effekt zur Folge haben, der von der Theorie bisher nicht erfaBt werden kann. 


7. Zusammenfassung. Es wurden die Strémungsvorginge in ebenen Tandem-Pumpengittern 
mit starker Umlenkung theoretisch und experimentell untersucht. 

Zunachst wurden ohne Beriicksichtigung der gegenseitigen Interferenz optimale Tandem- 
Gitteranordnungen fiir beliebig vorgegebene Geschwindigkeitsdreiecke untersucht, und zwar auf 
Grund von Ergebnissen an Einzelschaufelgittern mit Profilen NACA 0010, 4410 und 8410. Die 
Ergebnisse sind in Diagrammen dargestellt, die bei der Entwurfsaufgabe in der Praxis unmittelbar 
benutzt werden kénnen. Es ergibt sich weiterhin, daB in dem Bereich cy l/t > 1,2 das Tandem- 
Pumpengitter geringere Verluste liefert als das entsprechende Einzelgitter. 

Der Einflu8B der gegenseitigen Interferenz wurde an einem Tandem-Pumpengitter aus Profilen 
NACA 8410, das ohne InterferenzeinfluB eine Umlenkung von 50° beim Zustrémwinkel 6, = 140° 
verwirklicht, experimentell untersucht; Nachlauf- und Druckverteilungsmessungen wurden bei 
verschiedenen Abstanden und Versetzungen der beiden Gitter zueinander durchgefiihrt. 

Die MeBergebnisse zeigen, da sich die Verluste des Gitters I sinusférmig mit der Versetzung 
andern, was auf einen rein potentialtheoretischen Interferenzeinflu8 hindeutet, wahrend die Ver- 
luste des Gitters II etwa konstant bleiben und nur dann stark zunehmen, wenn sich die Schaufeln 
des Gitters II in den Dellen vom Gitter I befinden, was auf eine Interferenz der Grenzschichten 
beider Gitter schlieBen laBt. 

Beim iblichen Abstand von a/l = 0,15 bis 0,50 ist die Versetzung dann am giinstigsten, wenn 
sich die beiden Schaufeln wie bei einem Spaltfliigel mit einer Spaltweite von 1/, der Teilung an- 
ordnen, erstens weil der Gesamtverlust dabei am kleinsten ist (etwa 10% Kleiner als ohne Interferenz) 
und zweitens weil der Abstrémwinkel mit demjenigen ohne Interferenz sehr gut tibereinstimmt und 
daher keine Korrektur fiir den Abstrémwinkel nétig ist. Sind die optimalen Teilungsverhialtnisse 
der beiden Gitter verschieden, sollen die beiden Schaufeltiefen so bestimmt werden, da®B die beiden 
Gitter die gleiche Teilung besitzen. 

Um optimale Anordnungen von Tandem-Pumpengittern unter Beriicksichtigung der gegen- 
seitigen Interferenz theoretisch bestimmen zu kénnen, wurde ferner ein potentialtheoretisches 
Rechenverfahren zum Entwurf eines Tandem-Pumpengitters entwickelt (erste Hauptaufgabe). 
Die Rechenbeispiele ergeben, daB die Strémungsvorgange am Gitter I gut mit der Theorie erfaBbt 
werden kénnen, wahrend die bei Messungen beobachtete starke Anderung der Verluste des Gitters II 
nicht damit zu erklaren ist. 

Es wurde daher eingehend der InterferenzeinfluB einer Delle in der Zustrémung auf die Verluste 
der nachstehenden Schaufel, der fiir die Anderung der Verluste des Gitters II eine entscheidende 
Rolle spielt, an einem vereinfachten Modell untersucht. Durch eine Delle in der Zustrémung kénnen 
die Verluste an der Schaufel in einer verzégerten Strémung je nach der Geometrie und der Lage der 
auftreffenden Delle relativ zur Schaufel stark zunehmen. Bei maSigem Druckanstieg nehmen die 
Schaufelverluste durch eine Delle, die gerade die Schaufelmitte trifft, um etwa 35 bis 40% zu. 


(Eingegangen am 2. September 1958.) 
Anschrift des Verfassers: Dr.-Ing. Hideo Ohashi, Ichikawashi (Japan), Suganomachi 2—529, 
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A note on thermal stresses in hollow cylinders* 
By M. M. Stanisic and R. M. McKinley 


1. Introduction. R. Trostel' has recently given a solution for the stress field in hollow cylinders 
with axial and rotational symmetry subjected to a steady-state temperature field, and such that 
the axial strain component can be assumed constant. His solutions allow for a variation in the 
physical properties of the media, but are restricted to the case in which the Poisson ratio is one- 
half (y = 1/2). Approximate extensions of these solutions for general y are given in this paper which 
differ from the exact solutions by less than 3% in magnitude. Where higher accuracy is desired, 
a numerical method illustrated in the paper may be employed. 


2. Nomenclature. We denote 
T = T(r) = Temperature, °C, 
$i 


E = Young’s Modulus, kg/cm?, 
Ti = component, in the 1-coordinate direction, of the stress vector acting on a surface whose 
normal is in the j-coordinate direction, kg/cm?, 

Ce; = unit extension in the i-coordinate direction of an element initially parallel to the i- 
} coordinate direction, dimensionless, 
ay = Poisson’s ratio, dimensionless, 

o(T) = coefficient of thermal expansion, (°C)-}, 

To = external cylinder radius, cm, 

r; = internal radius of cylinder, cm, 

Pi Po = pressures imposed on the internal and external surfaces of the hollow cylinder respec- 


tively, kg/cm?, 


91(r) 
V(r) = integrals defined by equations (18)—(20), 
alr) 
y(r) = stress function (not in the Airy sense) defined by equation (14), 
T;, Tj = internal and external surface temperatures, respectively, of the hollow cylinder, °C, 
D = bine. the differentiation operator, 
dr 
6 =y (r + S| —y (r oy = central difference operator, 
6? = Venlo Va yn > 
h = space increment in the independent variable, 
€ = the shifting operator, 
ep(r) =p(r +h), 
A = the forward difference operator, 


Ap(r) =p(r+h)—p(), 
V = the backward difference operator, 
y Pir) = p(r)— p(r—h), 


= thermal conductivity, cal/sec — em? — (°C/cm). 


3. Approximate solutions of the Problem. The generalized Hooke’s law for isotropic media whose 
properties may vary with temperature gives the following relation between the normal elastic 
stresses, thermal stresses, and the corresponding strains (referred to a cylindrical coordinate system) : 


du 1 Sooke 
dr = 2, = E(T) [t, ae (Too a te)| Ele «(0) dO ? (1) 
u 1 ae 
rs = 6965 E(T) [too meal Css ee t,7)] is Up a(@) dO ’ (2) 
1 6=T P 
7 id E(T) [tg inns (tr a Tee)] cote, a() dQ . ( ) 


* M. M. Stanisié, Lectures in Mathematical Elasticity during summer semester 1958, Purdue University. 
1 R. Trostel, Ing.-Arch. 26 (1958) S. 134. 
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The conditions for static equilibrium reduce in the case of axial and rotational symmetry and 
zero body forces to 


d 1 
art re r (;, Tee) = 0. (4) 


The one independent equation of compatibility follows directly from the left-hand sides of 
equations (1) and (2): 
dego _ du u 


dpe dre er 


hence 
d 
r pe + (99 —@r) = 9. (5) 


It will be assumed in this paper that the axial component of strain is constant 
C52 — A, 


in which case equation (3) gives an expression for the axial stress component: 
@=T 
Baa Mtr a ere) aerial tau 0 al cee (6) 


If the ends of the cylinder are constrained from moving e,, = 0 = A, whereas, if the ends are free 
and are acted on by a force P, the constant A, is obtained from the condition of equilibrium 


(20rt ar =. (7) 
Trostel takes y = 1/2 and operates directly on equations (1), (2), (3), (4), and (6) to obtain the 
radial component of the displacement vector and the stress field. The resulting solutions are, 
however, not applicable to the more common cases in which » lies in the range of 1/3 to 1/4. Since 
the general problem appears not to have a solution in terms of quadratures alone, it is the purpose 
of this paper to determine corrections factors which will allow the use of Trostel’s equations with 
the same order of accuracy as usually found in experimental data for the variation of physical 
properties with temperature. For this purpose, it is more convenient to work through the com- 
patibility equation (5). Using equation (6) to eliminate t,, from equations (1) and (2) gives 


0=T 
ere = am [(1—>*) t,, —» (1 +9) to9 —9 A, E(T)] + (1 +) f o(0) dO, (8) 
E(T) a6 
6@=T 
em ea Il — cep (Lt). 9 ART) ey eae den (9) 
E(T) @=0 
1 , 
90 — err = iT (Coe tr) : (10) 
Noting that E(T) = E[T(r)], equation (9) then gives 
de 1 dE — y? 
SS E2 dr [1 y?) Tog _? (1 ai hed == : E : Gee = = ”) ae 
+ (1 +9) 0(T) 2 (11) 


From equation (4): 


ater 
r Via ir (Too a T,) 


80 that 
dege pall r dE r dt 0) 1 
= — Beg (1) 99 — 9 (1 +9) te] + A — 9) 5 22 — 9 (1 +9) 5 (99 —1,,) 


dT 
+ (1 +9) raf T) (12) 
Substituting equations (10) and (12) into the compatibility equation (5) 


dE v dte@e 1 
mer dC rear Erk ar + E (tog —T,,) + 


l—y 


rE%(T)“*=0. (13) 
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Suppose that there exists a function y(r) of at least class C2 such that 


Ty = | (14) 
then from equation (4) 
d 
Tee = + (15) 
and 
dtem —s d*y 
dr dr® 
follows 


dp = 
secs to) oa Sees ET 


Substituting these expressions into equation (13) and collecting terms we have 


dy rn dy v ; Sas 1 dT 
(CE) Gt Er EVE + (Ayre —E\P = par a(S (16) 
where 
de 
fe cairns 


It may be varified that equation (16) has no solution by simple quadratures by assuming that 
E(r) may be expressed as polynomials in r and examining the Riccati equation corresponding to (16). 
For the case vy = 1/2, equation (16) reduces to 


@ oa dT 
(E) Se +(E—rE) 24+ (E—k£E)t=— r Eo(T) (17) 
A particular solution of the reduced equation 
d*y 1 dy / —_ 


is seen by inspection to be yy =r. 
The general solution of equation (17) is then obtained by making the substitution 


V(r) = Yo V(r) = r Yr) 


which gives 
- d’*p 3 Ee AO es E dT 
at Farle ee Nice 
This is of the first order in dg/dr and may be solved directly: 
spe Ee ae aT E 
ee o(T) —- dr + ¢ =a 


by a second integration 
Bisa ead E 
Q Soe a(T) ar) dr + auf rar + ¢, 
where c, and c, are constants of integration. Finally we have 


ae ae 
E 
ya —2r[ E(fra(T) Farr tar fdr ter. 


From equations (14) and (15) follows 


E dT E 
tata —2[ a([rtatT) Gar) dr + a a Or + Cy, 
d E dT 10; dT 
too = Ge = —2[ a([ alt) G arlar—2 zs | rtalT) Gat 


E E 
+c adr + eG ze: 
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If the constants of integration are determined subject to the boundary conditions 


ais Upp iitee ae aD 
LHI Ly Rise ooaeall 
and the following notation, introduced by Trostel, is employed: 
eat)= | ar, (18) 
E aT 
727) == |" (T) a 
@=T @=T (19) 
pr) =r | 0(6)40—2+( | x(0) 10) ar, 
e=0 o=0 
9,(r) = lhe (rar) a or 
i (20) 


= [ Lean ar, 


one finally obtains for the stress field corresponding to vy = 1/2 


Ter = —Pi— (Po— P(A) + Beale) (Fey te (21) 


P1("o) Piro) Palo) 
E E 
= +r) J+ (= +94(0) 
= —e, ee Tees r __ P(r) ae p3(r) 22 
Teo = Pi (Po ra| V1 (ro ) A 2 Pa(To) Pu(ro) P3(To) ( ) 
O= 
St eee eee f x(0) dO. (23) 
These are the solutions obtained by Trostel. 
Rather than simplifying the general equation (16), let 
y =a solution of the general equation (16), 


vc1/2) = a solution for the case y = 1/2, eqn. (17) 
and examine the differential equation for the difference: 


rE [y" — yaya] + [EB —7 E’ rE’ — E|~ —[r E’ — E]*a@ 


me 1—2» 2 
= alg l ae (24) 


where primes denote differentiation with respect to r. 
Let } = py — yo), then equation (24) suggests the approximation 


(7 E) i’ + (E—r Ef + (EE) Bc 


r E°o(T) ie (25) 
Equation (25) is solved in the same manner as equation (17) to give 


yee fal a(T) Far) dr + ay 7 [aa ter. (26) 


l—vy 
From the definition of h it follows that 
Cory = Ceri +, (27) 
d 
(Too)y = Too)1/2 = (28) 


Since t,, is specified on the boundary, the constants of integration a, and a, are determined by the 
conditions 


r= 7% b=0, 
a AS hia. 


7 
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On substituting into equations (27) and (28), the approximations to the stress field for general 


are obtained as 
alee) (Py gute) (29) 


Sigicher sees pil’) 
Gr, = Pi (Po me d/ 1 ar YP (ro) P3(r)o} ’ 


Pilro) 
ae col + 9,(7) i 
Soe r pale) + P,(7) (30) 
too = Pi (Po i) Ae Puro) Pi(To) P3(To) 


Tze = 0 (T, + Toe) + Ay ij E fio ; (31) 


is 


The error in these approximations will be eae in the section on appliaction, however, since 
y S 1/2, equation (31) shows that the error in T,, is no greater the larger error in either 7,, or To9- 
The radial component of the displacement vector is obtained from equations (2) and (9): 


@=T(r) 
u(r) = ep [A —¥*) too — (+9) tre] — 9 Ay + (LE) 7 | (0) dO. (32) 
o=0 


There are two limiting cases in which equations (29)—(32) are exact. The first case for y = 1/2 
is obvious from construction. The second case is that of constant physical properties (denoted by 
Q-subscript) where the following reductions occur 


[a] 1G) = 1, (T;— 7) 


(this temperature field is deduced in the first 
in(75) application given later), 


[b] walt) > Bo | jer lee cat 


r? (3 
ri 
Oo EL, dT OX Ho (LT, — T;) 1 
[c] Qo(r) > =a ° [r > dr = 4 = [T dr, 
in) 


r es fae 
X9 Ey (T) — Tj) dr __ % Eg (Ty — T;) = (=) 
[d] ale) > iB | _ : : zy 
O=T(r) I Fe 
le] f o(@) dO = a T(r) =a | Ty — (To— Ti) 


L 
;. 
00 in a) 
E ri 


Equations (29), (30), and (31) become 


ie) fee 
Tr, = — Pi — (Po — Pi) = oa ; ey , (33) 
0 u 
ce eee oe) 
Too = — Pi— (Po Pi) a mpg = . (22 ’ (34) 
0 0 t 
m(f) 
Ty, = ¥ (Trp + Too) + Ay Ey — % Eo | To (Tos Ty) za , (35) 
To 
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which are exact solutions!. The radial component of the displacement vector, equation (32), 


becomes * 
In ( ) 


r 


To 

u(r) = 2 [1 — 9%) too —¥ (1 +») tye] — 9 Ay + (1 +9) 97] Ty — (To— T:) —— 
0 apa 
‘ 1 Gl 

4. Numerical solution for the stress function. For a reliable estimate of the error in the approxi- 


mations given previously, it will be necessary to have a solution of the complete equation (16) for 
the stress function y(r). The equation is in standard form as given below 


" 1 [ne , E’ 1 1 1 dT 
(bE +(e —Hee—rhane eam 


(36) 


r if: 


Since equation (37) is linear, it is amendable to several methods of numerical solution®. One parti- 
cular method which takes advantage of the magnitude of the coefficients in equation (37) will be 
developed in this section. 

For convenience, introduce the following notation: 


A E? ; » 
RGlie | eee (38) 9 
‘ dP B\t ] Ev 1 
P(r) r= |(e)—2 =| (39) 
~ E’ 1 1 
OO ltr e alee (40) 
I dT 
and make the change of variable in equation (37) 
y =hy er 5 
Oa iay =a See 
poye? — emg (2) 
E\1ue2 
Pree). iS) 
which reduces it to the equation 
¥” + f(r) ¥ = g(r) (44) 


where 
f(r) = (4 Qr) —2 P(r) — [POR, (45) 


) 


g(r) = [etre], 


1/2 (46) 
8) =(z) FO) 
and with boundary conditions: 
Ge Wet 
aoa: ada , (47) 
E(r,.)\—1/2 
Tsoi: j= Pita ep) . (48) 


Since equation (44) is linear, if y, and y, are any two linearly independent solutions of the reduced 


equation “ii 
De ret) Saerke (49) 
and Yq is a particular solution of the complete equation (44), then the general solution of equation 
(44) is 
YH ai Tee Ve + Yo: (50) 
Since in the numerical solution one condition is known, namely y at r,, there is only one independent 
constant and the general solution (50) may be written as 


Vi GN co. (51) 
F i Be for example Chi-Teh Wang, Applied Elasticity, p. 75, New York 1953, who takes Pi = Po S03 
ah = ee 
* W. E. Milne, Numerical Solution of Differential Equations, Chapter 7, New York 1953. 
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The scheme for numerical calculations is then to determine Yo such that it satisfies equation (44) 
and the initial boundary condition (47). 
Then y, is determined such that it satisfies equation (49) and the initial boundary condition 


r=7: y=0, (52) 
then the constant c, is determined by the condition that 
¥(Fo) = 41 (7) + Yo(To) (53) 


must satisfy the final boundary condition (48). : 
To solve equation (44), it will be reduced to a second order difference equation. The following 
symbolic identity relating the differentiation operator D the central difference operator 6 and 
spacing increment of the independent variable, h = Ar, may be established!: 


Ge reas: 
ee sinh (3) t (54) 
1 12 2 32 
D=;(6 m3 +5051 ) 
il /? 124732 2 
ea g.31° + 94-519 


aL eo ee a ae esa ie ee 
Dt =73(8 Te°390°' gosto : 
ee St ee ay 

y =7(6 ay +598) (55) 

Operating on equation (55) with 6* gives 
== 1 6 
dy =7a(OY 19 x) 

whereas carrying out the same operation on equation (44) gives 


Oy” + 0 [f(r) y] = & [7] - 


These last two results provide the relation 
1 1 
6 [2(")] — 8 Lf) ¥] =7a(O'y 75° »| | (56) 


It is desirable to eliminate 64 y from equation (56) to increase the order of the process. Noting that 
equations (44) and (55) both give expressions for y’’, we may obtain an expression for 64 y 


1 1 
 [g(r) —f0) y] = By — Hoty + Oy 
, 


1 
fy =—12h g(r) 4+ 12h f(r)y+lPy SOY 
Substituting this expression into equation (56) gives 


i 
h? 6? [g(r)] — A? 62 [f(n) y] =— 12H g(r) + SO) y + REY +58 y— 75 Oy, 


— 62 [12 y + hfe) y] — 12 bef (r) y = — HE [g(0)] — 12 g(r) + 550%, 


2 f(r h2 2 p(r) 68 
(1 ieee ‘\y + fle) y = Wat ot BOs (57) 
where 
O21) —= Seat 2 Ba a bei (58) 
For convenience, denote ; 
ae L oe os y (59) 


1 F. B. Hildebrand, Introduction to Numerical Analysis, p. 136, New York 1956. 
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or 
bie rep , 
de 12 
M(t) = Bt g(r) +7520 | (60) 
NG) (61) 


ee EO, 
3 12 
Then equation (57) may be written more compactly as 


68 
6? Y + N(r) Y = M(r) ——. . 
240 
For a reasonable choice of space increment, h, the stress function is of such a nature that it may be 


quite accurately approximated by polynomials of degree four or less so that the truncation error 
in dropping the term 6° y/240 can be considered to be zero. Noting also that 


6? y= Yn41— 2 Yat be P) 
we finally arrive at the difference equation by which equation (44) is to be replaced: 
Yni1=M,(r) + [2 —N,(7)] Yn — Yr—1- (62) 
The difference equation corresponding to the reduced equation (49) is obtained by setting M(r) = 0 
in equation (62) 
Yi41 = [2—N,(r)] Y,— Yn—1- (63) 
Equations (62) and (63) afford a high accuracy solution for the variable Y and, by the inverse 


transformations, for the stress function yp. 
If however, we choose the space increment such that 


h= 0.10 


then we are usually justified in making considerable simplifications in these equations and still 
retain four significant figure accuracy. In the first place 


h? f(r) = 5 - 10-5 
so that to five significant figures 
Y=y. 
We then notice that ' 
N(r) = 5-10. 


If N(r) is dropped, the propagated error at the n' step of the solution is seen directly from equation 


(62) not to exceed 
en Sm [N(r) ¥] Max 
which after 100-steps amounts to app roximately 
¢,, = 25 - 10>? (y = 100)" 
Moreover, in most cases with four significant figures 
M(r) = hh? g(r). 
With these simplifications, equations (62) and (63) appear as 
Yos1 = h? g(r) +29 —Jn—1, (64) 
Vrpn = 2Yn end (65) — 
In solving equation (65), it is convenient to take y, in some unit of 1, disregarding decimal point, 


since the solution may then be written down by inspection as proceeding by units of 1 (y, = 0 by 
equation (52), that is 


i af 
12.0 0 
12.1 0.1 
12.2 0.2 


12.3 0.3 
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The differentiation operator developed in equation (55) in terms of central differences is not 
convenient for use in calculating values of tg from a numerical tabulation of the stress function 
y. The operator can be derived in terms of forward or backward difference operators. If the neces- 
sary analyticity of a function is assumed and it is expanded in a Taylors series, one immediately 
establishes the symbolic identity relating the differentiation operator D the shifting operator 
€ p = p (r +h) and the base of the natural logarithms: 


e=ehD, Di = Ine; D=+Ine. 


If A is the forward difference operator, then 
Ap(r)=pir+h)—pe), A=e—1, e=1+44 
and 
D=—In(1+4), 


expanding the logarithmic term 


1 S m 1 n 
Ey ght apo eo) 
; (66) 
1 1 
D=7,(4 eoree ie tee Aas ), 
then 
= (sy —4 ary +4 Ary (61) 
to0= 5 (4¥ — > sale deny I 


Near the end of a tabulation, the corresponding formula in terns of the backward difference operator 
P(r) = p(r) — p (r —h) may be used: 

ey 

i a ~ ont+l1 


iain 
(68) 


1 1 1 
too= 7 (DV + Vets yt). 


To determine the error introduced by approximating the general differential equation to obtain the 
stress field given by equations (29)—(31), applications of the formulas derived earlier will now be 
made. 


5. Application 1. Consider the hollow cylinder of Fig. 1 sae 
This corresponds to the first example given by Trostel, who assumes the following laws for variation 
of physical properties: 


A(T) =4,—A, T, (70) 
o(T) =% +o,T, (71) 
E(T) = E, — E,T? (72) 
and gives the following constants for mild steel 
Ay= 0.12 [cal/cm? — sec — (°C/cm)], 
A, = 0.00007 [cal/cm? — sec — (°C?/cm)], 
Se = 12 x 10- RC hE 
= 0.001 x 10-5 ey ad 4 k ag Geometry of the hollow cylinder. 
eee ee PaGe ee ee ees 
E 13.2 [kg/cm? — °C?)]. : 


a) Temperature field. The steady state temperature field is easily established, on the as- 
sumption of law (70), to be 


rok |i—n)-[e-mI--) 


1 1 1 la fe) 
\ ro 


(73) 


16* 
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For the case of constant thermal conductivity A, > 0 and equation (73) reduces to 


(74) 


As pointed out by Trostel, the temperature fields given by equation (73) and (74) differ only slightly 
for the example under consideration. This may be seen by inspection of the first two columns in 
Table 1. Trostel capitalized on this fact to derive closed expressions for the integrals ¢,(r), y(r), 
and @p,(r) (it should be noted that the final expressions listed by Trostel are not complete). Actually 
with the exception of —,(r) which is indefinite, the integrals are about as easily solved by some ~ 
form of numerical quadrature. 


b) Stress field. Using the temperature function (74) and the defining equation (19), p,(r) can 
be expressed as 


“(T)— T;) 
g,() = ee (1 a8 T,) — ae in oe ) (75) 
2In 2) In a : 
Yo ro 

The stress function was determined numerically by use of the simplified equations (64) and 
(65) for the case of y = 1/4 (Trostel uses y = 1/2 in his solution). Equation (14) was then used to 
determine the radial component 7,, of the stress vector. Equations (15), (67), and (68) were em- 
ployed to calculate the rotational component of the stress vector Tog. Numerical details are 
contained in Appendix 1. The results are entered in columns four and five of Table 1. In the 
solution, the space increment was taken as h = 0.10 and 60 steps used to complete the calculation. 
In using the approximate solutions, equations (29) and (30), the definite integrals ¢y,(r) and @,(r) 
were solved by numerical quadrature and the expression (75) was employed for y,(r). The results of 
this second calculation are entered in columns six and eight of Table 1. Deviations from the numeri- 
cal results are given as a % error. Inspection of this error shows that the approximations consistently 


give a stress value which is about 1% too small in magnitude. This however is within the experimen- 
tal error of the constants employed in equations (70)—(72). 


Results assuming constant physical properties and using equations (33) and (34) are given in 
the last four columns of Table 1, with their corresponding error. 


Table 1. First Application 
r, = 12cm, rp = 18cm, T; = 350° C, Ty = 280° C, »y = 1/4, p; = pp = 0 


Temperature, °C Stress Field in kg/cm? 


Radius r|constant | variable Numerically Approximations by eqns. Constant property | Constant property, 
Ora conduc- | conduc- by eqns. 64 & 65 29 & 30 eqn. 34 eqn. 33 
tivity | tivity % 
Enq. 74 | Eqn. 73 TOO TOO TOO Error Gap Error 
: 0 pe 
12.5 : : 2) — 46.7|\— 7.9 
13.0 : faa : ' All .4|/— 81.9|— 5.5 
13.5 -7| 329.3) —1103 | — 595.8} — 592.5] —0.6/—110.0} —0.3] — 549.8; —7.7|—100.0|--—10.0 
14,0 | 323.4) 323.0] —123.3 | — 372.2) — 371.5] —0.2|—122.8} —0.3] — 343.6] —7.7|—113.8|— 8.7 
14,5 | 317.3} 316.9] —127.6 | — 128.5) — 127.8|.—0.4|—127.0} —0.5] — 118.1| —8.7]—118.8|— 6.9 
Ike) |) ay eles —124.4 |} + 81.8) + 79.7} —2.6|—124.0} —0.4] + 75.9 8.9 | —115.0. ges 
15.5 | 305.8] 305.4] —114.8 | + 271.2 269.8} —0.5 | —114.2] —0.5 257.8 | —5.0]—104.9|— 8.6 
16.0 | 300.3 — 99.7 | + 490.0 480.0 0|/— 99.7 0 436.4) —7.6]— 92.5)— 7.3 
16.5 | 295.0| 294.7] — 80.2 | + 636.5 635.1} —0.2|— 80.0} —0.2 587.6 | —7.8]— 73.5|— 8.3 
17.0 | 289.9 — 90.6} -- 805:5 801.0) —0.6|— 55.9; —1.1 739.5 | —8.2]— 51.2|— 9.5 
17.5 | 284.4) 284.3] — 29.8 | + 962.8 961.3; —0.2|— 29.5 | —0.7 881.9 | —8.4|— 27.8)— 6.7 
18.0 | 280.0} 280.0 0.0 | +1087.1 1086.0} —0.1 0.0 0.0 1018.9 | —6.5 0.0 0.0 


6. Application 2. We now consider an example which will provide a more severe test of the 
approximate solutions derived in the first of my former paper [5]. Consider a composite body 
made up of two concentric, hollow cylinders, denoted as cylinder (1) and cylinder (2), as shown in 
Fig. 2. If the inner cylinder (1) is taken to be a carbon steel and the outer cylinder (2) a low alloy 
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steel, then the pertinent physical properties vary with temperature as shown in F igures 3, 4, and 5 
in which the bracketed numbers identify the cylinders. 


a) Temperature field. In general for the steady state case, we wish to solve the problem 
ACI 3 Me (77) 
where //? is the Laplacian operator subject to the/two-point boundary conditions 
Pett, A sl E, (78) 
r=r: T=T,. (79) 
Expanding equation (77) 
Va. VT +HAV*T =0 


E(T) x10° [kg/ené] 


700 200-300 400 600 600 700°C 
Fig. 2. Geometry of two concentric hollow cylinders. Fig. 3. Young’s Modulus 
yp) = v2) = 1/4 
1 — a = 12 em 7 8 pi =0 
{ 7D — b= 16cm Tee = ee 
r® _ 6 = 16cm TT — T* 
{ rf) — ¢ = 20cm 700% pi=.0 


which reduces in the case of axial and rotational symmetry to 
di dT d?T Lede 
a a( )=o 


dr dr dr?! or dr 


} (Tag) = 9. 


or dr dr 
dT 
Th = Cis) 
dre 1 
dr qr’ 
fUgaae Cae (81) 
Rises ey 


1 At this point one may wish to introduce an explicit function for the thermal conductivity and solve the 
problem as 


fadT=c r+ e, 


however, the non-linear nature of the thermal conductivity in this specific example tends to make the 
algebra messy; hence, the iteration procedure developed above is employed. 
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ee 
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from equation (78) 


from equation (79) 


so that 


Introducing (82) and (83) into (81) we 


dr 
ee 
eer ear 


have 


Ten Teer yee 
dr dr 
ee Ar 
vt; vy 
dr 
Ar 
T— T; os Len 
T, — T; ets 


The last equation is immediately adopted to the composite body at hand by setting: 


Cylinder 
Ta=6 
fo = 0 
f ppd ih 
tad 


The desired temperature fields become 
for cylinder 1 


for cylinder 2 


1 Cylinder 2 
Tri 6 
[6 
f SR e bee 
T= T, 


4 Keb a T; ve 
(i 
: dr 
Ar 
Tr 
dr 
KO 
TO _ 7T* |, 
RET PE 


dr 
1@; 


b 


The common temperature T* is determined from the conditions 


—AQ (b) ( 


dT e: ‘ dT® 
dr )_.- — 12 ()( dr Is 


(82) 


(83) 


(84) 


(85) 


(86) 
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from which 
b 


dr 
| 
des Lie ots): b ; ¢ ' 
dr dr 
fast | a 
a b 


A(T) [cal/sec-cmé-(°C/cm)] 


(87) 


o« (7) x 10° [°C] 
150 


(1 


~ 140 + 


130 


- Temperature 


0 00 +~« 200 +~=«C SCO SC« OSC «OC 00599 200 300-400 +~«800+~=«600~=«*70O™C 


Fig. 4. Coefficient of Thermal Expansion Fig. 5. Thermal Conductivity 


Table 2. Second Application 


Stress Field in kg/cm? 


Approximations by eqn. Constant properties, eqns. 
(third iteration 29 & 30 33 & 34 
of eqn. 85, 86) % %, %, of, 

Thr Error TOO Error Thr Error L7oYo) Error 
es 

12.0 0 0 23.399 0 = 22.964 | —1.9 0 aa 20.709 | —11.5 
12.5 33.0 886.0 20.990)|| 875.7) —I1 20.491 |, —1.9 | 777.4 | —12.2 18.444 | —12.2 
13.0 67.5 1611.2 18.535/1582.0 | —1.8 18.232 | —1.6 |1443.0 | —16.8 16.398 | —11.5 
13.5 102.4 2197.1 16.266 |2148.1 | —2.2 16.087 | —1V.1)1917.7 | 27.2 14.529 | —10.7 
14.0 137.0 2661.5 14.274 |2606.5 | —2.2 14,132 | —1.0 |2353.5 | —11.6 12.768 | —10.5 
14.5 173.2 3029.3 12.099 |2968.7 | —2.0 12.002 | —0.8 |2693.7 | —11.0 11.176) —9.0 
15.0 211.5 3299.0 10.247 |3192.9 | —2.3 10.153 | —1.0 |2947.9 | —11.0 9.688 | — 5.6 
15.5 254.6 13493.8 9.188 |3433.1 | —1.8 - 9,202 0.1 |3136.6 | —10.0 8.311) — 9.5 
16.0 295.2 3621.0 7.657 |3621.0 0.0 1,632 |2—-0.3:/3288.1 | — 9.0 7.021) — 8.3 
16.0 295.2 3621.0 | —11.345 |3621.0 0.0 | —11.322 | —0.4 |3288.1 | — 9.0 | — 7.709 | —32.0 
16.5 345.1 3218.3 | —12.244 |3170.2 | —1.3 | —12.229 | —0.2 |2967.2| — 8.1 | — 9.186) —25.0 
17.0 389.5 2720.1 | —13.210 |2689.0 | —1.1 | —13.196 | —0.6|2538.4| — 6.8 | —11.350) —14.1 
725 435.2 2238.6 | —14.176 |2217.9| —1.0| —14.123 | —1.0|2150.1| — 4.0 | —12.024) —15.2 
18.0 481.8 1769.0 | —15.187 |1740.9 | —1.6| —14.946 | —2.3 |1743.2| — 1.5 | —13.317) —12.3 
18.5 Dal.7 1312.3 | —15.875 |1281.1 | —2.4| —15.781) —1.0 )1325.5 1.0 | —14.551| — 8.3 
19.0 584.4 881.6 | —16.292| 861.0} —2.2} —16.163| —1.3) 872.9 0.8 | —15.756| — 8.8 
ys) 642.7 444.1} —16.975 | 440.7} —0.8| —16.782| —2.0|} 440.3); — 1.0 | —16.884|) — 0.4 
20.0 700.0 0.0 | —18.018 0.0 0.0 | —17.860 | —1.6 0.0 0.0 | —17.967| — 0.1 
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The second law of the mean for an integral may be applied to equation (87) to obtain the initial 
estimate of 7* 
() 
In | — 
a 


,@ 


7 = T,4+(T)—T)) (88) 


AD ae 2 


Equation (88) gave an initial estimate T* = 300° C, and three subsequent iterations of equa- 
tions (85) and (86), starting with a linear gradient, gave T* = 295.2 °C, and the temperature field 
shown in column 2 of Table 2. 


b) Stress Field. If cylinders 1 and 2 are disengaged, then it becomes necessary to impose on 
the external surface of cylinder 1 a pressure, p, to prevent it from expanding, and at the same time, 
a pressure of equal magnitude but opposite direction on the internal surface of cylinder 2 to prevent 
it from contracting. The stress field then follows from the approximations (29) and (30) for cylinder 
1 by setting p) = —p. For this particular example p; = 0. The field for cylinder 2 is obtained by 
setting p; = — p. Also p) = 0. The equations are then 
for cylinder 1 


LD Pil’) 1 pir) Par) | 
ce =p (Fe) Be 1 —»y® P3(b) Go a 6} 2 6) 
E E 
— + 9,(r) 1 2a =\- P(r) 
. oe r Polr) + P3(r) 
ify = ‘| ae +7 ya Pall) ca ferme a 
for cylinder 2 
GS) 22 __ Alt) 1 Pir) __ alr) 

Trr =p (1 ey a 1 — y® pale) ne aon ¢ 0) 
E E 
= at) ; (= + lr) 

(2) weal: r me UU) ) 

Lat d ( Plc) + 1 — »® Pale) Pi(c) 3(¢) ee 


Since the temperature field in this example was established numerically, it is necessary to inte- 
grate the expression ¢,(r), equation (19), in like manner. The second form of the indefinite integral 


0 = Tr) 
is more convenient for this purpose. The integrand r f{ «(@) dO in the indefinite integral can be 


quite accurately represented by some form of polynomial approximation. In this particular example, 
the first four Gram polynomials of degree zero, one, two, and three, respectively, were employed. 

The interfacial pressure p is arbitrary until the tolerance between the two cylinders is specified. 
Denoting this tolerance by Au: 


Au = u)(b) — ul?)(b) 


where the radial component of the displacement vector is given by equation (32). For this exemple 
v1) = y) = 1/4 so that the value of axial strain component, A,, isimmaterial. The expression for the 
tolerance is a function of the unknown pressure p as is seen on evaluating equations (91) and (92) at 
r=b 


© = T* 
Au tdo(b) — r0(b) 1 1 
oe eee rou) Seley = Gy) eed) | alot he 

@=0 

(93) 

For the case of constant physical properties we have 
Au 760(b)  1(b) 1 1 
= 099) (Se 20 \—9 a +) p( ds — pn) + +9) a) TH. (04) 


For this example, take Ju = 0.10 cm, then equation (93) gives p = 3621 kg/cm?, while equation 
(94) gives p = 3288 kg/cm?. The complete stress field as calculated numerically, by the approximate 
extensions of Trostel’s solutions, and by the equations for constant physical properties is given in 


ao «yee 


‘ 
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Table 2. The average error introduced by the approximations in this example is seen to be larger 
than in the first application, but still less than 3° which is within the experimental accuracy of data 


over this wide a range. Again the approximations give stress values which are slightly too small in 
magnitude. 


7. Conclusions. The extensions of Trostel’s solutions derived in this paper may be employed 
tor general y over wide conditions with a resulting error less than 3%, the stress values being too 


small in magnitude by this amount. The error decreases as the variation of physical properties of 
the media decreases or as y—> 1/2. 


(Eingegangen am 6. Oktober 1958.) 


Anschrift der Verfasser: Prof. Dr. Ing., Dr. Phil. Milomir M. Stanisié, West Lafayette, Indiana (USA), Purdue 
University, Division of Engineering Sciences; Graduate Student R. M. McKinley, West Lafayette, Indiana (USA) 
Purdue University, Chemical Engineering. 
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Zur Theorie der Rotationsschalen 
vom Standpunkt numerischer Rechnungen 
Von E. Klingbeil 


1. Einleitung. Die folgenden Ausfiihrungen sind veranlaBt durch Uberlegungen! dariiber, wie 
sich statische Berechnungen von biegesteifen Rotationsschalen mit drehsymmetrischer Belastung 
auf elektronische Rechenautomaten programmieren lassen. 

H. Reifsner? behandelte speziell die Kugelschale. E. Meifner stellte erst fir Rotationsschalen 
konstanter Wandstirke®, spater auch verdnderlicher Wandstarke*, zwei gekoppelte Differential- 
gleichungen zweiter Ordnung auf, die sich zu einer Differentialgleichung vierter Ordnung zusammen- 
fassen lassen. Fiir spezielle Schalenformen gelingt es ihm, die Differentialgleichung vierter Ordnung 
auf eine komplexe Differentialgleichung zweiter Ordnung zu reduzieren. F’. Tolke® fuhrte spater 
unter Vernachlassigung sehr kleiner Terme in den Meifnerschen Gleichungen das Problem auch 
fiir beliebige Schalenformen auf eine komplexe Differentialgleichung zweiter Ordnung zuriick. Fiir 
praktische Berechnungen interessieren numerische Verfahren. Ein solches hat W. Lohmann® 
angegeben. Er fiihrt die Meifnersche Differentialgleichung vierter Ordnung unter Vernachlassigung 
kleiner Terme auf eine komplexe Riccatische Differentialgleichung erster Ordnung zuriick. Nach- 
dem diese auf numerischem Wege gelést ist, werden die vernachlassigten Terme durch Iterationen 
wieder mit beriicksichtigt. Einen neuen Weg beschritt H. Miinz’ mit Variationsrechnung. Auf dem 
Wege iiber die Hamiltonsche Theorie der kanonischen Differentialgleichungen entstehen vier Dif- 
ferentialgleichungen erster Ordnung, die unmittelbar auf numerischem Wege gelést werden kénnen. 

Eine Rotationsschale ist geometrisch durch ihren Meridian 
bestimmt. Die Form des Meridians l]48t sich auf verschiedene 
Weise durch Koordinaten festlegen. W. Lohmann und H. Miinz 
verwendeten rechtwinklige K oordinaten x, y. Fiir eine steile Schale 
(siehe Abb. 1) fiihren sie die Ordinate y, die gleichzeitig Breiten- 
kreishalbmesser ist, als Funktion der unabhingigen Verander- 
lichen x ein. Fiir die flache Schale dagegen benutzt Miinz den 
Breitenkreishalbmesser y als unabhangige Verdnderliche, um 
Singularitaéten der Ableitungen zu vermeiden. So braucht man fir 
eine beliebig gekriimmte Rotationsschale mindestens zwei ver- 
schiedene Verfahren, d.h. zwei verschiedene Systeme von 
Differentialgleichungen. F'. Wirthwein® hat in einer Dissertation einen mehrfach gekriimmten 
Flaschenboden behandelt. Er benutzt vier getrennte Verfahren fiir Kugelschale, steile Schale, 
flache Schale, Zylinderschale. Dadurch werden seine Berechnungen sehr umfangreich. 

Zur Programmierung des Schalenproblems fiir elektronische Rechenautomaten ist es aber 
zweckmafig, die einzelnen Bereiche der Schale mit ein und demselben Verfahren zu behandeln. 
Man braucht dann nur ein einziges Programm und spart auf diese Weise an Speicherzellen. 

Diese Vereinfachung wird durch die Einfiihrung der Meridianbogenlange s als unabhangige Ver- 
anderliche gewahrleistet. Die Funktion y(s) ist iiberall eindeutig und differentiierbar. Die Unter- 
scheidung zwischen flacher Schale und steiler Schale entfallt damit. AuBer F. Télke, der bereits 
mit der Bogenlange gearbeitet hat, bedienen sich alle unten angefiihrten Autoren anderer 
Koordinaten. Ich lege zunachst dar, wie sich die Gleichungen der genannten Verfahren schreiben, 
wenn die Meridianbogenlange als unabhangige Veranderliche verwendet wird. Dann sage ich noch 
einiges tiber die Lésungsmethoden, iiberlasse jedoch dem Leser, fiir praktische Berechnungen oder 
Programmierungen unter den Verfahren auszuwahlen. 

1 Den Herren Prof. Dr. A. Walther (Darmstadt), Dozent Dr.-Ing. habil. W. Lohmann (Aachen) und Dr. 
H. Bottenbruch (Darmstadt) danke ich fiir Anregungen und Aussprachen. 

* H. ReiBner, Spannungen in Kugelschalen, Miiller-Breslau-Festschrift, S$. 181. Leipzig 1912. 

3 E. MeiBner, Phys. Zs. 14 (1913) S. 343. 

* E. MeiBner, Vjschr. Naturforsch. Ges. Ziirich, 6 (1915) S. 23. 

° F. Télke, Ing.-Archiv, 9 (1938) S. 282. 

6 W. Lohmann, Ing.-Arch. 6 (1935) S. 338, 

7 H. Miinz, Ing.-Arch., 19 (1951) S. 103 und 255. 


* F. Wirthwein, Elastizitatstheoretische Untersuchung eines mehrfach gekriimmten Flaschenbodens. 
Dissertation TH Aachen 1955, 


Abb. 1. Koordinaten 
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A. Mechanik der biegesteifen Rotationsschale im Koordinatensystem (Ss, y) 
2. Voraussetzungen. Wir setzen foleende Annahmen fest: 
a) Es gelten die Voraussetzungen der einfachen Balkenstatik. 


b) Es wird angenommen, daf die Schalennormale auch nach der Verformung noch senkrecht 
auf der Schale steht, d. h. da keine Querkraftverformungen auftreten. 


c) Die Wandstarke der Schale ist klein gegen die Kriimmungsradien. Daher reicht es aus, die 


geometrischen Beziehungen, die Verformungen und Verschiebungen am Element der Schalen- 
mittelflache zu betrachten. 


d) Es gilt das Hookesche Gesetz. 


3. Definitionen. Wir definieren folgende GréBen: 


a) Geometrische GréBen (Abb. 2) 
s = Meridianbogenlange, 
y = Halbmesser des Breitenkreises, 
r, = Hauptkriimmungshalbmesser der Schale in Meridianrichtung, 
r, = Hauptkriimmungshalbmesser der Schale in Breitenkreisrichtung, 
gy = Winkel zwischen der Schalennormalen und der Rotationsachse. 


Die HilfsgréB8en 9, r,, r, lassen sich durch folgende Beziehungen in y und s ausdriicken (s. Abb. 2). 
cosp=y', sing =/1—y"?, 
PS pea aon 
aaa Ta Vien 
wobei Striche Ableitungen nach s bedeuten. 


a 


b) SchnittgréBen: 


Ns = Normalkraft je Langeneinheit eines Schnittes ? = konst., (mit ? als Breitenkreiswinkel) ; 
positiv, wenn sie die Schale auf Zug beansprucht, 


N, = Normalkraft je Langeneinheit eines Schnittes s = konst., positiv, wenn sie die Schale auf 
Zug beansprucht, 

Q = Querkraft je Langeneinheit eines Schnittes s = konst. (Vorzeichen nach Abb. 2), 

My, = Biegemoment je Langeneinheit eines Schnittes } = konst., positiv, wenn es die Innen- 
flache der Schale auf Zug beansprucht, - 

M, = Biegemoment je Langeneinheit eines Schnittes s = konst., positiv, wenn es beim Umlauf 
mit s im Uhrzeigersinn Zugbeanspruchung in der Schaleninnenflache erzeugt. 


c) Verschiebungen und Verformungen: 
v = Verschiebung in s-Richtung (positiv in s-Richtung), 


w = Verschiebung senkrecht zur s-Richtung (umlauft man mit s die Schale im Uhrzeigersinn, so 
zeigt das positive w aus der Schale heraus), 

&, = Meridiandehnung, 

€, = Breitenkreisdehnung, 

y = Verdrehung der Meridiantangente (positiv bei Abnahme des Winkels 9), 

#, = Abnahme der Meridiankriimmung, 

% 3 = Abnahme der Breitenkreiskriimmung. 


Nsds 
Mg ds 


+2 


6(M, 
Mya asad 
yd d+ LVasae 


aly 
gya 2 Wacan 


Abb. 2. Bezeichnungen Abb. 3. Gleichgewicht am Schalenelement 
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4, Gleichgewichtsbedingungen. Aus Abb. 3 lassen sich leicht die Gleichgewichtsbedingungen 
ableiten. Man betrachtet das Gleichgewicht der Krafte in Richtung der Meridiantangente und in 
Richtung der Schalennormalen, auBerdem das Momentengleichgewicht um eine Breitenkreis- 
tangente. Nach Eliminieren der HilfsgréBen —, r,, rz bekommt man die drei Gleichgewichtsbedin- 
gungen in der Form 


(NS) NG oe OE poy 07 


(fae 
— /f I 
Oy) + NoVl—y'* —N, ae — Pay =O () 
OL) Maia Core 


5. Das Elastizititsgesetz. Nach dem Hookeschen Gesetz ergibt sich ohne Beriicksichtigung 
der Querkraftverformungen 


N,=D(e, +29) ee eye 
N,=D(é,+ v¢,); M, = K (#, + v %,) 


mit h = Schalenstarke und 


(II) 


Eh he 
estes K= 0 * 


~~ 22 2 


6. Die Forminderungsbedingungen. Die Herleitung dieser Gleichungen sei hier nicht angegeben, 
sondern hierfiir sowie fiir das bessere Verstandnis der beiden Ziffern 4 und 5 verweise ich auf das 
Lehrbuch von Fliigge!, in dem sich ahnliche Herleitungen befinden: 


gate, m= 
Ls 
vy’ +wyl—y? t 
oe oe. L JON mane, (111) 
Lie eae: 


B. Die Differentialgleichungen der biegesteifen Rotationsschale 
im Koordinatensystem (Ss, y) 


7. Die Meipnerschen Differentialgleichungen mit der Bogenlange als unabhangige Veranderliche. 
Entsprechend dem Wege von Meifner fiihrt man die Variablen 


= pearly 
See ete cas 


ein. Dann ergeben sich durch geschickte Umformungen aus den Gleichgewichtsbedingungen (I), 
dem Elastizitatsgesetz (II) und den Formanderungsbedingungen (III) die zwei Meifnerschen 
Differentialgleichungen mit der Bogenlange als unabhangige Veradnderliche. Auch hierbei verweise 
ich auf das Lehrbuch von Fliigge, in dem die Herleitung, allerdings in anderen Koordinaten, kurz 
beschrieben ist. Mit dem Linearoperator L 


1 mr hy”? 
Te) Oe ee 
Oe qa vitae 
und den Koeffizienten 
00) = ae mer (3h'y’ +hy"), 
1 ur , , iad , / 
O(s) =~ [2 (hy thy) —o (hy —h yh 


yl—y 


1 W, Fligge, Statik und Dynamik der Schalen, 2. Auflage. Berlin 1957. 
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_ergeben sich die Differentialgleichungen 


KU) +9u=~C—")y, 


LV)+ 2Q2V=—EU-+4@. 
Das Belastungsglied ® hat die Form 


=F wo) :: a) ots Chee oe } 


mit 


Beye Psy 


AO (ji—y2)3 faye 


Darin bedeutet F das Integral 
=f p.yy' ds— f p,yV1—y? ds. 


In Anlehnung an die Arbeit von Lohmann! vereinfacht sich das Belastungsglied © im Falle des 
konstanten Innendrucks p zu 


2) = +0 
mit 
p2Y Wl Ben 1 ii 
10) = san (Gr —yP” yet) tPA 
ees vy” px 1 ED 
00) = a9n (Gray yo) Pa 


Z(s) = px (y? — 6) - 


8. Die Differentialgleichung von Télke. Wie in Ziff. 1 erwahnt, hat F’. Télke in seiner Ab- 
handlung bereits in folgender Weise mit der Bogenlange gearbeitet: 
Fiihrt man in den Meifnerschen Gleichungen mit der Bogenlange als unabhingige Veranderliche 
die Transformation 


12 (1 — 2?) Lith 
y= ER V5" van VSé 
aus, und schreibt man das Belastungsglied um mit 
’ ap 
YP se Saag ras ’ 
so entstehen mit den Koeffizienten 
3 /y/\2 Ly” h’ 1 yw Lk’ 
als) = 3 (7 Basse we vr y peed es eres gs 
Ph 1 kh” 
«’(s) = vr + ie TA ie & 


eee ee 
unter Vernachlassigung von «’ die Ditiecoaualeleictanven von Tolke 
eal os 
t—at+pn=—t//2 | 
Sie lassen sich komplex zusammenfassen: 
x” + (iB —a) X =2//% : 


1 Siehe FuBnote 6 von Seite 242, 


\ : r 
246  E. Klingbeil: Theorie der Rotationsschalen vom Standpunkt numerischer Rechnungen _Ingenieur-Archiv _ 


Kennt man die Lésungen fiir X 
Xa) =X +i X,, Xe =X +i Xy, Xp) = Xp, + 1X,» 

so ergeben sich die allgemeinen Lésungen fiir ¢ und 7 zu 

¢ =— X,,— K, X, + K, X, — K, X3 + Ky X, 

n= + Xp, + K, Xy + Ky X, + Ky X3 + Ky X,. 

9. Die Miinzschen Schalengleichungen im Koordinatensystem (s, y). a) Das Variations- 

problem. Der folgende Weg findet sich in anderen Koordinaten ausfiihrlicher bei H. Miinz. Die 
potentielle Energie der belasteten Schale 


v=20{ B (N, €, + No &s + M,%, + My %s) — (p, v + p, w)| ¥ ds 
muB nach dem bekannten Satz der Mechanik ein Minimum sein: 
oll == (0) 


Diese Forderung steht in Analogie zu dem Hamiltonschen Prinzip der Dynamik: 
t, 
6 | Idt = 0% 
to 


Daher 1aBt sich die Hamiltonsche Theorie anwenden. 

Jetzt wird der Integrand L mit Hilfe des Elastizitatsgesetzes (II) und der Formanderungs- | 
bedingungen (III) durch v, w, y und deren Ableitungen v’, w’, y’ ausgedriickt. Da v, w, y nicht von- 
einander unabhangig sind, muf iiber ihre Vertraglichkeitsbedingung ein ,,Lagrange-Faktor“ A 


eingefiihrt werden: 


L=>{|D be ope aes ae? ee \( eo 
—y'2 Ly y 
vy te N=7V aly Za 2(Xy 
+( y Be Arm) Metre 3 5 


Jt 


tay (w' papemey eaten) 


Mit den Impulsen 


Po= Gr =D y(¥ — pete] + oy Fully) =r, 
aas 

Pw ow’ af =—yQ, 
oL , ; 

Prep ao Sed re) =yM, 


entsteht aus der Formel 
Hl (P.» Pw Px Vv, W, x) = Pv v aig Pw w’ is Px x! aay L (v’, w', oa v,W, Xx A) 
die Hamiltonsche Funktion 
v ye w 47 


es ee AP Riad ALT. 
Pw yl—y? ee jl—y? 


1 Py \2 Pian oe ena 2 ES ah 
+ hp [By 2» Boy + wy =y9)— 1 (oy + T= 


eat 


P,\2 Faken 
+K |) SO eayaanel —0—*) 2y'l +y Po +P.) 


Nach der Vorschrift 


a » _. OH 
Oy OPy td 

poe ee nOH 
ow’ tropa: 

Pi rie re aes 


; 


* 


XXVH. Band 1959 E, Klingbeil: Theorie der Rotationsschalen vom Standpunkt numerischer Rechnungen 247 


ergeben sich die kanonischen Differentialgleichungen 


rite a ' rears VP) 
ee ee oP ey ot Ly) —p,y, 


fin Bs Mae Bess i } qo) VI? 
Se a he) OY we E99) — py, 
2 Cope sae ae 
Pie Piigee i Pe ot BL oe) yo 
y 
; WN = x77 yl—y? Pv 
v ypv— = Ww —_ 
3 g mse ery )+D5> 
Nea 
DSO Senge 
4 eres 
7 , ae Px 
i Py Ky ~ 


b) Ein bekanntes Integral. Durch Eliminieren von N, aus den beiden ersten Gleichgewichts- 


bedingungen entsteht ein unmittelbar quadrierbares Integral, das weiter oben schon auftrat und 
mit F bezeichnet wurde: 


(Ny Vl—y? + Oxy’)! =pyy’ —peyVl—y? 


oder 
(p,/1—y? —p,¥')’ =p.xy’ —p,yVl—y® 
Mit der ,,erzeugenden Funktion“ 
V =(p,/1—y’* — pay’) a + (Py + Pwo V1 —Y') G2 + Py Ws 


entsteht eine kanonische ete mit deren Hilfe man die ee des Integrals ausnutzen 
kann: 


ov ov Te Serio: , 
Pate Pell) Pad Uae eg i) uae ey 
OV OV. 7 if 
Pa = Gq = Po’ + poVl—y?, we a tay 
OV 
"fe AS aa A= ap, 2° 
Mit den Abkiirzungen ‘ 
K, =p, /l—y? —By', K, =p,y +p./1—y? 
ergibt sich daraus eine neue Hamiltonsche Funktion 
PLyl— ys pay) 
H, = H—= = 4(p,j1—y? 2— py’) +7 {D i Heese 


soe Lei : ; , 
—2yq, et IE (v4) gi] + K |(2) 294 By 


— (1—»?) y” a|| + y (Ky a + Ke %)- 


Hieraus entstehen wieder sechs kanonische Differentialgleichungen. Da zwei sofort zu integrieren 
sind, fiihrt das Problem jetzt nur noch auf ein System von vier gekoppelten Differentialgleichungen. 
Hat man also die erste kanonische Gleichung 


Pi=—yK, 


integriert, so ist jetzt das System von vier Differentialgleichungen zu lésen: 


p2 —? Po py ae = te (ha?) D—yK,, 


¥ 
Pig? Ps ce Pall y's yA pyy + (l=—-*) K = 955 
; 0 fee dicen oF I lon 
+9 G7 = gVl—y? + a} D : y’ 
, Ve 1 
4s TO Ga sy PP 
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Die letzte kanonische Differentialgleichung la8t sich wieder einfach integrieren. Sie lautet 


; i 72 +p y’ l1—y’? l1—y’? 
pepe) - : l y Ge =, 
c) Dimensionslose GréBen. Mit einer mittleren Wandstarke h) und einem mittleren Breiten- 


kreishalbmesser y, der Schale fihrt man dimensionslose GréSen ein 
= = ele ee 
Gs Yow . So 4 hy ? Son 
Jetzt bedeuten Striche Differentiationen nach ¢. Wahlt man die neuen Veranderlichen P,, P2, Ps, 


Q,, Oz, 03 und verwendet die erzeugende Funktion 


Y i aan 
y120—>»*) y 1 
V= E hi y, Yo P11 + y Pa Oe oe 7 PQs| , 
so ergibt sich nach der Vorschrift 
iar le 1 _ wv 
A biome a Ehgy, 2X ~ apy 


eine weitere kanonische Transformation 
Eh? 


4 
Dee q, = 12 (L—+*) 49» 
yi2 ( —») 


Sn Pe =jR~— 12(1—*) Qa, 
y12 dQ — ?) 
Eh (en ee il 
3 = 7 0 Ps» qs = 12 (l1— »*) — Q,. 
12 (1 — 2?) " 
Hieraus erhalt man durch Transformation von H, die neue Hamiltonsche Funktion 
HI. OV 1 See a7 D 
Hy = ig — = OlPanvl—e* — 1@) 
1 1a » en Ta: ' 
wed) ae eo C ke Nite) 5 Osh al Ue eat 
yizd—v) 1 yi2Q—2) 1 0 i ee 
= 2A ae 2A cat Saal eee AX 


12 0 —2) : 
es [Kia + Ko7Qs) + 7 (PaQe + P; Qs) « 


Man errechnet daraus auf bekannte Weise wieder sechs kanonische Differentialgleichungen. Zwei 
davon sind direkt zu integrieren. Eine sofort, die andere nach Auflésen des verbleibenden Systems 
von vier Differentialgleichungen. Die direkt zu integrierenden Gleichungen lauten 


1 
Pi =——Fp— eK, 


, 1— , yl—o? = aaa In} f yl— oe”? 
== A P21 64 22? ho ee 0 ee 
Qi ae ( iV Ca, 2 0’) meee U ne Qz - 


Es verbleibt das System von vier Differentialgleichungen erster Ordnung fiir die Variablen P,, 

Ps, Qe, Qs: 

Se 
A 


a llth @ ela sie i j 
ad ard bed red Mg mya oer 


4 / esa Oo ee tS 
@—(“_»£}a, = oy > (\PaVl—e? + Pro’ n) + @Vl—e®, 


Ot? (Ase ba Yost 2X2) 
= on 1 E? 


Q 
mee 


7 Q j 
MES ER ats y2d—r) 1 

es GRA EA 7 eee 
Q; e A) 7 oy Pa: 


Das sind die Miinzschen Schalengleichungen in dimensionsloser Form, hier aber mit der »»bezogenen 
Bogenlange“ ¢ = s/y) als unabhangige Veranderliche. 
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C. Numerische Lésungsmethoden 


10. Die Integration der Miinzschen Differentialgleichungen. a) Lésungen des homogenen 
Problems. Zunachst betrachtet man das homogene System von vier Differentialgleichungen 
erster Ordnung mit K, = K, = P, = 0 als Anfangswertproblem. Man integriert dieses System 
von vier Differentialgleichungen viermal mit einem der bekannten numerischen Verfahren. Und 
zwar gibt man dafiir die vier Quadrupel von Anfangsbedingungen vor: 


| Q, P, Qs Ps 


|S pee Meta ( Sle & ae () 


ee, Ob 8 Oy 
Die so gefundenen Lésungen des homogenen Problems sollen heiBen: 
eam 
a) alien 
2 2 3 
b) Eine Lésung des inhomogenen Problems. Wie H. Miinz} auf Seite 255 ff. schreibt, 
errechnen sich aus den obigen Lésungen des homogenen Problems die Partikularlésungen Q8, P9, 
Q3, P3 des inhomogenen Problems. Fiir die Programmierung ist jedoch folgendes Verfahren geeig- 
neter. 


Man lést einfach das inhomogene System von vier Differentialgleichungen mit der Anfangs- 


bedingung 5 
1 Oe Te 
2 2 3 3 
0 0530 


c) Das Randwertproblem. Da es sich um ein lineares Randwertproblem handeit, hat eine 
allgemeine Lésung die Form 
Poa PP Ge PhaeC, PPC, PP 4c, Pe. 
Durch die Randbedingungen lassen sich die Konstanten C, bestimmen. 
Aus Griinden der Genauigkeit ist zu empfehlen, eine langere Schale in mehrere Bereiche zu 
unterteilen. Aus Rand- und Ubergangshedingungen folgt dann ein lineares Gleichungssystem 
fiir die Integrationskonstanten der einzelnen Bereiche. 


11. Die Integration der Meifsnerschen Differentialgleichungen. Das soeben angegebene Ver- 
fahren 14Bt sich analog auch zur Integration der Meifnerschen Differentialgleichungen heranziehen, 
da es sich hierbei ebenfalls um ein lineares Randwertproblem handelt. 

AuBerdem gibt es noch das bereits in Ziff. 1 erwahnte Verfahren von W. Lohmann. Dabei 
wird zunachst das homogene Problem unter Vernachlassigung von @ und {2 auf eine Riccatische 
Differentialgleichung zuriickgefiihrt. Nachdem diese numerisch integriert ist, errechnen sich die 
Lésungen des homogenen Problems, die unter Beriicksichtigung von © und Q verbessert werden. 
Eine Partikularlésung des inhomogenen Problems findet man schlieBlich nach der Methode der 
»,Variation der Konstanten‘“‘. Die Programmierung fiir elektronische Rechenautomaten erfordert 
einen héheren Arbeitsaufwand als man zunadchst voraussieht. Deshalb ist fiir Rechenautomaten 
mehr zu empfehlen, das System von Differentialgleichungen zweiter Ordnung, wie soeben beschrie- 
ben, direkt auf numerischem Wege anzugehen. 


12. Der Bereich der Achsennahe. Im Punkte y = 0 versagen die hier angegebenen numerischen 
Verfahren, da y bei den maBgebenden Gleichungen im Nenner auftritt. Aus diesem Grunde empfiehlt 
es sich, im Bereich der Achsenniahe die Schale durch eine Platte (‘Tangentialebene) oder durch eine 
Kugelschale (Kugelkalotte) zu ersetzen. F. Wirthwein fiihrt bei seinen Berechnungen des Flaschen- 
bodens eine Kugelschale in Achsennahe ein. Er vereinfacht die dabei auftretende Differential- 
gleichung, indem er den Bereich der Kugelschale als sehr klein annimmt. 


(Eingegangen am 16. Oktober 1958.) 


Anschrift des Verfassers: Eberhard Klingbeil, Frankfurt a. M., Bergerweg am Heiligenstock. 
Aus dem Institut fiir Praktische Mathematik der Technischen Hochschule in Darmstadt, Prof. Dr. A. Walther. 


1 Siehe FuBnote 7 von Seite 242. 
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Anwendung der Deltamatrizen auf inhomogene Probleme 


Von E. Pestel 


Herrn Folke K. G. Odqvist zur Vollendung seines 60. Lebensjahres 


Das Verfahren der Delta- oder Determinantenmatrizen, das von Fuhrke' fiir die Untersuchung 
von Eigenschwingungen von Staben und Rahmentragwerken entwickelt wurde, laBt sich in ein- 
facher Weise auf die Behandlung erzwungener Schwingungen und statischer Probleme erweitern. 
In verschiedenen Arbeiten? wurde gezeigt, daB sich bei der Untersuchung der letztgenannten Auf- 
gaben die gewéhnlichen Ubertragungsmatrizen um eine Spalte und eine Zeile vergréBern. Im Falle 
des ebenen Biegeproblems ist dann die sogenannte erweiterte Ubertragungsmatrix eine 5 X 5-Matrix. 
Um alle Unterdeterminanten zweiter Ordnung dieser erweiterten Ubertragungsmatrix in das von 
Fuhrke angegebene Schema? der Deltamatrix einzuordnen, ist das Lexikon fiir die Anordnung der | 
Unterdeterminanten wie folgt zu erweitern: 


1 | 2 
ime 


3 
1,4 


4 
2,3 


sect. a a 


Die Deltamatrix (1) der erweiterten Ubertragungsmatrix 


| G, | & | 
GC, | & 


(1) 


wird damit eine 10 x 10-Matrix, die sich in vier Untermatrizen einteilen laBt. Die Tatsache, daB 
die zwei Untermatrizen ©, und ©, in den letzten vier Zeilen der erweiterten Deltamatrix eine Null- 
matrix bzw. der homogene Teil 1 der erweiterten Ubertragungsmatrix I sind, folgt daraus, da} 


die letzte Zeile der erweiterten Ubertragungsmatrix (2) 


Sv 


| 


Uy1 Ujo-Ujg Uig Uy 
Ue, Use Ue3 Unoa Us 
Uz, Ugg Ug3 Uz4 U3 
U4, Ugg = Ugg U44 U4 

0 0 0 0 il 


vier Nullen und eine Eins enthilt. 
Die Untermatrix (©, die die ersten sechs Zeilen und die ersten sechs Spalten umfaBt, ist identisch — 
mit der Deltamatrix der freien Biegeschwingung, wahrend die Untermatrix ©, den EinfluB der 


Belastung wiedergibt. 
Wenn man — im Gegensatz zu Fuhrke — die Probleme ,,von links nach rechts“ 
dann entsprechen die Spalten 1, 2, 5 und 6 der erweiterten Deltamatrix (3) 


L.S. —- |R-s a 4 (—») w. M Q 
ia UieeOrs tUya: Ui. Useio ae a ee ee 
Uny Ugg Usg Usy Usg Ung | Unz Usg Use Uno 
Us; Us -Uzs. Use U5. Use «| > Usy, 2Ug, 0 Usg Us, 
Uy, Ug, Uys Use Uss Une Usz Ugg Ugg Uno 
Us: Us, Usg Usa Usy Use Us, Usg Us, Uso 
Us: Use Uss Use Uss Use Us, Uss Uso Uso 

0 0 0 0 0 0 

0, "1055.0 Ogee Olea G u 
0 0 0 0 0 0 
0 0 0 0 0 0 


Eds Re Ing.-Arch, 23 (1955) S. 329, 
® vgl.z. B.S. Falk, Ing.-Arch. 24 (1956) S. 
336 


3 H. Fuhrke, a. a. O. S.3 


(2) 


durchrechnet, 


(3) 


216 und E. Pestel und G. Schumpich, Schiffstechnik 4 (1957) S. 55. 
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_bestimmten am linken Ende des Stabes vorhandenen Randbedingungen, wahrend die Zeilen 1, 2, 5 

und 6 jeweils bestimmten Randbedingungen am rechten Balkenende entsprechen. Wie aus der 
Bildung der Elemente der Deltamatrizen gemaf dem obigen Lexikon ersichtlich ist, sind die Spalten 

_ 7, 8, 9 und 0 der erweiterten Deltamatrix (3) den Anfangswerten am linken Balkenende zugeordnet. 
Diese Anfangswerte, von denen infolge der Randbedingungen am linken Balkenende jeweils nur 
zwei unbekannt sind, lassen sich dann als der Quotient zweier Zahlen ermitteln. 


Abb. 1. 


Die Anwendung der erweiterten Deltamatrix auf Probleme erzwungener Biegeschwingungen 
von Staben oder entsprechender statischer Probleme 14Bt sich am einfachsten an einem Beispiel 
_ zeigen. Zu diesem Zwecke behandeln wir die in Abb. 1 dargestellte Aufgabe. 
Zunichst stellen wir die drei gewohnlichen erweiterten Ubertragungsmatrizen zusammen: 


a) Ubertragungsmatrix fiir das Feld J zwischen den Punkten J und 2: 


east a ee | ean Ol 
S 2 Ed, 6EI, 24E I, 
hy i | q: hi 
“ t je 7a Oa a ey 1 a 
q: ii 
F>P* 6 1 Ee —43 
0 0 0 1 —q:], 
0 0 0 0 1 


b) Punktmatrix fiir die Federstiitze in Punkt 2: 


1 0 0 0 0 | 
0 i 0 0 0 
0 0 1 0 M (5) 
—C 0 0 1 0 
Or 0s 0 A ia 


c) Ubertragungsmatrix fiir das Feld 2 zwischen den Punkten 2 und 3: 


jee Z : <7 | 

2E I, 6E I, OE I, 
0 I - ET ay Aa 

Ean, 

2 2 | 2 (6) 

0 0 1 i 3 /P 5G! 
0 0 0 LS — Pp 
0 0 0 0 1 N. 


Bevor wir die entsprechenden Deltamatrizen aufstellen, ist es niitzlich, die Auf lagerbedingungen 
am linken und rechten Balkenende zu beriicksichtigen. Da die linke Balkenseite eingespannt ist, 
_ brauchen wir entsprechend (3) von dem homogenen Teil der erweiterten Deltamatrix fiir das Feld 1 


hele 
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nur die sechste Spalte. Im Hinblick darauf, daB wegen der linksseitigen Auflagerbedingung Ww, 
und y, Null sind, d.h. also nur Moment und Querkraft von Null verschieden sind, benétigen wir 
auBerdem nur noch die neunte und zehnte Spalte der erweiterten Deltamatrix fiir Feld 1. Ent- 
sprechend der bekannten Regel fiir die Bestimmung der Elemente der Deltamatrizen ergeben sich 
unter Zuhilfenahme des obigen Lexikons die sechste, neunte und zehnte Spalte der erweiterten 
Deltamatrix fiir Feld 1, wie folgt: 


Pgh big Oe hand Uae eae ger oes 
12 (E 1)? 24 (E I,)? 144 (E I)? 
i 5 qi ie Bes 
aw ons 24 EI, 24 EI, 
R q la 3qit 
pag 8 Bi eT, 24 EI, 
Ii oe ii ey Cie 
SEI) a EM IS 6 
Ly gli Cove ed i ; | 
TEE mest, wae 8 (7) 
qi | 
: i ke Ea 
0 Ii i 
ra a 6E I, 
0 ki R 
Ba 2h 
0 1 L, 
0 0 ] 


Die erweiterte Deltamatrix fiir die Federstiitze, an deren Ort auch ein duBeres Biegemoment 
angreift, wird mit Hilfe des eingangs angegebenen Lexikons auf Grund von (5) bestimmt: | 
a 


1-20 e000 Cui Meret 
0. Haz 0: 2O.A 00 Sed OenOR EO 
Gas Tt 0s 0 ne 0 G20 0,050 
0280-2 0 le 02880 0 M 0 0 
Co st0rat 0 se cO Sean 3p eeO | 
08, Ce 07 0 Seen CM0 0-M (8)| 
Oo Orn Oe 0m 20 Troe Oe s0 
0: 30 g028. 0s. 0 erO Tp teu dite 21h 
te viens 0-0 sO 

| 0¢°30°.0% (02 0.80 ain a OmeeO meg 


Da infolge der Auflagerbedingungen am rechten Balkenende y, und Q, Null sind, braucht von) 
der erweiterten Deltamatrix fiir das Feld 2 nur die fiinfte Zeile berechnet zu werden, die folgendes: 
Aussehen hat: 


a9 


(9) 


P(R—@ 


Fiir die weitere Rechnung, die zwar auch im allgemeinen Zahlen durchgefiihrt werden kénnte,| 
wollen wir fiir die Daten des in Abb. 1 dargestellten Tragwerkes folgende Zahlenwerte verwenden:| 
q=l1, M=4, C=1, P=3, =2, (1235 ¢d=2,7ks)— oe EJ, = 1/24 
Indem wir diese Werte in die erweiterten Deltamatrizen einsetzen und diese entsprechend dem| 
iiblichen Multiplikationsschema von Matrizen anordnen, erhalten wir folgendes Rechenschema: 


4 ad 
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(6) (9) (10) 2 


Sat ean cee 
Aeeeee 1-5 Bae 
eR oe = 3 26 
Fees o 53 
Z Sy a eee ay 
a era —3 
Gnas 5 
GL 9s.8 6 
Ob igwees 3 
ome sae | 1 

Mees 0 E-0,) 0520 Pe TT [Pe eat es an Fs (20) 
Seats 260 0.0 0 NT ae aa cea ha lee / 
Bae oi ai Ys G7 074-0 0 202) S20 Pale Ma / 
0) APR stay eens beers eee Ghee M0 0 (eee / 
Mee 050-1 <0 On tO etree 0 Aiveatys bea 5 
Peerios Or 00. «04. eR hie ar ee) so ae 10 
je 0 oO Oe 6 hae SiS NIE feepot eee / 
pee 602 0, 02, 0 (Vics bag eae 0 38 3 6 
aor 0-0 Oe 0 Merb haa haan hee 3 
Seon aD 0 0 2a 0d 0 1 hy 
(Su Sea es (ee ee eee Ee ee ee ee 


Mit Hilfe der untersten drei Zahlenwerte, die wir nach Durchfiihrung der Matrizenmultiplikation 
erhalten, lassen sich dann M, und Q, wie folgt bestimmen: 


35 39 
Beso bette 36: 


M, = 
Als Rechenregel ist hierbei beriicksichtigt, daB Zahlenwerte aus den Spalten 7 und 9 mit negativem 
Vorzeichen eingefiihrt werden miissen, eine Regel, die der bei der Lésung inhomogener linearer 
Gleichungen mit Hilfe von Determinanten zu beachtenden Regel entspricht. 

Nachdem nunmehr der Zustandsvektor am linken Balkenende bekannt ist, kénnen wir die 
Zustandsvektoren an der elastischen Zwischenstiitze und am rechten Balkenende mit Hilfe der 
gewohnlichen Ubertragungsmatrizen (4), (5) und (6) wie folgt errechnen: 

0 |j=—w, 


0 at A 
eae 


= ae == Q, 
1 aa 


Oe Oiag0 Ae 0.) a) 1 |= (11) 
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i alts 0 tN ee R 
1 0 0 0 0. SoeTt => W3 
294 
hears! 0 0s ac) —35 |\=vF 
41 | aap 
108 “ 
—l 0 0 1 0 aecry =Q3 
0 ON SO SPO tae he henmerl el e 
(= al 786 
Nes RP eae t= <3) 0 |= 
67 
O10 “brine as eo +4, |=M, 
O A0p.05 0 FAL Bee OM ane 


Durch die Verwendung von Deltamatrizen wird somit auch bei inhomogenen Problemen der Vor- 
teil erzielt, daB die am Anfang der Rechnung eingefiihrten zwei Unbekannten nicht als Quotient 
zweier Determinanten wie im Falle der Verwendung gewohnlicher Ubertragungsmatrizen berechnet 
werden miissen, sondern sich als Quotienten zweier Zahlen ergeben. Damit wird auch die Gefahr 
ernster numerischer Fehler vermieden, die sich aus der Tatsache ergeben kénnen, da haufig die 
Nennerdeterminante bei der Anwendung des Verfahrens der gewéhnlichen Ubertragungsmatrizen 
eine kleine Differenz groBer Zahlen wird. 


(Eingegangen am 22. Oktober 1958.) 
Anschrift des Verfassers: Prof. Dr.-Ing. E. Pestel, Hannover, Am Welfengarten 1. 
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Zur Theorie des raumlich gekriimmten Stabes* 
Von G. Kimmel 


1. Einleitung. Der raumlich gekriimmte Stab war bereits oft Gegenstand theoretischer Unter- 
suchungen. Diese Arbeiten gehen jedoch zumeist von Voraussetzungen aus, die dem zu unter- 
suchenden, mehr oder weniger speziellen Problem geniigen. Die Kriimmungshalbmesser werden in 
der Regel als sehr gros gegeniiber den Querschnittskoordinaten des Stabes angenommen, so daB die 


Spannungen iiber die Biegesteifigkeit des geraden Tragers gleichen Querschnitts ermittelt werden 


kénnen. In der vorliegenden Arbeit wird diese Voraussetzung nicht gemacht. Es wird gezeigt, daB es, 


-aufbauend auf den Uberlegungen von F. K. G. Odqvist' fiir den Kreisringtrager, d. h. unter exakter 


Formulierung der Bernoullischen Hypothese, auch fiir einen raumlich gekriimmten Stab méglich 


ist, die Dehnung und somit den Normalspannungszustand jeder einzelnen Stabfaser in Abhangigkeit 


von der Belastung anzugeben. 


2. Die Geometrie der Raumkurve. Die Stabachse des Tragers sei in einem rechtsgewundenen 
kartesischen Koordinatensystem (x,, xj, x3) durch die Parameterdarstellung 


x, = %,(t) (k = 1, 2,3) (1) 
-gegeben. Das Bogenelement ds der Raumkurve erscheint dann in der Form? 
Pe Cee 
a3 =|. a, (2) 


Bezeichnet man mit x die erste Kriimmung und mit t die zweite Kriimmung oder Windung der 
Raumkurve, wobei 


d?x);, d2xp, ies 
M4 == aye = VX Xe, (3.1) 
und 
lA tt ttt 
__ bifk Xy Xj Xe (3.2) 
Serre 
BP 


ist, so folgen die Vektoren des normierten begleitenden Dreibeins der Raumkurve zu 
, Ib ae 
T= 475 Hy T;, By = €:j% T; Hi, - (4) 


Dabei bezeichnen T; den Tangentenvektor, H; den Hauptnormalenvektor und B; den Binormalen- 


_ vektor der Raumkurve. Mit (4) lauten die Frenetschen Formeln 


1—72H,, H,;=—xT,+7B;, B, =—vH,;. (5) 


3. Zusammenhinge zwischen der verformten und unyerformten Raumkurve. Ein Punkt P der 
Stabachse unterliege dem Verschiebungsvektor r,. Bezeichnet man mit u,v und w die Betrage der 
Komponenten des Verschiebungsvektors in Richtung der Tangente, Hauptnormalen und Bi- 
normalen, so kann r, (Abb. 1) in der Form 

r, =u T, +0 H, + w B, (6) 


geschrieben werden. Der Tangentenvektor T,, der verformten Stabachse folgt dann mit (4), (5) 
und (6) zu , 
T,=T,(1+w —x0)+H, («+0 —rw) +B, (c0+w'). (7) 


Man sieht sofort, daB oa keinen Einheitsvektor darstellt, denn seine Norm ist, wenn die Ver- 
schiebungen als klein vorausgesetzt werden, 


T, T, =1+2w—2x0, (8.1) 


* Gekiirzte Fassung einer von der Bergakademie Freiberg genehmigten Dissertation. Referenten: Prof. 
Dr.-Ing. habil. D. Riidiger, Prof. Dr.-Ing. A. Kneschke. 

1 F. K.G. Odqvist, Ing.-Arch. 22 (1954) S. 98. 

2 Es wird hier und im Folgenden vorausgesetzt, da in Ausdriicken, die zwei gleiche Indizes enthalten, tiber 
diese Indizes entsprechend den drei Dimensionen des Raumes von 1 bis 3 zu summieren ist. 
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und sein Betrag ergibt sich zu 


\ 


VT, T,=1+w —x0. (8.2) | 
Der Vektor f;, der Verdrehung des Stabelementes bezogen auf die unverformte Stabachse (Abb. 2) 
14Bt sich analog zu (6) in der Form - 
B.=O0T, +p A + B, (9) 


é 


Abb. 1. Verschiebungsvektor des Punktes P. 


Abb. 2. Verdrehungsvektor des Stabelementes. 


darstellen, wenn unter #, y und @ die Betrage der Komponenten von f,, in Richtung der Tangente, 


Hauptnormalen und Binormalen verstanden werden. Um eine Beziehung zwischen den Drehwinkeln ~ 
g und py und den Verschiebungen u, v und w herzuleiten, wird der Winkel y zwischen dem Tangenten- — 


vektor 7;, der unverformten Stabachse und dem Tangentenvektor T;, der verformten Stabachse | 
ermittelt. Mit | 
sn | 
gilt dann 
__ eijk Tj T; 
\TpTp 


i 


Da andererseits jedoch auch 


¥,=yp H, + ¢ B; (10.2) 


ree 
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- ist, so folet 


ope 
T, Tp 
Unter Verwendung von (7) und (8.2) folgt mit den bekannten Beziehungen 
B, = 64,7, H,, T; = &4; Hy, B,, HH, = ¢,,, B, T, 


- unmittelbar 
y H;+ 9 B =—(t0 +4 w’) H, + (wu + vo! — rw) B,. 
_ Damit lauten die gesuchten Beziehungen 


y=—(tv+w), g=(xu+v —tw). (11) 


___ 4. Die Verschiebungen jedes einzelnen Querschnittspunktes. Bisher wurde der Trager nur durch 
_ seine Stabachse beschrieben. Nun soll der Verschiebungszustand jedes einzelnen Querschnitts- 

punktes dargestellt werden unter der Annahme, daf sich der Querschnitt langs s nicht andert. 
_Jeder Punkt Q des Tragers wird (Abb. 3) in der Form 


q 
x, = %, + 9, = % + H, + 7 B, (12) 
dargestellt. 


Abb. 3. Verschiebungsvektor des Punktes Q. 


Unter der Voraussetzung der Bernoullischen Hypothese vom Ebenbleiben der Querschnitte, 
kann man den Tragerquerschnitt als starres Gebilde auffassen, dessen Bewegung durch eine Trans- 
lation r;, und eine Rotation f;, festgelegt ist. Mit Hilfe des Drehtensors' D,, gilt fiir den Vektor q, 
(Abb. 3) folglich die Beziehung 

% = Did - (13.1) 
Da die Drehung f), als VerformungsgréfBe klein ist, erscheint mit 
snp » Bp, cos fp & 1 

der Drehtensor in der Form 

Dy = O41 — Exim Pm « (13.2) 
Damit lautet Gleichung (13.1) 

Ik = Uk — Ekim Pm N° (13.3) 
Das skalare Produkt T,, q, verschwindet, das heibt T,, und q;, sind orthogonale Vektoren. Jeder 
Querschnitt, der vor der Verformung senkrecht auf dem Tangentenvektor steht, steht demnach 


auch nach der Verformung senkrecht auf dem Tangentenvektor 7, der verformten Stabachse. 
Der Verschiebungsvektor |, folgt nach Abb. 3 zu 


h, =, + 4k — Me - (14.1) 


1 Vgl.z. B. A. Duschek u. A. Hochrainer, Grundziige der Tensorrechnung in analytischer Darstellung. Bd. 1. 
5. 78, Wien 1948. 
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Unter Beriicksichtigung von (13.3), (14.1), (6), (9) und (12) ergibt sich 
= —Ep +7) T, + & 79) H, + ww +€ 9) B,. (14.2) 


Wird der Verschiebungsvektor |, durch seine Komponenten u,, v, und w, in Richtung der Tangente, 
Hauptnormalen und Binormalen der unverformten Stabachse dargestellt, so folgen die Betrage 
dieser Komponenten zu 


u(é,n) =u—Eptny, v(é,7) =o —79, w(é,n) =w+es. (15) 


Fiir € = 7 = 0 sind die Verschiebungen nach (15) mit den Verschiebungen u, v und w der Stabachse 
identisch. Beachtet man in (12), daB x,, H;,, und B, Funktionen der Bogenlange s sind, so ist 


x,(3) = «,(8) + € H,(8) +7 B,(8) (16.1) 


die Parameterdarstellung der Raumkurve einer Stabfaser § = konst., 7 = konst.. 
Das Langenelement ds dieser Faser folgt daraus zu 


ds =/(l—“é? + 72 (@ + 7%) d8. (16.2) 
Nach der Verformung geht ds iiber in 


= uw —Enptiny —xv+xn8 
Meee Es ERENT 
+ (1 2 ea 
ae Os) , 1/2 
tre(ipet | ds . (17) 
Die Dehnung ¢ dieser Faser betragt bekanntlich 

_ ds —ds 
Menta pea 


Mit (16.2) und (17) erscheint die Dehnung der Faser § = konst., 7 = konst. des raumlich ge- 
kriimmten Stabes in der Form 


1 Oo; ° , oOo; ° od Oy Ul Ua 
= G—x5it Ped ® “xu +é&(—o —xw +0 +Povitw)+&(xg +7) 


typ +xb—txu—rtV’ +? w) +2 t (# —xy)+én (—xy' —#20417x)]. (18.1) 


E 


Die Verdrehungen ¢ und y lassen sich daraus eliminieren, und man erhalt 


1 
eee (a) 


+E (OPW tune tuo! —xtw —xr wit) 4+n(—270 —7o— ww" 4x8 
—txb+ Pw) +7r(i+uto+uw)+éen(Qxty +uvo+tuw’ 2d 
+t?u—txw)]. (18.2) 


€ [uw —xo + (—2ew —xu—v"’ 42th +7 wt ed 47%d) 


In (18) ist mit t = 0 und x = konst. die von F. K. G. Odqvist! auf Grund der Bernoullischen Hypo- 
these aufgestellte Dehnung ¢, unter Beriicksichtigung der von A. Bjérklund? vorgenommenen 
Korrekturen, enthalten. 


5. Die Gleichgewichtsbedingungen. An einem Stabelement der Lange ds greifen an den Schnitt- 
ufern je drei Schnittkrafte und Schnittmomente an, die den wirkenden Spannungen Aquivalent 
sind. Schwerpunkt und Schubmittelpunkt des Querschnitts fallen zusammen. 

Die Belastungskomponenten in Richtung der Tangente, Hauptnormalen und Binormalen sind 
Pr» Py und pg. Bezeichnet man mit my, mg und mg stetig angreifende Momentenbelastungen um 


1 F. K. G. Odqvist, Ing.-Arch. 22 (1954) S. 98. 
2 A. Bjorklund, Ing.-Arch. 23 (1955) S. 421. 
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die Achsen des begleitenden Dreibeins, so kénnen die sechs Gleichgewichtsbedingungen nach H. G. 
Bechert' in der folgenden Form angeschrieben werden: 


N’—x Qu + pr =0, (19.1) 

Qu +%N—tQ3+ py =0, (19.2) 

. Qe +tQn + pp=0, (19.3) 
. My; —x My -+m,=0, (20.1) 
My + * Mp—t Mp—Q,3 + mg =0, (20.2) 

Mz; +tMy+Q_ +m, =0. (20.3) 


Myr diy 


Abb. 4. SchnittgréBen am Stabelement. 


6. Aufstellung der Differentialgleichungen des Verschiebungszustandes. Die Gleichung (16.1) 


q 
stellt die Raumkurve jeder beliebigen Faser des Stabes dar. Der Tangentenvektor T', an eine solche 
Raumkurve ist gegeben durch 


q 
T, = (1—x&) T,—ty H, +7é B,. (21.1) 
Er ist wie T,, in (7) kein Einheitsvektor, denn seine Norm folgt zu 
q 4 
1,1, =(—x2 +2 PR +8). (21.2) 


q 
Da die Dehnung ¢ in die Richtung von T,, fallt, kann die Dehnung jeder Stabfaser vektoriell in der 
Form 


Tk 
——— (22) 
V Tp Tp 
geschrieben werden. Bei Giltigkeit des Hookeschen Gesetzes ordnet die Dehnung ¢, nach (22) 
jedem Querschnittspunkt den Spannungsvektor 
Oo, = E E;, (23) 


zu. Die Zerlegung des Spannungsvektors nach den Richtungen des begleitenden Dreibeins liefert 
die Normalspannung o7 und die Schubspannungen oy und og 


ef — & 


Or —— E co es fe 6) = = 5 
Y= ea Ty 4 £8) 
ISS iT ee o4 
OF Vie ears ( ) 
tE 
OB 


=e : 
VG — x €)? + 0 (4? + &) 
1 H. G. Bechert. Zur Statik raumlich gekriimmter Trager. Diss. Karlsruhe 1954. 
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Die Darstellung der Spannungen durch die ihnen aquivalenten Schnittkrafte erfolgt auf be- 
kannte Weise zu 


F F F 


P 
Das Torsionsmoment M, setzt sich aus zwei Anteilen zusammen. Der primare Anteil M7, resultiert 
s 


aus den iiblichen de St. Venantschen Torsionsschubspannungen. Den sekundadren Anteil My er- 
zeugen die in (24) berechneten Schubspannungen oy und og. Es ergibt sich somit 


P s 
My = My + Mp, 


12 
F 


Mit (18), (24) und (25) sind die Langskraft N und die Schnittmomente My und Mg nach er- 
folgter Integration iiber die Querschnittsflache F’ als Funktion der Verschiebungen U, v, w der Stab- 


achse und der Verdrehung 9 um die Stabachse bekannt. Auch das sekundare Torsionsmoment My 


Pp 
ist in dieser Form gegeben. Der Zusammenhang zwischen dem Torsionsmoment M7 und dem Ver- 
schiebungszustand ist zundchst noch unbekannt. Dieser Zusammenhang ist jedoch zur eindeutigen 
Berechnung der zehn Unbekannten Funktionen 


u, v, Ww, 0, IN, Ou, Qz; M,, My, Mz 


notwendig, denn es stehen mit den sechs Gleichgewichtsbedingungen (19) und (20) und den drei 
Gleichungen (25) nur neun Gleichungen zur Verfiigung. Die fehlende zehnte Gleichung wird auf 


Grund der Bernoullischen Hypothese aus der Anderung des Drillwinkels um die Stabachse bei 


Fortschreiten um ds gewonnen. Ein Querschnitt verdreht sich um §,. Die Anderung des Dreh- 
winkels bei Fortschreiten um ds ist 


dp, = Bi ds. (26.1) 

Mit (9) und bei Benutzung der Frenetschen Formeln (5) ergibt sich die Winkelanderung zu 
dp, = [0 —xy) Ty + (y +2919) Hy + try) Bld. (26.2) 
Die tangentiale Komponente von df, stellt somit die Anderung des Drillwinkels dar. Fiir diese gilt 
(# —xy)ds. (27) 


Diese Winkelanderung am Stabelement ds mu mit dem primaren Torsionsmoment in Beziehung 


stehen, da die sekunddren Schubspannungen keinen Einflu8 auf den Drillwinkel haben. Sie ist 


formal mit der Winkelanderung am Kreisringtrager identisch. Fiir den Zusammenhang mit dem 
Torsionsmoment wird 


12 
; M 
oe —ny = ae (28) 


festgelegt, wobei in erster Annadherung die Torsionssteifigkeit GK von der Theorie gerade Stabe 
gleichen Querschnitts tibernommen wird. Damit folgt 


Mp = CK (O +r txt) — found + pope dk. (29) 


Gleichung (29) stellt die noch fehlende zehnte Gleichung dar. 


Damit sind der Spannungs- und Verschiebungszustand des raumlich gekriimmten Tragers ein- 
deutig gegeben. Fiihrt man fiir die in (25) und (29) nach Einsetzen von (24) und (18) auftretenden 
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Integrale tiber die Querschnittsflache F folgende Abkiirzungen ein 


i: dF os é 
J, al x €)2 + 72 (n? E2)]3/2? J, {d x éyP + 12 (72 £2)]3/2 dk 
Fr ' F 
ie = dF lpr A oe Sa ROL J 
8 [a — xP + (PL EB 4 | (l= &? + 2 (+ &)93/2 
F 
ie = «a dF = dF 
i Cea pieced Cpe ip | ais (las) woe 
F 
Py a COG, 2 Ave &4 d 
: = eee (4? + &))3/2 ee a [dx éP+ 7? + A)? ae 
F F 
- : (30) 
=| coat e@tmm (awe + eee 
F 
oa En En 
a CS ees! Cet) ee ao pis [CL — #6)? 08? + £") 8/2 ae 
F 
a Th Bn 
Ja =| oayeearteye sh [7+ ete 
F 


dF , 


Ju= | [(l — x €)? : 2 re si maple 


so erscheinen die Schnittkrafte und Schnittmomente als Funktionen der vier F ormanderungs- 
groBen der Stabachse in der Form 


N=Ei(—x# Jy +2un Ig— Ig —ut Jy $227 Jy — OT Sy) 
+E v (—% +? Jp +2 Jgp—v3 Jg—uv? J, tut? Iy~—P Ig —0' Jog $20 4 Jy —220' Spe) 
+ Ew(t' J,—2e0 Jz +20 Ig +02 Jy —20? x Sy +7? 2x? Sie) 
+ EO (Sy —2% Jy +2 Jy) + Ew (Jy —3% Jp + 32? Jz — 23 Js) 
ee J 0 J) 7 2 J, 3 Se ee Jo 
THOT, + I,—2Js—x Je) + Ev” (— J, + 2% J, —% Js) 
+ EW” (— Jn + 2% Jy — 2 Sep) > (31.1) 
My, = Et (—xt J, +221 JIg—@tIy—Hv Jy $2uv! Syg—ve x! Jus) 
+ Ev(—v Jyt2uv Ig 0' Jy Sg +t? Sy + 22? Sy — 23 Sg — 4? Sg FeO Syg— 2? 8? Sy) 
+ Ew (etIg,—240P Ig ted, +7 Jy 2H Sy +220’ Sy) 
+E Ox (Jg—2% Je +22 Jz) + EW (Jy —3 x Jy +32? Sy, — 28 Jy) 
tEv 2¢(—AtuJd,— J7) + Ew t (2 Jy — 3% Sy +4 Ig + 27 Sy — 2? Sis) 
+E OT (Sy + Jy3 —% Jy —% Sys) + Ev" (— Jy + 2% Syn, — 2? Sy) 
Bw (— J, + 24 Je—2 J) » (31.2) 


M, = Eu(e Jg—2xx Jp +22 x! Ig tut Jy — 2277 Sy + 2° t Jy) 
+ Eo (%J,-—V? Jg—2 27 Jn +28 Ig +4 Ig —4 0 Ig teed, +r Sy 2nd Sag +270’ Jy) 
+Ew(—v Jn +2uv JI;—20 Jg— Jy $+ 220? Sg — 2? 0? Sha) 
+ EO x (—Jy + 2% Jy —2 Jy) + Ew (—J, +3 Jg—3 %* Js + x? Je) 
+ BY 24 (Jy —2 4 Jy +2? Jy) + Ew' t (—2.5, + 3% 5,—% Je + 22 Jz — x? Js) 
+E Wt (—JIs—Je +2 Jz +% J) + Ev" (Jp —2 x Js + 27 Je) 
+ Ew" Jy —2% Jp, + Ju), (31.3) 
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My, = Et (—x' Js —u Je tun Jotun’ Jpg —2t Jy —ut Jy +t Sy +t Shs) 
$ Eben AS bOI PT TOI tO Ty tet Sy +t Sy—0 Jug Se 
/ ij; GK 
+ xT Ju +#T Jas + 4 | 
+ Ewt(t Jett Je—uv Jg—uet Ig tl Sy +B Sig —2 0? Jy — 2 0? Sis) 
+ E@xt (Sy + Jig —% Ju —% Sys) + Ev! 20 (— Ja — Js + % Su + % Ses) 
+ Ew t(J3 + J,—2% Js —2u% Ig +2? J, +27 Jo) 
LEW (QPS IAI nt Sy tet dy tat de tet ds +g | 


L Ee’ t(—JIg—Je +2 J, + Js) 


FERRO +A de L 72 J, ba | 


+ Ew" t (— Jz — Jy) - (31.4) 


Zusammen mit den Gleichgewichtsbedingungen (19) und (20) stellen die Gleichungen (31) die 
Differentialgleichungen des Spannungs- und Verschiebungszustandes eines Tragers dar, dessen 
Stabachse eine beliebige Raumkurve beschreibt. Die Gleichgewichtsbedingungen sind fiir jeden 
Belastungsfall integrierbar. Beim Randwertproblem zweiter Art (statisch bestimmter Fall) kann 
der Schnittkraftverlauf eindeutig angegeben werden. Die Gleichungen (31) liefern den dazugehéri- 
gen Verschiebungszustand. 

Beim Randwertproblem erster oder dritter Art (statisch unbestimmter Fall) ist der Schnitt- 
kraftverlauf bis auf gewisse Integrationskonstanten bestimmt. Mit dem entsprechenden Verschie- 
bungszustand nach (31) sind samtliche vorgegebenen Randbedingungen zu erfiillen. 


7. Anwendung der allgemeinen Theorie auf einen Trager, dessen Stabachse eine gewéhnliche 
Schraubenlinie beschreibt. Mit 9 als Drehwinkel und a als Halbmesser im Grundkreis lautet die 
Parameterdarstellung der Schraubenlinie 


x == @ c0s-0, Xo = ESI Oe Xs == C 0 (32) 


Das Bogenelement folgt mit 


Cc 
zu 
ds =ay/1+ ad. (34) 
Die Kriimmung ergibt sich in der Form ; 
if - 
und die Windung folgt zu 
@ 


Kriimmung und Windung sind konstante Werte langs der Bogenlange s. Die Gleichgewichts- 
bedingungen (19) lassen sich mit der Abkiirzung 


1 


2 
i 1+ w? 


(37) 
durch Elimination von N und Q, auf eine Differentialgleichung zweiter Ordnung mit konstanten 
Koeffizienten fiir Q, zurickfiihren. Sie lautet 


d? ted 
aot + On = 9 (pr— > GE + Pa), (38) 


und ihre Lésung ist 


Q y Q 
0, =e | (pp — @ pz) sin (90 —A) ta— | Pa cos (9 —A) dA +a (N + wQz) sin 0 Oe cose, (39.1) 
b 0 
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Mit der Lésung fiir die Querkraft Qy sind auch die beiden Schnittkrafte Q; und N bestimmt. Sie 
werden 


2 2 
Pe 2020 | (pr — © pz) cos (9 —A) dA + a0 J py sin (9 —2) da— aa J (w pr + pa) 
3 


+ w0?(—N + @Qz;) cos 9 —wa Og sing +0270, + oN, (39.2) 


ef (re) ps) 008 (0 dha | py'tin (9) A owe | (opp + pp) da 
0 f 
+ a2 (N — w Qz) cos 0 + 4 Oy sing + wa? (Qg + aN). (39.3) 


Dabei sind die Integrationskonstanten so gewahlt, daB die Schnittkrafte in Abhangigkeit der an 


der Stelle @ = 0 wirkenden Krafte N, ae und On erscheinen. Die Gleichgewichtsbedingungen (20) 
lassen sich durch Elimination von My und Mg, auf eine Differentialgleichung zweiter Ordnung 
mit konstanten Koeffizienten fiir M,, zuriickfiihren. Sie Jautet 


Ma a 1 dmg 
und ihre Lésung ist 
Q Q 


My =a (— a pz + a? mp — w a? mg) sin (o —A) a—z{ My COS (e —A) dh 
0 


0 


a 
—2@ ieee sin (A — 1) du sin (9 —A) dd 
00 


oa 
+S ° | [ Pe cos 2 — 1) du sin (9 —2) di + a0 (sin g + wg c0s 0) Op 
0 0 


—awosing Oy + 4a (sing —e cos g) N +. sing (w My — My) + cos 9 My. (41.1) 


Mit der Lésung fiir das Biegemoment My sind auch das Torsionsmoment M7 und das Biege- 
moment My berechnet: 


ee [2 a0 ie ha (Vw) pp = mop ow mal cos (0 —A) dA 
6 
ef @ a@ © Py—myz) sin (9 —A) ieee | [—2 aw pp—a (l—o?) pg —w? mp — mg] dd 
0 d 


A oA 
+2 aa wtf (pp — © pp) sin (A—u) du cos (9 —A) dA— 2 aos {Pa cos (A— 1) du cos (@ —A) da 
00 


+ ao? [wo sin 9 + (1 — w2) (1 — cos 0)] Qg + awa (—sing + @ cos Q) Oy 


+ ao a? [2 (1 — cos 0) — @ sin 9] N + wa? (l1—cos@) Mz + a? (cos 0 + w?) M,-+Mysine, 
(41.2) 


2 
My == | [2 (1 —0%) pr —2 aa? o py + oo mp — co a? my] eos (9 —2) dd 
fr 


4 af [a(1—o) py + wmg]sin(@—A) dd + a0 f [—a(1 —o) pp + 240 pp—Comg—ma] dA 
0 0 


A 
—2ara cof f (pp — @ pp) sin (A— p) dus cos (g —A) da 
06 


° 


A 
452 a c® { f py.coa (A —p) da cos (0 —A) dA + aw o? [— w @ sin g — 2 (1 — cos @)] Qp 
00 


+ aa? [o? 9 sin 9 + (1 —~ w2) (1 — cos g)] N + aa (— w? 0 cos 9 — sin Q) Oy 
+ wa? (1 —cos 9) Mp +02 (1 + @* cos g) Mp —wasing M,. (41.3) 
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Die Integrationskonstanten sind wiederum so eingefiihrt, daB die Momente in Abhangigkeit der 
an der Stelle 9 = 0 wirkenden Momente M7, My und Mg auftreten. 
Damit sind simtliche Schnittkrafte als Funktion des Belastungszustandes und der Randkrafte 
bzw. Randmomente bekannt. Fir den statisch bestimmt gestiitzten Trager sind die Randwerte 
N, Qy, Qp, My, My und Mz (42) 
aus den Gleichgewichtsbedingungen am Gesamttragwerk bekannt. lm statisch unbestimmten 
Falle kénnen diese Randwerte nur im Zusammenhang mit dem Verschiebungszustand ermittelt 
werden. Der Zusammenhang zwischen den Schnittkraften und den Verschiebungen ist durch (31) 


hergestellt. 

Der weiteren Berechnung wird ein symmetrischer Tragerquerschnitt zugrunde gelegt. Seine 
Symmetrieachse falle mit der Richtung der Hauptnormalen zusammen. Dann verschwinden die 
Integrale Jy)... Jy; (30). Die Gleichungen (31) erscheinen dann mit den dimensionslosen Ab- 
kirzungen 


A, =i [x J, + (7? + 2x?) Jy — x (0? + t) Js et? Ig — 2" 0" Jel 5 


a 


Ay = (U1 — 3% Jy + 3%? Jy —H* Js) » 
A, = 31 (2 Jp—3% Jy +2 Sg + Jg— 2? Je)» 
Ay =3t (Js + Ip—% Ig —% Jp) » 

A, ="; (— Jy $22 Sg — JI); 


Ag = aut (—Iy + 2x J,—# Jy), 


Ay = gle yp — (2 +228) Ig te OO +7) Is—n PI pee Si, 

A, = (— Jy + 3% Jg—3 7 Jy + 8 Je), (43) 
Ay ="F. (2 Jp + 3% Ig —H% Ie + Jy — 1? Jp) » 

Ay = + (—Js— Ig + % Sq + Se) » 

An = 3: (Ig—2% Sg +2 Ss) 

A, = = —x#J,—“J,+ (+7?) J, 4+? + 7°) Je putty put Sy tue ‘ 


Ay = =F (Jp + Jg— 2% Ig — 2% Je + 2? Sy + 0 Ss) » 


Ay = BPI +P Lath bath tee), 

Pa a2Pa +P te hte), 

in der Form 
N=ak(43 ! 4, L 1 ! Aa Asi) (44.1) 
My=@ EA,[i— ob — 39 ! is ! pore (44.2) 
ioe es (A, | Aaa = F Aeo Gp + An Ges) (44.3) 
My =o E(Ayyi + Ase + Au Ge + Au 9 Ge + Aue Gee): (44.4) | 


i ee 
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Nach Elimination von u und @ ergeben sich mit den dimensionslosen Abkiirzungen 
A,twoaA 
BA, Fo A, SL EIine gas 
1 2B 
B, = — ——, (43 + B, A), ee aa 
Ai, + wa? A 
Bs aaa By = Ay, + Bs A, 
1 B 
B, ee (Aj3 + Bs Ap) , Bg = wo? a. ) (45) 
B, = » Ay + Ay + B; (w A, + Ay) , 
1 
By == — 2 Ay + a? Ay + B; (—2 A, +0 A;)], 
By = 0 A, + A, + B, (@ A, + As) , 
1 

Byy = =3[— 20 Ag + a7 Ay =F B, (—2 A, + a As)] ? 

zwei Differentialgleichungen fiir > und w. Sie lauten 
° dw dy dw dMy B B, , M 
Pe Ba Pat 1 aa Bia a a do ee ate a (46) 
° dw dv aw dMy B viel 
Ei By + EG By + EG By + EG B= ey) ae 


Die simultanen Differentialgleichungen (46) lassen sich durch Elimination von w auf eine 
Differentialgleichung vierter Ordnung mit konstanten Koeffizienten fiir v zuriickfiihren. Man erhalt 
unter Beriicksichtigung der Beziehungen, die zwischen den einzelnen Koeffizienten B nach (45) 
bestehen, die Differentialgleichung 


d‘y dv ° dMy K CNG Kear? ae = @Mr K, sh 
pag aa ee = do a do? a 0? Pasi 24 M es do? at+M. 
(a) 
Die in (47) eingefitihrten Abkiirzungen K haben die nachstehende Bedeutung: 
6207.08 K _ B,B,— B,B 
Sa A, . 2° Bs By — B; Byz’ 
=, 3s Bo — BB, ne eee LB, 
3G (B,By—BBa)' “4B, By By By’ re 
K, ey A By n= Bs oe Oe. 
5 Bs; By— B, Bis : B, By— B, By,’ 
es By 
S is B,; By— B, B,, 
Die Lésung der Differentialgleichung (46) ergibt sich zu 
° 1 @ : 
Sag we | Ki {Ma (o —A) sin (0 —A) dd 
2 
+ [Na (K, — K,) + Mz (K,— K;) + My (K, — K7)] (@ — A) cos (9 —A) da 
0 
g . 
+ fla (Ky + Ky) + My (Ky + Ke) + Mp (Ky + K,)] sin (0 —A) di 
0 
+ % cosa +a(—ty + ww + ag) +49) sing + w aa (w Po — Vp) @ Cos E 
—aoyosine. (49.1) 
Fubrt man noch die Abkiirzungen 
A ee 2 as aA was 
ey ot ad,’ Ap watd,? (50) 


18 
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ein, so werden die anderen VerformungsgréBen erhalten. 


1) == — on KajeMre (o —A) cos (9 —A) da 
0 
a ea EN (K, — Ky) + Mg (K,— Ks) + My (Ke — K,)] (@ —A) sin (@ —A) a 
a $ ; (l— w*)a & 
ee A) da = JING Boome Bg toa [1 —cos (9 —A)] dA 
i @ 
Beach Bor Sesto ete [1 — cos (9 —A) da]} —o a8 sing 
+ wo? [—ty + w wy + 4G +4 Do] cos 0 + aa? (w py — Hp) (— w* @ sin @ + Cos Q) 
— aa yy (w? @ cos 0 + sin g) + 0? [w ty + wy) —2aw Gy +4 (1 —o*) A], (49.2) 
l 2 D 
OP eer ys — Ky J Ma ein (0 —A) Ao as) + M, (K,— K;) + My (K,—K;)] 
x Gosi(o A) di-ea (Ki 541) f LN a kya Ry ede he] ae 
+a [1 a? (I — 08) K, — Ky] J [Na Ky + Mg Ky + My Kj] da 
K. Q 
Tote Be theo AO ae Ug) 208 Oa Ny Oak a + Wo) (49.3) 
° 1 
b= oR —a Ky f My (9 —2) c0s (9 2) da 


as af [N a (K, — K;) + Mg (K,— Ks) + My (Kg — K,)] (@ —A) sin (9 —A) dh 


2 0 
+a K, J Maren oA) dA 2: foLN 6 Be MB sei Maal ete Cea 
Ee ia Dar 


a w? 


4 [IN aK, +My Ks + My K,] dA 
1 
+ xo? Ll — w?) (l— a? K,) — Ka | Wax, + M, K,+ M, Kaa 


@ 
K. ° . ° e 
+ Me, | Nad +05)sine—o8[—ty + oth +09 +00 0] cos @ 
0 


+ a @ a (0 9 — Bo) (Q sin @ + cos Q) — a W a YH (— Q cos. + sin @) 
+ @ a? [w Uy + Wy — 240M + 4 (1 — ow?) DB] +4 (wy +H) » (49.4) 


@ 
i 
yp =saaE|— Ki | Marcos (o —A) da 
s 0 
a i [N a (K, — K,) + My (K,— Ks) + My (K, — K,)] sin (g —A) da 
+ & (@ Pp — Vo) sin @ + yy cos, | (49.5) 
e @ 
C3 : 
? =soniK [ Me sin (g@ —A) dA + | [N a (K,— K;) + Mg (Ky —K;) + Mp (K,— K,)] 


(Ky 


om leah + SY [pat pK K, + My K,] dA 


1 
[1 —a? (1—w*) K,—K,] | [N a K, + My K,+M,K,]da 


tality | Nal tos (om—Mome—onmsing +et(od, + Po) « (49.6) | 
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Damit ist auch der statisch unbestimmte Fall gelést. Die Gleichungen (42) und (49) enthalten 
die den Integrationskonstanten entsprechenden zwélf Randwerte 


° ° ° ° °o ° ° ° ° 
IN, Qn Qp M,, My, Mz, Up, Vo, Wo Das Wor Po- 


Zu ihrer Berechnung stehen im allgemeinen Falle (es sind keine Randkrafte oder Randmomente 
bekannt) die sechs Gleichungen (49) zur Verfiigung. Der Trager ist sechsfach statisch unbestimmt. 
Der Stiitzungsart entsprechend kénnen jedoch schon verschiedene Schnittkrafte am Rande 9 = 0 
bekannt sein. Dann erniedrigt sich der Grad der statischen Unbestimmtheit um die Anzahl der 
bekannten Randwerte. 

Sind die unbekannten Schnittkrafte bzw. Verschicbungen bekannt, so folgt die Spannung 
o = 0(s, &,7) an jeder beliebigen Stelle des Stabes nach (24) mit (18). 


8. Zusammenfassung. Fiir einen Trager, dessen Stabachse nach einer beliebigen Raumkurve 
verlauft, werden die Beziehungen zwischen den Verformungen und den Dehnungen abgeleitet. Die 
Kriimmung und Windung jeder einzelnen Faser des Stabes sind dabei beriicksichtigt. Dem Ver- 
formungszustand wird konsequent die Bernoullische Hypothese ebenbleibender Querschnitte zu- 
grunde gelegt und der Querkrafteinflu®B vernachlassigt. Das Problem der Berechnung eines beliebig 
gekriimmten Stabes ist somit auf die Integration eines Systems von sechs linearen Differential- 
gleichungen zuriickgefiihrt. Fiir die gewéhnliche Schraubenlinie gelingt es, das Differential- 
gleichungssystem durch eine Differentialgleichung vierter Ordnung mit konstanten Koeffizienten 
zu ersetzen. Der Spannungs- und Verschiebungszustand kann als Funktion der Belastung explizit 
angegeben werden. Der iiber dem Querschnitt aufgetragene Spannungshiigel stellt keine Ebene, 
sondern eine gekriimmte Flache dar. 


(Eingegangen am 5. November 1958.) 
Anschrift des Verfassers: Dr. Giinter Kdmmel, Freiberg (Sachsen), Beuststr. 7. 
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Der Wellenwiderstand von Schiffen auf flachem Wasser 
Von F. Kolberg 


1. Einleitung. Seit Lambs! Untersuchung des Kinflusses der Zahigkeit auf freie Wasserwellen 
sind von verschiedenen Forschern (Havelock, Wigley, Inui) Arbeiten erschienen, welche den Zahig- 
keitseinfluB auf den Wellenwiderstand von Schiffen zu beriicksichtigen versuchen. Die dort ge- 
machten Vorschlige zur Erfassung der Zahigkeitseinfliisse sind halbempirischer Natur. Erst 
Srétensky® gelang zunachst fiir unendlich tiefes Wasser bei Linearisierung der zum Problem ge- 
hérigen Differentialgleichungen und Randbedingungen unter der Annahme, daB die freie Ober- 
flache keine Tangentialspannungen aufweist, die Angabe von Formeln zur Berechnung des Wellen- 
widerstandes und der Wasseroberflache unter Beriicksichtigung der Zahigkeit bei der stationéren 
Bewegung von Schiffskérpern an der Wasseroberfliche. Im folgenden beschaftigen wir uns mit 
der Erweiterung der Formeln von Srétensky auf flaches Wasser, wobei wir unter Annahme der 
Srétenskyschen Grundkonzeption auf die Herleitung der Randbedingungen an der freien Oberflache 
sowie die Aufstellung der Lésung des betreffenden Differentialgleichungssystems ausfihrlich 
eingehen. 

Bei unserer Herleitung der Formeln fiir Flachwasser sind zusatzlich zu Srétensky die Rand- 
bedingungen am Flachwasserboden zu beriicksichtigen. Wir verlangen, da am Flachwasserboden 
die Komponenten der Stérungsgeschwindigkeiten verschwinden. 

Einerseits untersuchen wir die Wellenbewegung, die durch Normaldriicke auf die freie Ober- 
flache (Hognersches Drucksystem) erzeugt wird, und berechnen den Widerstand dieses Normal- 
drucksystems sowohl im ebenen wie im dreidimensionalen Fall; zum anderen ersetzen wir das 
Schiff durch eine Quellsenkenbelegung der Mittschiffsebene und berechnen analog zu Michell den 
Druckwiderstand. Zu beachten ist hierbei, daB die gefundene Formel fiir den Wellenwiderstand 
die Bedingungen fiir zihe Fliissigkeiten, die sich auf die Schiffsoberflache selbst beziehen, nicht 
beriicksichtigt. 


2. Differentialgleichung des Problems und Lésungsmethode. Wir beziehen uns auf ein fest mit 
dem Schiff verbundenes kartesisches Koordinatensystem (x, y, z), dessen x-Achse in Bewegungs- 
richtung und dessen z-Achse senkrecht nach oben weise. Ist ¢ die konstante Geschwindigkeit des 
Schiffes, und sind u, v, w die Stérungsgeschwindigkeiten in x-, y-, z-Richtung, und ist p’ der Druck, 
so wird unser Problem beschrieben durch die Navier-Stokesschen Differentialgleichungen und die 
Kontinuitatsbedingung : 


Ou Ou Ou 1 dp’ 
| eed 
Cae a) ate tas I ee 0 Ox vy Au, 
Ov Ov Ov 1 dp’ 
(¢+tujs + v5 Maer cpe aera sora ell 
(1) 
Ow Ow Ow 1 op’ 
| 
(e+ u) = 4 Or ots es sg an ee 
Ou, dv Ow 
Ox a: dy of Oz Ue ; (2) 
Linearisierung dieses Systems ergibt 
ou Y = imop! 
Diets peas ace 
Ov y ie 1 op’ : 
Ow sey | ee ee 
Ox ocdz  c?’ 
Ou , ov. dw 
ox ayer a (2’) 


1 H, Lamb, Hydrodynamics, 6. Aufl. 1932. 


* L.-N. Srétensky, Sur la résistance due aux vagues d’un fluide visqueux, P di i 
behaviour of ships in a seaway, Wageningen (1957), puiemlacaiars se coe 
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Die Lésung unseres Problems setzen wir nun in der Form 


U=U+U, V=%y+%, w=w,+ 


an, WO Uy, V,, W,, p dem Gleichungssystem 


W., Pp’  =—p—ogz 


wane, Avi 0, 
wane, Aw == 0 
+ St + =o, Ap = 0 


geniigen sollen. Fiihrt man ein Potential ® ein und setzt 


aD 
Pee’. eva ’ 
wo also A® = 0 ist, so erhalt man nach (3) fiir die Geschwindigkeiten 
6 _ . 0@ 
Uy —— Sr oa 1 ay? 


Die Geschwindigkeitskomponenten uz, v,, w, miissen dann dem Differentialgleichungssystem 


Quy y 


pet ed VA 

__” Ay, = 0, 

oe Aw, = 0, 
ae ta + ae 


genugen. 


(3) 


(4) 


(5) 


3. Die Randbedingungen an der freien Oberflache und am Flachwasserboden. Beachtet man die 


Formeln 


; Ow Ou dw 
| oe +2Qvs, Pex = 09(55 + 3) 


Ov | 


Pry =0%(5 + 5 


(7) 


so lauten die Randbedingungen an der freien Oberflache bei unseren Voraussetzungen tiber das 
Verschwinden der Tangentialspannungen an der freien Oberflache, wenn T(x, y) das auf die freie 
Oberflache wirkende Normaldrucksystem und z = ¢(x, y) die Gleichung der freien Oberflache ist, 


Dew. Ou Ow 


D [pez — T(x, y)] =f 
Dt x 


fur Zz == (x, ¥). 
der Gleichungen (7) und (1’) schlieBlich 


Diese oe ae 


Ov Ow 


ae 


ou Y 2% Ow g Ih tl 2 = 
Ox ¢ duct Cc Oxdz r oe Qc Ox eee FY 
oder 
2v Ow g a lid : a * 
eget aca ee 


(8) 


Die erste dieser Gleichungen ergibt bei Linearisierung unter Beriicksichtigung 


(9) 


(10) 


(11) 
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Aus (10) und (11) folgt durch Integration nach x fiir die Wasseroberflache: 


ts) =(-0+N-Lz| (12) 


wobei beachtet wurde, daB fiir x > co die Werte von p, T und dw/éz gegen Null streben. 
Am Flachwasserboden z = —/h sind die Randbedingungen 
Us) OD = y w=0 (13) 


zu erfillen. 


Zusammenfassend haben wir also als Randbedingungen 


ou Ow Ov Ow 
Be oaeol? top doy age 


ae firs = 0, (14) 
2v Ow g a 
“¢ Ox Oz a ras Seprlte Aah 
wo fiir ein’ Michellsches Schiff T = 0 zu setzen ist, und 
i=. 0; v= 0; w=—0 tire -=—h, (15) 


4, Darstellung der Lésung, Berechnung des Widerstandes und Aufstellung der Gleichung fiir die 
freie Wasseroberflache. 


a) Michellsches Schiff. Ersetzt man das Schiff durch eine Quellsenkenbelegung der Dichte 
o(x, z) auf der Mittschiffsebene, so wird man in bekannter Weise! fiir das Potential ® auf den 
Ansatz 


P(x, y, 2) = a eerrose lia Gnas og at fae [ 
(x, 9,2) ra] [oC (5+) : tag | [aed pat f | 
exp {i K [(x —) cos @ + y sin O]} {F, Goj [K (2 + h)] + F, Sin [K (2 + A) dK 


gefiihrt, wo 
i= («—2? +y? + @—0}?, rh = (x —é2? + y? + (2 +2h4+0) 


ist. Fir z—¢€ > 0 und z+2h-+€> 0 haben wir hierfiir die Integraldarstellung 
1 - - 5 
@ =95 | [oe dg dt [40 | exp i K [(x —£) cos @ + y sin O} 
S —x 0 


x Qexp[—K@ +h] Gof [Kh +0] + ROoj [KE +H] + %Sin[Ke+m]aK. (16) 


Partikulare Lésungen der Differentialgleichung 


0 2 2 2 
FSET ip a0 mit We eee 


Ox c Oy? 
sind 
fi = exp [i K (x cos O 4+ y sin O)] Coj [m (z + h)], 
fz = exp [i K (x cos O + y sin O)] Gin [m (z + h)] 
mit 


ic KcosO 
, : 


m2 — K2 | 


1 J. K. Lunde, Trans. Soc. Nav. Arch. Mar. Eng., 59 (1951) S25—76. 
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Fiir die Geschwindigkeitskomponenten u, v, w erhalten wir jetzt unter Beriicksichtigung der 
Kontinuitatsbedingung fiir z—¢< 0 und z+2h+6> 0 den Ansatz 


aoe ‘dO { exp { K [le —£) o0 sin 
=a | fot) ae at [40 | p{i K [(x—) cos@ + y sin@}} 


x {—i K cos @ [exp [K (z—)] + exp[—K (2 +2h + 0)] + F,Gof [K ( +h] 
+ Fy Gin [K ( + h)]] + Fy Gof [m (@ + h)] + Fy Sin [me + Ay} aK, 


vay | faeeratat [40 fom [i K{(@—£) cos @ + y sin @}] 


x {—i K sin 0 {exp [K (« —£)] + exp[-— K & + 244 0)] + F, Gof [K (« + h)] (17) 
+ F, Gin [K (z + h)]} + F, Cof [m (z + h)] + Fy Gin [m (2 + h)]} dK , 


wy, | [ol60) dé dt [40 [exp [i K {(« —&) cosO + y sin O}] 
S —2 0 
x {K {— exp [K (2 —2)] + exp[—K (z+ 2h + ¢)] —F, Gin [K (z + A)] 
— F, Gof [K (« + h)]} = ((F, cos © + F, sin @} Gin [m ( + h)] 


+ {F, cos O + F, sin 9} Coj [m (z + h)])} dk. 


Mittels der Randbedingungen (15) am Flachwasserboden kénnen wir drei der im Ansatz (17) 
_ auftretenden Funktionen F,, F,, F;, F,, F;, Fg eliminieren, und es ergibt sich fiir z —¢ > 0 und 
 z+2h+¢>0 fir die Geschwindigkeitskomponenten u, v, w der Lésungsansatz 


i =a | foe ted [40 fom [i K {(x —£) cos 9 + y sin O}] 


x {21 K cos O {exp [— K (h + £)] Gof [m (& + h)] — Gof [Kh +9] 
x exp [— K (z + h)]} + i K cos O {— oj [K (z + h)] + Coj [m ( + h)]} 
x F* —i K cos O Gin [K (z + h)] F? + Gin [m (z + h)] FF} dK, 


5 =a | [eo tsa [a0 fom [i K {(« —£) cos @ + y sin OF] 


x {2i K sin O {exp [— K (h + £)] Gof [m (z + h)] —Gof [K (h + 0)] 


x exp [— K (@ + h)]} +i K sin O {— Gof [K (z + h)] + Gof [m@ + A)]} (17’) 
x F¥ + {—iK sin Gin [K (@ + h)] + =" Gin [m@ + h)]} 


cos O 


x ES Gin © 


Sin [m (z+ h)] Fo} dk, 


he = aa | [eo ted [a0 fom [i K {(x —£) cos O + y sin O}] 


x {2 K (Gof [K (hb + O)] exp [— K (¢ +h)] +S exp [—K (h +0] 


x Gin [m @ + A)]} + K{— Sin [K& + A] += Gin [m (2 + h)]} 


x F* + K {—Gof [K (z + h)] + Gof [m (« + h)]} Fr} dK. 
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Durch die drei Randbedingungen (14) an der Wasseroberflache sind die hier auftretenden 
Funktionen F*, F3, Ff bestimmt. Damit kennen wir alle interessierenden GréBen. Nach (12) er- 
halten wir die Gleichung der freien Oberflache zu 


(x,y) =— | | o(€,£) dé dt [ 10 [ex [i K {(x —£) cos +-y sin O}] 


X [Co} [K (h 4D] 'm {20 K Gof [K h] — (2 K + ic cos @) Go} [mh]} 
_ Gin [K (h +O) K{Q2ym Gin [Kh] —QvK +iccosO)Gin[mA]}]=aK, (18) 
mit 
A =4yK?m{29K +iccosO} + g{K oj [K h] Gin [mh] —m Gin [K h] Co} [m h]} 
+ K{2vK +iccos@}2{K Gin [K h] Gin [m h] — m Gof [K h] Co} [m h]} 
— 49? K2 m{ K ©oj [K h] Coj [m h] — m Gin [K h] Sin [mh]} . (19) 
Der Druckwiderstand ist hier nach der Formel 


R =—2f [ple 0,2) 5 de de 
s 


zu berechnen. 
Beachtet man, daB nach (4) 
p(x, 0, 2) = ahs = aa| {oe de dt | a0 | exp [i K (« —&) cos O]ic K cosO 
Ss —7 0 


x 


y=0 


x {exp [— K (@ + £)] + exp [— K @ + 2h +£)] + FY of [K @ + h)] + FF Gin [K & + h)]} dk 
(20) 


ist, und daB man o(x, z) = =e hat, wenn y = y(x, z) die Gleichung der Schiffsoberflache ist, 


so erhalt man fiir den Druckwiderstand schlieBlich 
ba ce) 
R =—4eIm [eos dO | {m4 (U2 + V2) 2K BUUI+VI)+KC(P4+P)} 24K, 
—x 0 


21 
wobei wir folgende Abkiirzungen eingefiihrt haben: 
A = g pj [K h] oj [mh] + K {27 K + ic cos O} Gin [K h] Go} [m h] “ 

— 4? K? m Gof [K h] Gin [m h] , 
B = 29Km{2vK +iccosQ} + g Wo} [K h] Gin [mh] 

+ K {2K +iccos@} Gin [K h] Gin [m K] — 47? K? m oj [K h] Cof [m h] , 
C = g Gin [Kh] Gin [mh] + K{2y K + ic cos O}? Coj [K h] Gin [mh] 

— 492 K? m Gin [K h] Co} [m h], 

A = 4y K?m{2vK + iccosQ} +g {K oj [K h] Gin [mh] — m Gin [K h] Gof [m h]} 
4K {2K +iccos@}{K Gin [K h] Gin [m h] — m Gof [K h] oj [m h]} 
—4y? K?m{K oj [K h] Gof [m h] — m Gin [K h] Gin [m h]} , 

U +i V = [f Gof [K (t+ 2)] exp [i K cos 8 x] a(x, 2) de de, 


I+iJ =ff Gin [K (h + 2)] exp [i K cos @ x] o(x, z) dx dz. 
s 


(22) 


b) Hognersches Drucksystem, dreidimensionaler Fall. Das auf die Wasseroberflache 
wirkende Drucksystem sei dargestellt in der Form 


Teg) a T(E, ) dS dy | | exp [i {k (E—x) +1 (n—y)}] dk dl. (23) 
Nees 


= 
is 
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Fir das Potential © wird man auf den Ansatz 
1s ” E 
Dex.952) = gal | Tn) dé dn| [exp fiCkE—») +1 —y)} 
Ss = 


x F, Gof [m (z + h)] + F, Gin [m ( + h)]} dk dl (24) 
gefiihrt, wo m? = k? + [? ist. Partikulare Lésungen der Differentialgleichung 


ooelae ay af @ 
Oe ee er Og ea a2 


sind 
‘ fi = exp[—i{kx+ly}]GojIn(e+h)], fo =exp[—i{kx +ly}] Gin [n(@ +h)] 
2 Sit ee 


Fir die Geschwindigkeitskomponenten u, v, w erhalt man unter Beriicksichtigung der Kontinuitats- 


_ bedingung den Ansatz 


1 4 ; ; 
mal {MEn aan | [ow [i {h(E — x) + Uy —y)}] {ib (F, Cof fm (@ + B)] 

+ F, Sin [m (z + h)]) + F; Cof [n (e + h)] + F, Gin [n (z + h)]} dk dl, 
v= yaa | Tn) aE dy| [exp [i(k (E— x) +1) —y) J] {i 1(F, Cof [m (e + hy] 

are —o 

+ F, Sin [m (z + h)]) + F, Gof [n &@ + A)] + Fy Gin [n (z + h)]} dk dl, (25) 
(a TDG) aban (it exp (eke ox = 

i | (A) ean] | pli{k(E—*) +1) —y)}] 

[=m (FF Sin im + HY) + Fe Cof fm (@ + HN) + (FB + He) Sin Im & +H] 

+ (FF tS Fo [n (© +h) dk dl. 


Wieder kénnen wir mittels der Randbedingungen (15) am Flachwasserboden drei der im An- 
satz (25) auftretenden Funktionen F,, F,, F;, F',, F's, Ff eliminieren, und es ergibt sich fur die 
Geschwindigkeitskomponenten u, v, w der Lésungsansatz 


w= zal [ Mn adn exp FRE —%) +1 yy} {ik Gof Im +H] 
ed pe 


— Gof [n ( + h)]) F* + ik Sin [m (@ + h)] FX + Gin [n (@ + h)] F*} dk dl, 


va ae [Mem aea | fon [i(k (E—x) +1 —y)}] 


x | 1 (of [m (z + h)] —Goj [n (@ + h)]) F* ie 
+ il(Gin [m (z + h)] — 22 Gin [n (z + h)]) FE —+ Sin [n ( + h)] Fs deal, 


(ce) 


ay Mee i —x = 
w= ae] [72 a | p [i{k (E—2) +1 (y —9)}] 


x {rm (— Gin [im ( +] + —_ 


Gin [n (z + h)]) F? 


+ m (— Gof [m (z + h)] + Cof [n (z + h)]) FF) dk dl. 
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Durch die drei Randbedingungen (14) an der freien Oberflache sind die in diesem Ansatz auftreten- 
den Funktionen F*, F*, F* bestimmt. Die interessierenden GréBen sind damit samtlich bekannt. 
Nach (12) erhalten wir die Gleichung der freien Wasseroberflache zu 


c(e, 9) =A, | [tenaean| [ow [(k(E—x) +1 —y)}] 
S —o 


X (m Coj [m h] Gin [n h] — n Sin [m h] Cof [n h]) F dk dl (26) 
mit ee 


m—h+ PP, ee == he - 


und der Abkiirzung : 

A =4y m2? n{2ym? —ike}+{2ym?— i kc}? {m Gin [mh] Gin [n h] — n Gof [mh] Coj [n h]} 
— 4% m3 n {m Coj [m h] Cof [n h] — n Gin [m h] Sin [n h]} + gm {m Co} [m h] Sin [n h] 
—n Gin [m h] Co} [n h]} . (26’) 


Der Druckwiderstand ist nach der Formel 
R =—| [Tn 5 ay 
Ss 
zu berechnen, wo C(x, y) aus (26) zu entnehmen ist. Damit wird 


1 (pion 
R=—, Im [fe V2) (m Gof [m h] Gin [nh] 


— n Gin [m h] Cof [nh] =” dk dl, (27) 
wenn wieder die Abkirzungen 
U+iV=ff T(x, y) exp [i (kx + Ly)] dx dy (28) 
Ss 


eingefiihrt werden. 


c) Hognersches Drucksystem, ebener Fall. Hier ist v =0 und die Abhangigkeit von y 
entfallt. Ist das auf die Wasseroberflache wirkende Drucksystem dargestellt in der Form 


T(x) =5— i T(é) dé | exp [i k (x —£)] dh, 


so wird man fiir das Potential ® auf den Ansatz 


D(x, 2) = za | TE) a | exp [i k (x —2)] { F, ©of [k (= + h)] + F, Gin [k (z + h)]} dk 
Ss —a 


gefiihrt. Partikulare Lésungen der Differentialgleichung 


any ne: _ af ay 
sind f 
f= exp (ikx)Cof[m(e+h)], fy —exp (ix) Gin [m(e + f)] 
mit m? = k? + soe Fiir die Geschwindigkeitskomponenten u und w machen wir unter Beriick- 


sichtigung der Kontinuitaétsbedingung den Ansatz 
uy | TE) as [exp [ik —8)] {ik (Gof + M)] + F Sin [h ( +m) 
s —o 
+ F, oj [m (z + h)] + F, Gin [m (z + h)]} dk, 
co (29) 
- =r,| 7) as | exp [i k (w —€)] {— k (F, Gin [A (@ + A)] + Fy Gof [ke + h))) 
s —o 


a 22 Sin [m (z + h)]) + Fy Cof [m (z + h)]} dk. 
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Mittels der Randbedingungen u = 0 und w=0 am Flachwasserboden z = —h kénnen zwei 
der im Ansatz (29) auftretenden Funktionen F,, F,, F,, F', eliminiert werden, und wir erhalten 
fiir die Geschwindigkeitskomponenten u, w den Lésungsansatz 


uaz | TE) a [exp [ik (© — 8) (6k (— Gof Th ( + )] + Cof [m & + B))) FY 
+ i(—k Gin [& (e +h] + m Gin [m (+ h))) FE} ak, 


aes (29") 
oe a T(é) as | exp [i k (x —é)]\k(— Gin [hk (z + h)] + = Gin [m (z + h)]) F} 
Ss —oo : 


+h (—Gof fh @ + AY] + Gof fm + H]) Fe ae 
Durch die beiden Oberflachenbedingungen 


Ou . Ow 
Bin tae 


2 Ow rae | 
ee ar toe Pe 


die 77 == (0 


sind die Funktionen F* und F} bestimmt. Damit sind alle interessierenden Gréfen bekannt. Die 
Gleichung der Wasseroberflache ergibt sich nach (12) zu 


ee) = 54, [Te ae[ew ike pp oy OTM aie Ml Cele, (30) 
mit ve 
A =4ivk?m(c—2ivk) +k (ce —2ivk)? {m Gof [kh] Coj [mh] — k Sin [k h] Sin [m h]} 
— g {m Gin [k h] Coj [m h] — k Gof [k h] Sin [m h]} — 4 v? k? m {k oj [k A] Co} [m h] 
—m Gin [k h] Gin [m h]} (31) 
und 
m? = he 4 2% 


Der Druckwiderstand ist nun nach der Formel 


dt) 
Re i T(x) & de 
s 
zu berechnen, wo ¢(x) aus (30) zu entnehmen ist. Damit wird also 
1 oa 7 k 
yon | [UP + V2] (m Gi [HH] Gof [mH] — Gof fH] Sin fm hl) Fak (82) 


U +i Vi= fi T(x) exp [tk x] dx (33) 
5 


gesetzt ist. , ; 
Wir bemerken noch, daB die von Srétensky fiir unendlich tiefes Wasser aufgestellten Formeln 
aus den hier angegebenen beim Grenziibergang h —> co hervorgehen. 


(Eingegangen am 12. November 1958.) 
Anschrift des Verfassers: Dr. F'. Kolberg, Wiirselen b. Aachen, Kaiserstr. 69. 
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Theory of flow through profiles with large camber and thickness 
arranged in cascade with small pitch chord ratio 


By R. Oba 


1. Introduction. With speeding up and the increased power of axial flow hydraulic machinaries, 
the blade elements in these impeller and guide vane etc. tend to consist of the profiles with not 
small camber and thickness, and also the pitch chord ratio of these straight cascades tend to di- 
minish. Especially, it is the case with gas turbines. For the theoretical studies on the flow through 
such a straight cascade the reports have been published by W. H. Isay', M. Shirakura? and 
H. Murai3, etc.*. In the first two investigators’ solutions, however, the numerical calculations 
are very complicated and H. Murai’s solution seems to be not quite satisfactory for us, either. 

In the former theories® of the cascade that consists of the profiles with not so large camber 
and thickness, even if the pitch chord ratio is as small as to be used in the Kaplan turbine, there 
are difficulties® for the poor convergency of the infinite series to be used and so on. It is well known 
that the theories® 7 & % 10,11 with assumption of large pitch chord ratio can not be applied to such 
a cascade??. 

Accordingly, in this paper, the author proposed an analytical method for the solution on the 
flow through the straight cascade with small pitch chord ratio, in which the numerical calculations 
are comparatively easy not only in the case of the profiles with large camber and thickness as 
gas turbines, but, in those with small camber and thickness as axial flow turbines, -pumps and 
compressors. 

Now, depending upon the method in common use, to solve the potential flow through the 
straight cascade is reduced to transforming the appointed cascade into a circle or two or more. 
In these cases, if we transform the cascade profiles into a circle to simplify the calculation, there 
arises a difficult point that the accuracy of the calculation becomes deteriorated, as M. Shirakura 
pointed it out, because of the midway points of the two neighbouring profiles are transformed 
into the point at infinity, and because of the corresponding points of the ones at infinity in the 
z-plane draw near a circle in the ¢-plane. H. Murai has avoided it skilfully. 

In this paper, to more simplify the calculation, the appointed cascade is transformed into the 
one of the profiles with small camber by the Howell’s transformation and by the inverse trans- 
formation, and so it becomes possible to transform the cascade profiles into the figure most re- 
sembling a circle by the next transformation; thus we can remove the above mentioned difficulties. 
Through such analysis, the numerical calculations become very easy for various cascade arren- 
gements. And also it is ascertained that a part of the former theories?, is the special case of this one. 

For the numerical examples, in the gas turbine cascade used by F. Numachi and T. Kurokawa14 
and by H. Hausenblas’ in their experiments, as H. Murai has calculated, the surface pressure dis- 


1 W. H. Isay, Z. angew. Math. Mech. 33 (1953) S. 397. 

 M. Shirakura, Trans. Jap. Soc. Mech. Engrs. 18 (1952) p. 16; M. Shirakura, Trans. Jap. Soc. Mech. 
Engrs. 21 (1955) p. 577. 

° H. Murai, Z. angew. Math. Mech., 35 (1955) S. 48 and Rep. Inst. High Sp. Mech., Tohoku Univ. 7 (1956) p. 15. 

4G. S. Samoilovich, Prikladnaya Matematica i Mekhanica 14 (1950); C. H. Wu and C. A. Brown, J. Aeron. 
Sci. 19 (1952) p. 183. 

° E. Kasahara, Trans. Jap. Soc. Mech. Engrs. 20 (1954) p. 577; T. Kawasaki, J. Jap. Soc. App. Mech. 5 
(1952) p. 3; H. Kikuchi, Rep. Inst. High Sp. Mech. Tohoku Univ. 6 (1956) p. 103. 
fesse eG Trans. Jap. Soc. Mech. Engrs. 17 (1951) p. 57 and Mem. Inst. High Sp. Mech., Tohoku Univ. 2 

p. 93. 
? G. Kamimoto, Trans. Jap. Soc. Mech. Engrs. Vol. 17 (1951) p. 58. 


8 H. Murai, Trans. Jap. Soc. Mech. Engrs. 17 (1951) p.'52 and Mem. Inst. High Sp. Mech. Tohoku Univ. 2 ; 


(1950) p. 83. 

° H. Kikuchi, Mem. Inst. High Sp. Mech. Tohoku Univ. 10 (1953) p. 85. 

10 B. Fujimoto and K. Hirose, Mem. Faculty Engineering Kyoto Univ. 12 (1950) p. 20. 

1 R. A. Howell, Phil. Mag. 39 (1948) p. 913. 

12 In the analysis for the cascade that consists of the profiles with small camber and thickness, the numerical 
calculations based upon the above mentioned theories from footnotes } 2:3 are not always easy. = 

13 See footnotes 6, 7, 11. 

14 F, Numachi and T. Kurokowa, Trans. Jap. Soc. Mech. Engrs. 15 (1950) p. é . Hi 
Mech. Tohoku Univ. 5 (1951) p. 39. : P Darblinense cs, ee 

15 H, Hausenblas, Ing.-Arch., 19 (1951) S. 75. 
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tributions and other factors are calculated, and the good agreement between the calculated value 


and the experimented one is proved. In these calculations, the simplicity of the numerical cal- 
culation is also confirmed. 


2. Symbols. The following symbols are used in this paper: 


t/1 pitch chord ratio, 

wt 
y = mere stagger angle, 
) cascade angle, 


V_,%—_ velocity of flow before passing profiles in a cascade, angle between velocity V. 


and profile chord, ‘s 
Vins Cite velocity of flow after passing profiles in a cascade, angle between velocity V,, and 
profile chord, 
Va geometrical mean of velocities V_., and Vx, incidence angle, 

a circulation around a single profile, 

Ca lift coefficient of a single profile, 

Op zero lift angle, 

v surface velocity along the profile, 

W complex potential, 

z =x-+iy physical plane, 

: = Q exp lg, 

Cn = Oy eXpt Py =fy +1 Hy mapping function transformed the profile into a unit circle. 


The suffix N represents the number of the mapping plane, the suffices — oo and oo respectively 
signify infinite front and infinite back in the z-plane, and represent the corresponding points in 
the mapping plane. Besides, not indicating the profile chord expressly, we can take it following 
our own direction in our own way. 


3. Conformal mappings and mapping functions. As shown in Fig. 1, the cascade at issue is an 
arbitrary cascade? and consists of profiles with large camber, thickness and has a small pitch chord 
ratio. Such a cascade is transformed into a single figure in the ¢,-plane by the well known mapping 
function? 
ei? 


(4 == tanh fe z— | where z,, complex constant . (1) 
Then the single figure in the ¢,-plane is transformed into the one in the ¢,-plane in Fig. 2 by 
the transformation of coordinate, equation (2). Besides, we transformed it inversely into the ones in 
the ¢,-plane by the similar mapping function of equation (1), equation (2). 
6, = D,¢, + Dy where Dy, D, = k, e'%, complex constant , (2) 
: ¢; = tanh‘ ¢,. (3) 
In these transformations, determining D,, D, suitably, we can transform the cascade in the 
¢,-plane into the one that consists of profiles with small camber and has not so small a pitch chrod 
ratio. Besides, in an accelerated cascade as shown in Fig. 3, we must see to that the upper front 
and the lower back parts of the profile approach the profile chord to be determined in the following 
equation (4), and the points ¢,.,, ¢;. approach thesaid parts of the profile and keep away from 
the profile as far as possible. In a retarded cascade, it is just opposite to the above. 
Then we transform the mapped cascade of pitch chord ratio t,, of stagger angle y, into the 
single figure in the ¢,-plane by the mapping functions (3): 


6, = Dz, DP. where D,, D, = kz e'%, complex constant, (4) 
any of ty i”g fs + ii, g5 1 ob m Ge j 
be =a (#7log B= ™ + em log eae (5) 


The leading edges of the profiles in the C,-plane get on with the mid-points v of the radius of 
curvature. And, D,, D, must be determined, as one mid-chord point agrees with the origin of 
the ¢,-plane, as the eotile chord gets on the &,-axis, and as the chord length is unity. 


1 If the camber and thickness are not so large, and if the pitch chord ratio is not so small, we can calculate 
it by the R. A. Howell’s theory. 
2 See footnote 2 and 11 of page 276. 
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Well, even if the gas turbine cascade, the profiles in the z-plane have customarily a thin thick- 
ness on the camber line like an aeroplane wing. Then, in the calculation of the corresponding 


Fig. 2. f,-plane. Fig. 5. €,-plane. 


Fig. 6. ¢-plane. 


Fig. 3. ¢,-plane. 


Fig. 1...6. Physical plane and mapping planes. 


| 
points by the transformation in equation (4), we can resort to a thin profile theory like S. Abe’s! and 
G. Kamimoto’s*. Then, in the first approximation, even if pitch chord ratio be small, adopting the 
next relations, equations (6,), (6,), we can get the good approximation except the extreme neigh- 


1 See footnote 6 of page 276. 
* See footnote 7 of page 276. 
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bourhood of the leading edges: 


Se: 2m cos g; ; 2 m sin 9; . 

&= = (cos y, tanh—! ee + sin y, tan ea + ty sin yy, (61) 
F 2mst, . ; 

ns =~ Sin (Ps — Ps) log 05 + n ty cos yy (62) 


where 


s=/1 +m! + 2m? cos2y,, 
f= 1+ mt— 2m? cos29;, 
ee ee eee eee 


Further, if we denote the radius of curvature of the leading edges 0, the first approximate points 
of the nearby leading edges are calculated by the relation 


, x 6 (1 — m4)? 
(i, — yy = SE) (7) 
2t,m Ys? 
where parameters m, ~;7 must be determinated by the two relations 
1— rm? 
CE TS Oe eye ele 2 (8) 
1 2 2m cos p : 2m sin p 
olen ies —s ST col 5T 8 
— (cos vq tanh ron + sin y, tan ae : (84) 


For given t,, y4, these equations (8,), (8) are solved, and the results are shown in Fig. 7. 


™ 
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Fig. 7. Pitch chord ratio t,, stagger angle y, and cascade parameter m. 


To calculate more accurately the points in the ¢,-plane, we put the first approximate values 93, 
s into the next relations, equations (9,), (9,) with which the equation (5) was solved, and calculate 
the corresponding values &4, 74 of €,, 7, and so calculate the differences A&,, An, from the appointed 
values &,, 14: 


oo a [(T, + T,) cosy, + (T; + Ty) sin y,4] + 1 ty sin yy, (91) 
m = 5 (T, — T,) sin y, + (T; — T,) cos y4] + nt, cos y 4 (92) 
where 5 2 oem cost, 
kes Qs ™ COS | ae oD ee 
aE rary sian Le REST pm of * 
_12Q,msing er ay ies HY 
i, = tan Ser ae P= tan 02 — m2 
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Then the corrections Ao,, Ap; of 055 Ps are calculated by the relations 


tT ©. 
Aes = A (Nia bs— Nea a) » (10,) 
, 1 
03 Ags = A (Nias + Noa §4) (10,) 
where d 2 atte : : 
sae 2 A62 Oy Seeks fs i , 
N, = mis cos Va cos Qs ( : oy Gs ae gs je Hs sin Va sin Ps (" Fe Gs = en - as ) ’ 
j — me p72 oie we ee. x Ail Pde} /9 2 
N, = = sin 74 COS Ps ( ! fa Qs &s a wd cos 4 sin #5 a Qs =e Qs fe x } : 
f' = (08 + my? — (2 m5 cos.y5)?, ff” = (1 + m? 93”)? — (2 mg; cos g;) , 
Age ING ING. 


In the same way, we can further proceed to the successive approximate calculations. But, mostly 
the second approximate values are sufficient for our use except the points of the nearby leading edges. 

The above transformations from equation (1) to equation (5) are essential, in such a case of the 
arbitrary profiles with large camber, thickness and small pitch chord ratio. If the profiles have 
not so large.a camber, thickness and not so small a pitch chord ratio or if the absolute value of the 
stagger angle is small, we can proceed easily to the calculation neglecting the transformation, 
equations (1), (2), (3)', or equations (2), (3), (4), (5)? respectively. But, even if the camber and the 
thickness are small, if the pitch chord ratio is too small, the transformation of the equation (5) 
with the transformation, equations (1), (2), (3) will be made more easily. 

If the camber and the thickness are not so large, and if the pitch chord ratio is not so small, 
the required figure in the ¢;-plane is quite a circle. But, if the cascade is like the one of degree as 
shown in Fig. 1, we can not avoid the figure being slightly unlike a circle. 

Therefore, generally, the figure in the ¢,-plane must be transformed into a figure to make it 
resemble a circle closely by the mapping functions 


Cs = D6, + D, where D,, D; = ks e'%*, complex constant, (11) 
7 7 ; a E a : J where T, real constant. (12) 
? 6 


If the figure in the ¢,-plane is not quite a resemblance to a circle, we can transform it again into 
the closed curve bearing a striking resemblance to a circle by reiterating the above two trans- 
formations. 

Finally, by the well known T. Theodorsen — I. E. Garrick’s transformation®, the closed 
curve bearing a striking resemblance to a circle is transformed into a unit circle, whose center is 
on the origin in the ¢-plane: 


Cs = D,¢, + Dg where Dy, D, = k, e'%, complex constant , (13) 
co C. ; 
Ge == Gexp SS a where C, = A, +i B,, A,, B,, real constant . (14) 
n=0 


where D,, D, must be determined as the closed curve in the ¢,-plane border on sufficiently a unit 
circle whose center is on the origin. 


4. Complex potential and the hydrodynamic characteristics of the profile in a cascade. Generally, 
the complex potential corresponding to the z-plane and taking the unit circle for one of the stream 
line in the ¢-plane is given by the formula 

NS emcees G1. iD_ Gees iD cia 
W =, (05 + Io le Slog tat = oe ea ste (15 
PE WS delaras Be oak SEES 2a 2 eS a tir a ie bo 
where M is the strength of sink and source, J". and J’, are the circulation respectively. On the 
other hand, because of the circulation in the unit circle is equal to the one of the profile, next 
condition is to be satisfied : 


Ps ere. (16) 


* In such a case, the theories of footnotes 6 and 7 of page 276 are special cases of this one. 
* In such a case, the theory of footnote 11 of page 276 is a special case of this one. 


3 T. Theodorsen and T. E.Garrick, General Potential Th f Arbit i i 
Peosne mention) eory 0 : rary Wing Scetion, N. A.C. A. 
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Besides, as the flow conditions at infinity, the next relations must be satisfied: 


dW 


oP eel ge ety (17) 


z=+0 


In the results, the necessary ee potential is ine as follows: 


WBE (een ng tee a0) tg MESH) HE gg CEE 
za a Pees) =e es cai 
This agrees with the one given . H. Murai for aes s sake. 
And the circulation J” being left as undetermined constant is determinated by the Kutta-Jou- 
kowsky’s hypothesis. If the allies edge rounds out and is not pointed, we calculate I’ as dW/d¢ =0 
at py = @y. Then the circulation J" is determinated as follows: 


1b! _ K_~sm@+y+%*)—K, ae 
tae K_,. K+ K,, cos x, + K,, cos x_ 


(19) 
where 

K-« exp t#-0 =C-0 —Cp, Ko exp 1 H%0 = Co —fbp. 
Then the hydrodynamic characteristics of the profiles in a cascade can be calculated Eee MS oi: 
That is, the lift coefficient is determined as follows: 


Die : 
C= y= I sin (a —a,)- (20) 
The inclination of the lift coefficient curve I and the zero lift angle a, in this equation are cal- 
~culated as follows: 


2 ko K —— 
EBs) t \K2 + K2—2 K,, K_. COS (Xo, ~ Xoo) (21) 
l Kes Koo Ko 5, cos #5, 44k, cose 
K_., sinz,, — K,, sin x_ 
eS po th ED ee 22 
On y van Ko = cosuin = S cos | (22) 


Then the velocity in the physical plane, z-plane, is given as follows: 
dW dWd 


Especially, calculating the velocity distribution on the surface of the profile, it is 
7 = 1-02 Qs- Q- 5 (24) 
where 
2 [F, sin (« + y) + F, cos (« + y)] + Ti 
Q, = F, 7”) 
ae Os co ee | Ooo — may 
Cee eee re? be) a (P — Yo) 
ie : 2 Ll 
Fy = 2 sin (p_,, — P) — == sin (p— 9.) — EB sin (P—9_.)» 
1 gt meal eT 
Doi Cate Om Ae ie cos (p — 9.) —* cos (P — P_ wo) » 
2 A 1 
Fy = [E222 —2 008 (p —9..)| |S — 2 005 —9)] 
, ee 
oo 2 00 2 
Uae male — J} (nA, cosng +n B,sinng)| + [2 4, sinngy—n B, cos ng) ; 
Qs n=1 n=1 
1 

1+ ef —2 of cos2y,]2 

Oaks 1+ ef —2 2 cos2g,| ” 
77 aes 
0,= mt ky Vf" [1 + 05 —2 4 cos 2 (~;— Ysr)] ? , 
20s Go ty 


One 1+ of —2 of cos 2 q | 
ore 1+ of — 2 93 cos 2 gy y 
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5. The limit of the application of this theory. The limit of the application of such a theory depend- 
ing upon the required accuray and the trouble of the calculation, in short, to decide this are the 
degree of that the points at infinity in the z-plane draw near the unit circle in the ¢-plane, the con- 
vergency of the Fourier series to be used and the accuracy. The corresponding of the former is m1 
as shown in Fig. 7. By such m, depending upon the required accuracy or the camber, thickness 
of the profile, it is proper roughly that m = 0.70 is the limit of the application for most of the 
former theories. Once over this limit, abruptly, either the accuracy of the calculation drops, or 
the convergency of the Fourier series detoriorates. In the author’s theory, such difficulties are 
removed with less compricate numerical calculations. 

The more m draws near to unity, the more the points in the ¢-plane corresponding to the ones 
to take at regular intervals on the profiles in the z-plane are in close order. On the other hand, 
for the accuracy of the calculation, at least, the Fourier series needs to be calculated precisely near 
the close points. In such conditions, because of the good convergency of the Fourier series in the 
author’s theory, and because of the points to be in close order, it becomes easily to keep good ac- 
curacy by means of a few corrections for the Fourier series by the ordinal simple method. There- 
fore, the limit of the application of the author’s theory must be extended considerably. 

In the rear numerical examples, even if the profiles are built up with large camber and thickness, 
even if m is 0.97, correcting the Fourier constants from the method of summation (36 points method) 
as we above mentioned, the author can calculate them with good accuracy. 


6. Examples of the numerical calculation. 


a) The experimented profile by F. Numachi and T. Kurokawa. For the cascade of 
guide vane blades of gas turbines experimented by F’. Numachi and T. Kurokawa, that the profile 
form is shown in Fig. 8 and is tableted by Table 1, the lift coefficient and the surface pressure 
distribution on the profile in the case? of y = — 45.61°, « = 45.61° are calculated. 


Fig. 8. The profile experimented by F. Numachi — T. Kurokawa. 


Table 1. The dimension of the profile from F. Numachi and T. Kurokawa 


x/I 0 0.01722 0.04990 0.07485 0.09980 0.14970 0.19960 
Yoll 0.01871 0.06257 0.10025 0.12525 | 0.14691 | 0.18268 0.21088 
yyll 0.01871 0 0.01826 = 0.03373 0.04795 0.04241 0.09266 
x/l 0.24950 0.29940 | 0.34930 | 0.39920 0.44910 0.49900 0.54890 
Yoll 0.23663 0.25584 0.26796 | 0.27275 0.27031 0.25983 0.24266 
Yyll 0.10918 0.12255 0.13249 0.13912 | 0.14222 0.14276 0.14037 
x/l 0.59880 0.64870 0.69860 0.79840 0.89820 0.96307 1.00000 
Yoll 0.22250 7.19855 0.17260 | 0.12031 0.06806 0.03363 0 

Yyll 0.13418 0.12445 0.11158 | 0.08109 0.04271 0.01497 0 


In the analysis, for the polygonal line being on the upper surface near the trailing edge of the 
profile as shown in Fig. 8, the author uses the nonpolygonal line to round out slightly in the same 
way as H. Murai did. As the value of the calculated lift coefficient by this theory is 2.08, it differs 
from the above mentioned experimented value 2.02 by 3°% and from the calculated value by 
H. Murai 2.10 by —1%. The surface pressure distribution is calculated by the pressure coefficient 


1 Even if the case with the transformations, equations (1), (2), (3 
the standard roughly. = eke) 


® See footnotes 5, 6, 7, 8, 9, 10 of page 276. 
° H. Murai had taken up this for the numerical example of his theory. 


(4); being calculated by t/l and y, m is 
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1 : 
(P — P_,) ee o V2. referring to the velocity at front infinity V_. as the report by F. Numachi 


and T. Kurokawa shows. The result is shown in Fig. 9 to compare it with the experimented value, 
and with the calculated value by H. Murai. The full line shows the author’s calculated value 


and the broken one shows the H. Murai’s calculated one, then, the two agree with each other 
fairly well. Besides, they seem 


to agree with the experimented 
value considerably. 


b) The experimented 
profile by H. Hausenblas. 
For the cascade B,,! among 
the cascade of blades of 
gas turbines experimented by 
Hi. Hausenblas, that the profile 
form and the cascade arrange- 
ment are shown in Fig. 10, 


the lift coefficient curve and Author's cal. 
its surface pressure distribution aes Murais cal. 
on the profile are calculated. ° — Numachi-Kurokawa’s exp. 


First, it is shown in Fig. 11 | 
to compare with the calculated 
and the experimented lift co- 
efficient curve. In Fig. 11, the 
axis of abscissas takes « and 
t= f, — B_,.» and the axis of OP Re ae 50 = SP AGO 0 BR a 
ordinate takes C, I[t. In this Fig. 9. Comparison between the calculated surface pressure distribution and the experimented 
figure, the full line shows the one by F. Numachi — T. Kurokawa. 

author’s calculated value, the 
dotted and dash line shows the 
H. Murai’s calculated value 
and the broken line shows the 
experimented value, then the 
former two agree with each 
other fairly well, and also seem 
to agree with the third. Next, 
in Fig. 12, with the surface 
pressure distribution in the 
case B_., = 37° that the drag 
lift ratio is the minimum, the 
calculated one and the experi- 
mented one are compared. The 


fv | by | v |v 


axis of abscissas takes the pres- 
sure coefficient referring to the | 445° | 32.0° |-22.07°| 0.708 
velocity V and the axis of 
: : ft | aft | af | a,/t | oft 
he chord-wise | { 
ordinate takes the c nao 0078 


| 03/5 | 0.374 | O260: 


distance of the profile drawing 
in the lower side. In this figure, Fig. 10. The profile and the cascade arrangement experimented by H. Hausenblas. 
the author’s calculated one in 

full line and the H. Murai’s calculated one in dotted and dash line agree with each other fairly 
well and also they seem to agree with H. Hausenblas’ experimented one, too. 


Well, as the Fourier analysis in these numerical calculations (1), (2) was based on the ordinal 
method of summation (36 points method) altogether, the correction as above mentioned was 
made only in the nearby parts of £.,¢— w- In either case, the transformation of the equation (5) 
was tentatively calculated to the third approximation, and it was found that the second approxi- 
mation agrees with the exact one up to the fourth significant figures. Therefore, the author guesses 
it that the second approximation for the calculation will be practically sufficient in most cases. 


1 The same as the foregoing footnote. 
19* 
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7. Conclusion. 


1. A theoretical method is proposed, in which the flow through the straight cascade built up 
the profiles with large camber and thickness and with small pitch chord ratio, and especially the 
surface pressure distribution and the lift coeffi- 
cient curve can be calculated by the compara- 


70 . . . 
tively easy numerical calculation. 


2. By the numerical calculations for the 
_cascade of the blades of gas turbine experi- 
mented by F.Numachi and T. Kurokawa, 
and by H. Hausenblas, it is confirmed that 
the author’s calculated results of the lift 
coefficient and the surface pressure distri- 
bution agree with the experimented ones 
and the H. Murai’s calculated ones, and that 
the numerical calculations are easier than the 


H. Murai’s. 


Author's cal, In conclusion, the author expresses his 
—-— Murai's cal, cordial gratitude to Prof. F. Numachi, Fa- 
——— Hausenblas'exp. culty of Engineering, Tohoku University, for 
his generous guidance throughout this study, 
g« and atthe same time, to Assist. Prof. H. Murat 
’ -/0° 0° 70° 20° go°. of the Institute of High Speed Mechanics, 
: baa ES ser or er MPS ETF ar = Tohoku University, for his helpful suggestions 
: : o) and his offer of the numerous data in the 
Fig. 11. Comparison between the calculated lift coefficient curve and : é 
the experimented one by H. Hausenblas. numerical calculations. 


60 


40 


Author's cal. 
a —-— Murai’ cal. 


—— — Mausenblas'exp. 


Fig. 12. Comparison between the calculated surface pressure distribution and the experimented one by H. Hausenblas. 


(Eingegangen am 15, November 1958.) 
Anschrift des Verfassers: R. Oba, Sendai (Japan), Tohoku University, Institute of High Speed Mechanics. _| 
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ZUR INFORMATION 


Mehr als 50 Jahre sind vergangen, seit J. HELMSLEV seine ,,Neue Begriindung der ebenen Geometrie”’ ver- 
éffentlichte, die M. DEHN 1926 als den héchsten Punkt bezeichnete, den die moderne Mathematik iiber 
EUKLID hinausgehend in der Begriindung der Klementargeometrie erreicht hat. Seitdem ist das Operieren 
mit den Spiegelungen ausgebaut worden, die einerseits die geometrischen Objektive reprdsentieren und 
andererseits die Bewegungsgruppe erzeugen. Das Ziel des Buches ist, durch eine neue Darstellung der 
Probleme der ebenen metrischen Geometrie jedem Interessierten den Zugang zu der neueren Entwicklung 
zu erschlieBen. 

Die ebene metrische Geometrie wird rein gruppentheoretisch aus wenigen einfachen Gesetzen iiber Spiege- 
lungen entwickelt. Die Gesetze gelten sowohl in der euklidiscken als in den nichteuklidischen Geometrien; 
sie implizieren weder Anordenbarkeit noch freie Beweglichkeit. Geometrische Satze werden durch Rechnen 
mit Spiegelungen bewiesen. : 

Es wird der Zusammenhang der ebenen metrischen Geometrie mit der Theorie der projektiv-metrischen 
Ebenen und der dreidimensionalen metrischen Vektorraume hergestellt, deren Bewegungsgruppen gleich- 


falls durch Spiegelungsgesetze beschrieben werden. 
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Zweite, neubearbeitete Auflage. Mit 163 Abbildungen. XI, 209 Seiten Gr.-°. 1959. 
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